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EXECUTIVE  SUMMARY 

A  ship  service  fuel  cell  power  system  has  been  explored  in  a  multidisciplinary  and 
multidimensional  effort  to  improve  understanding  the  dynamic  behaviors  of  individual 
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components  and  the  performance  of  the  system.  Models  of  individual  components  of  the 
system  are  developed  in  the  Virtual  Test  Bed  format,  which  then  allows  easy  reconfigura¬ 
tion  of  the  system  and  allows  nearly  unlimited  studies  of  system  tradeoffs.  Developed 
models  include  those  for  the  fuel  cell  stack,  compressors  and  expanders,  methanol  re¬ 
former  (as  a  precursor  to  strategic  fuels  reformer),  hydride  bed  (as  an  alternative  hydro¬ 
gen  storage  technology)  and  control  systems.  These,  together  with  other  models  devel¬ 
oped  under  other  research  programs,  have  been  combined  into  a  system  simulation  t  at 
allows  to  study  the  first-order  system  response. 

At  the  component  level,  research  focused  on 

•  characterizing  and  modeling  the  dynamics  of  proton  exchange  membrane 
(PEM)  fuel  cell  stacks.  The  work  was  performed  on  both  the  standalone 
model  and  the  VTB  model. 

•  characterizing  and  modeling  the  dynamics  of  air  compressor  in  VTB 

•  characterizing  and  modeling  the  mathanol  reformer  in  VTB 

•  characterizing  and  modeling  metal  hydride  bed  in  both  standalone  model 
and  VTB  model 

•  Modeling  a  feedback  control  system  using  embedded  MatLab  object 

At  the  system  level,  a  representive  shipboard  fuel  cell  power  system  was  assem¬ 
bled  and  simulated  in  VTB  environment  for  study  of  the  system  transients. 

Specific  accomplishments  of  the  studies  included 

•  PEM  Fuel  Cells  Standalone  Model 

o  A  test  stand  was  assembled  for  the  purpose  of  characterising  the 
transient  response  of  PEM  fuel  cell  systems 

o  Half-cell  data  was  gathered  as  input  for  formulation  of  a  detailed 
fuel  cell  model. 

o  A  mathematical  model  for  the  effect  of  adding  an  air  bleed  to  a  re¬ 
formate  containing  monodioxide  in  the  anode  feed  was  developed. 
The  model  can  very  well  predict  the  oxidation  rate  of  CO  and  the 
results  are  consistent  with  existing  experimental  data. 

•  PEM  Fuel  Cells  VTB  Model 

o  The  system  level  model  was  developed  to  account  for  the  interac¬ 
tive  reversible  and  irreversible  potentials  with  dependencies  of 
pressure,  temperature,  fuel  and  air  consumptions,  and  heat  produc¬ 
tion. 

•  Compressor 

Several  system  level  models  of  different  level  of  details  were  de¬ 
veloped  for  different  research  foci. 

A  general  centrifigual  compressor  was  developed. 
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•  Hydrogen  Storage 

o  A  large  hydride  bed  was  instrumented  and  experimentally  charac¬ 
terised  during  charge  and  discharge  processes. 

o  A  detailed  model  of  the  hydride  bed  was  formulated. 

o  A  system-level  model  of  the  hydride  bed  was  developed. 

•  Mathanol  Reformer 

A  CSTR  model  of  the  methanol  reformer  was  developed. 

•  Feedback  Control  System 

A  MatLab/Simulink  object  for  the  feedback  control  system  was 
developed  based  on  the  PID  control  algorithm. 

•  Shipboard  Fuel  Cell  Power  Plant  System 

o  The  dynamic  system  is  assembled  in  VTB  computational  environ¬ 
ment.  The  system  primalrily  consists  of  a  PEM  fuel  cell  stack  that 
can  be  easily  reconfigured  in  terms  of  voltage  and  power  levels,  a 
steam  methanol  reformer  for  fuel  processing,  or  alternatively  a  hy¬ 
dride  bed  for  fuel  supply,  a  compressor  driving  air  to  the  fuel  cell. 
The  compressor  is  driven  by  a  motor,  which  uses  the  system  output 
power  once  it  starts. 
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INTRODUCTION 

This  report  summarizes  progress  in  modeling  of  shipboard  fuel  cell  power  plant 
project  conducted  at  the  University  of  South  Carolina  as  a  collaborative  effort  between 
the  Departments  of  Chemical  Engineering  and  Electrical  Engineering.  The  work  had 
three  significant  foci.  The  first  was  enhancing  our  knowledge  of  the  physical  perform¬ 
ance  of  individual  electrochemical  and  mechanical  power  components  through  the  devel¬ 
opment  of  detailed  mathematical  models.  The  second  was  integrating  the  knowledge  ob¬ 
tained  by  standalone  models  into  VTB  simulation  models  of  individual  components.  The 
third  was  to  prototype  a  representive  fuel  cell  power  system  in  VTB  environment  for 
study  of  dynamic  performance. 

A  significant  focus  of  the  work  was  on  integrating  the  knowledge  of  chemical  en¬ 
gineers  into  dynamic  simulation  models  that  could  be  used  to  study  the  electrical  per¬ 
formance  of  fuel  cell  power  systems.  This  report  refers  to  the  construction  of  various 
types  of  simulation  models.  First  are  stand-alone  models  that  are  generally  used  for  de¬ 
tailed  studies  of  the  internal  dynamics  of  a  particular  component,  often  including  detailed 
representations  of  the  physics,  up  to  three-dimensional  transport  and  reaction  processes. 
Second  are  system-level  models  which  generally  have  lower  complexity  since  they  are 
intended  to  represent  mainly  the  characteristics  observable  at  the  terminals  of  the  device. 
These  models  are  useful  for  the  study  of  systems  that  incorporate  these  devices.  The 
process  of  model  construction  generally  starts  with  physics-based  models  that  attempt  to 
fully  describe  experimental  data,  then  progresses  to  simplification  of  those  models  into 
forms  that  run  fast  enough  to  be  used  effectively  in  system  level  studies.  At  this  level, 
resolution  of  some  of  physical  effects  is  generally  lost. 

This  report  is  arranged  as  follows.  In  the  next  five  sections,  we  summarize  the 
model  development  of  major  components,  such  as  PEM  fuel  cells,  centrifigual  compres¬ 
sor,  methanol  reformer,  hydride  bed  and  feedback  control  system.  Then  a  representative 
fuel  cell  power  system  is  described.  Finally,  we  describe  the  project  deliverables  and  the 
publications  that  resulted  from  this  work.  The  appendices  contain  information  related  to 
developed  models  and  copies  of  the  publications,  and  the  student  education  information. 


FUEL  CELL  TECHNOLOGIES 

Fuel-Cell  Test  Station 

Experimental  infrastructure  and  procedures  were  developed  to  measure  the  physi¬ 
cal  parameters  and  verify  fuel-cell  performance  models.  A  “Dual  Channel  Fuel  Cell  Test 
Station”  was  purchased.  Each  channel  consists  of  a  single  mass  flow  controller  on  the 
cathode  channel  to  regulate  air/02  flow  and  three  mass  flow  controllers  on  the  anode 
channel  to  regulate  the  flow  of  H2/reformate/N2,  and  to  bleed  in  a  calculated  amount  of 
carbon  monoxide  (CO)  and  air.  A  four-way  valve  allows  connection  of  different  gases  to 
the  anode  and  cathode  main  streams.  Each  channel  also  contains  temperature  controllers 
for  the  anode  and  cathode  humidifier  bottles,  anode  and  cathode  gas  lines  from  the  hu¬ 
midifier  bottle  to  the  cell,  and  the  fuel  cell  itself.  The  fuel  cell  is  connected  to  the  anode 
and  cathode  gas  inlets  and  outlets  and  to  an  electronic  load  bank.  The  load  bank  allows  to 
electrically  load  the  cells  with  constant  current,  constant  voltage,  sweeping  voltage,  or 
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pulse  voltage.  All  the  mass  flow  and  temperature  controllers,  and  the  electronic  load  are 
computer  controlled. 


.  Characterize  dynamic  performance 


Experimental  Data 

Figure  1  shows  the  performance  of  a  typical  fuel  cell  obtained  at  70  °C  at  ambient 
pressure  by  scanning  the  voltage  of  the  fuel  cell  from  1.0  V  to  0.3v  at  a05  V  increments 
and  resting  30  s  at  each  voltage.  The  initial  drop  of  the  curve  (from  0.95  0.8  V)  results 

from  limitation  by  the  reaction  kinetics,  predominantly  the  oxygen  reduction  reaction 
(ORR).  The  linear  drop  between  0.75  and  0.4  is  controlled  by  ohmic  drop  in  the  cata  yst 

layer  and  the  membrane. 


Figure  1.  Performance  of  a  10  cm^  fuel  cell  using  H2  and  air  at  70  °C  at 
ambient  pressure 

This  curve  can  also  be  used  to  qualitatively  analyze  the  transient  performance  of 
the  fuel  cell.  The  fuel  cell  is  typically  operated  at  0.6  V  with  a  current  density  of  0.45 
A/cm'  with  the  humidity  at  equilibrium.  If  the  cell  voltage  is  increased  to  0.8  V,  the  ac¬ 
companying  reduction  in  current  to  0.1  A/cm  would  produce  less  water  at  the  cathode, 
which  would  dry  the  cathode  and  the  membrane.  Oppositely,  decreasing  the  voltage  to 
0  4V  would  cause  the  cathode  to  flood,  but  dry  out  the  anode.  The  ohmic  resistance  at 
the  anode  would  increase,  as  would  the  mass  transfer  resistance  at  the  cathode.  Charac¬ 
terising  these  transients  will  be  the  focus  of  Phase  II  work. 


Fuel  Cell  Models 

A  series  of  fundamental  steady-state  and  transient  models  for  fuel  cells  were  de¬ 
veloped.  Some  of  the  models  are  used  to  extract  physical  parameters  from  experimenta 
data,  and  others  are  used  to  predict  the  performance  of  a  fuel  cell. 

Performance  Model 

A  typical  cell  model  is  divided  into  seven  regions:  flow  channels,  gas-diffusion 
layers,  and  catalyst  layers  on  both  the  anode  and  cathode  sides,  and  the  separator.  In  the 
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flow  channels,  there  are  no  reactions  and  the  dominate  transport  mechanism  is  convec¬ 
tion.  In  the  gas-diffusion  layer,  the  dominate  transport  mechanism  is  diffusion  and  there¬ 
fore  the  Stefan-Maxwell  equations  apply.  In  the  catalyst  layer,  both  diffusion  and  reaction 
are  involved.  At  steady  state,  the  time  derivatives  are  zero  and  diffusion  is  balanced  by 
the  reaction  rate.  To  date,  these  equations  have  been  incorporated  into  the  numerical  al¬ 
gorithm,  but  the  performance  curves  have  not  yet  been  obtained. 

Generating  Model  Data 

Model  data  is  generated  from  measurements  on  a  half-cell  as  shown  in  Figure  2. 
The  left  chamber  is  initially  filled  with  nitrogen  gas  and  then  switched  to  oxygen  gas. 
The  oxygen  diffuses  through  the  Nafion®  117  membrane  and  is  reduced  at  platinum  sites. 
The  right  chamber  is  filled  with  an  electrolyte  solution  of  H2SO4  and  acts  as  a  transport 
medium  for  ions  moving  from  the  Pt-working  electrode  to  the  Pt-counter  electrode.  An 
Ag-AgCl  reference  electrode  is  used  as  reference  in  all  experimentation. 

Pt-working 


Figure  2.  The  half-cell  experimental  apparatus  includes  two  chambers 
separated  by  a  Nafion®  117  membrane  platinized  on  one  side, 
shows  the  measured  current  as  a  function  of  time. 

Figure  3  shows  the  measured  current  as  a  function  of  time.  A  potential  of +0.1  V 
was  applied  between  the  platinized  membrane  and  the  Ag/AgCl  electrode,  the  current 
was  recorded  until  a  steady-state  was  reached  (-0.070  pAmps  for  Nafion  1 17  at  25  °C). 
The  experiment  was  then  begun  by  replaced  the  nitrogen  with  oxygen.  The  current  was 
recorded  until  a  new  steady  state  (the  limiting  current,  i*)  was  reached.  For  Nafion  1 17, 
the  experiment  required  six  to  ten  minutes  to  reach  i*.  and  data  were  taken  every  0.50 
seconds  for  all  trials.  After  the  background  current  was  subtracted  from  the  total  meas¬ 
ured  current,  ioo  was  found  to  vary  from  80  to  275  pA  for  temperatures  of  10  °C  to  80  °C, 
respectively. 
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Figure  4  Oxygen  reduction  current  vs.  time  for  Nafion®  117  membrane^at  40  C. 
compares  data  for  an  individual  test  (40  °C)  of  oxygen  diffusion  across  Nafion"  117  and 
the  numerical  fit  using  method  of  least  squares.  In  all  cases,  the  equation  fit  the  data 
well  with  the  margin  of  error  between  the  data  set  and  the  numerical  fit  being  less  than 
2%.  For  higher  temperatures,  the  current  typically  continues  to  increase  slightly  with 

time. 


Figure  4.  Comparison  of  the  experimental  data  for  oxygen  diffusion  in  a 
Nafion®  117  membrane  at  40  °C  to  least  squares  fit. 


System -LEVEL  Model 

The  VTB  system  level  model  of  the  PEM  fuel  cell  was  enhanced  to  include  ef¬ 
fects  such  as;  1)  gas  inlet  pressure  effect  on  voltage-current  characteristics;  2)  Tlierma 
effect  on  voltage-current  characteristics;  and  3)  natural  thermal  connection  for  heat  re- 
moval  While  the  model  contains  the  necessary  complexity  for  system  evaluation,  it  was 
designed  with  ease  of  use  and  high-speed  simulation  in  mind.  Examples  of  the  model 
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output  are  shown  in  Figure  1,  where  the  model  is  operated  at  several  values  of  gas  inlet 
pressure  ranging  from  150  Pa  to  300  Pa  and  at  several  temperatures  ranging  from  -10  to 
+60  C.  The  model  was  validated  against  manufacturer’s  data,  as  shown  in  Figure  2.  Fur¬ 
ther  details  of  the  characteristic  equations  are  available  in  the  fuel  cell  model  help  file. 
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Figure  1.  Pressure  and  temperature  dependence  of  fuel  cell  voltage  and 
current  characteristics. 


Figure  2.  Comparison  between  V-I  characteristics  of  model  and  manufac¬ 
turer’s  data  shows  that  the  model  is  valid. 

HYDROGEN  STORAGE  TECHNOLOGIES 

Characterization  of  hydride  storage 

Development  of  a  mathematical  model  for  the  discharge  of  hydrogen  from  a  metal 
hydride  bed  requires  definition  of  the  thermodynamic  relationship  that  describes  the  equi¬ 
librium  loading  of  hydrogen  in  the  metal  hydride  as  a  function  of  both  pressure  and  tem¬ 
perature.  Two  adsorption  isotherm  relationships  were  developed  in  this  project.  The  first 
one  is  based  on  a  composite  Langmuir  model  and  the  second  one  is  based  on  a  modified 
Virial  model. 

The  prediction  of  the  heat  of  adsorption/reaction  from  these  models  is  critical  to 
the  accurate  modeling  of  the  discharge  behavior  of  hydrogen  from  a  metal  hydride  bed 
since  it  controls  the  temperatures  in  the  bed.  The  models  predict  different  values  for  the 
heat  of  adsoqption/reaction  since  the  composite  Langmuir  model  has  two  values  depend¬ 
ing  on  whether  the  loading  is  above  or  below  the  switch  point,  and  the  modified  Virial 
model  has  a  loading  and  temperature  dependent  function.  Experimental  calorimetric  data 
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are  needed  to  verify  the  heat  of  adsorption/reaction  predicted  from  these  two  different 
thermodynamic  models.  This  issue  will  be  addressed  in  HAPS  II  by  constructing  a  dosing 
calorimeter  that  can  be  used  to  simultaneously  measure  the  loading  and  heat  of  adsorp¬ 
tion/reaction  of  hydrogen  on  metal  hydrides  of  interest  at  different  temperatures.  Pre¬ 
liminary  results  of  this  kind  of  analysis  are  given  below  in  the  System  Level  Model  sec¬ 
tion. 

A  Metal  Hydride  Hydrogen  Storage  Test  Facility,  shown  in  Figure  7,  was  de¬ 
signed  and  constructed  for  the  purpose  of  parameterizing  and  validating  the  model.  It 
consists  of  a  Savannah  River  Technology  Center  metal  hydride  bed  containing  26.0  kg  of 
Lmi.o6Ni4.96Alo.o4  and  flow  control  elements.  Six  runs  were  carried  out  at  hydrogen  flows 
varying  from  5  to  40  SLPM. 

Atmosplieric  venf 


Figure  7.  Schematic  of  the  Metal  Hydride  Hydrogen  Storage  Test  Facility. 

The  complexity  of  this  experimental  system  represents  a  significant  challenge  to  any 
modeling  effort.  Many  more  experiments  need  to  be  carried  out  under  different  initial 
and  process  conditions  to  properly  validate  the  various  assumptions  used  in  the  mathe¬ 
matical  model,  and  to  verify  the  thermodynamic  and  transport  properties  .  This  task  will 
be  carried  out  in  HAPS  II. 

Detailed  Model 

The  model  that  describes  the  dynamic  discharge  of  hydrogen  from  a  metal  hy¬ 
dride  bed  has  two  components:  1)  general  equations  that  describe  the  dynamics  mside 
the  metal  hydride  bed,  together  with  different  assumptions  that  the  user  can  utilize  to 
simplify  the  general  equations  according  to  the  mathematical  capabilities  of  the  simulator 
and  2)  equations  that  describe  the  dynamics  involved  in  a  gas  pressure  regulator  placed 
downstream  from  the  metal  hydride  bed  and  a  ready  to  use  simplified  model  for  initial 
use  in  the  virtual  test  bed  (VTB)  has  been  created  and  will  be  implemented  in  HAPS  II. 
Based  on  the  above  developments,  a  layout  of  the  system  described  mathematically  is 
presented  in  Figure  8.  It  consists  of  a  cylindrical  metal  hydride  bed  connected  in  series 
with  a  gas  pressure  regulator  downstream  to  provide  hydrogen  at  a  relatively  constant 
pressure  Pfc  to,  for  example,  a  fuel  cell. 
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Figure  8.  Schematic  of  the  system  described  mathematically  including  the 
metal  hydride  bed  and  an  inline  gas  regulator,  and  typical  profiles  of  the 
main  variables  involved  in  the  discharge  of  hydrogen  from  the  bed. 

The  complexity  of  the  complete  model  is  user  specified  and  depends  on  the  capa¬ 
bility  of  the  simulator.  The  most  rigorous  model  accounts  for  various  mass  and  heat 
transfer  resistances  in  both  the  axial  and  radial  directions  and  accounts  for  internal  heat 
exchange  with  a  circulating  fluid.  The  simplest  model  ignores  resistances  and  isothermal 
local  equilibrium  conditions  prevail.  These  bounding  models  represent  extreme  limits 
and  can  be  used  for  actual  process  design  and  hypothetical  process  feasibility,  respec¬ 
tively.  They  can  also  be  used  to  judge  complexity  of  a  particular  H2  discharge  process  by 
comparing  with  experimental  data. 

Results  from  various  models  utilizing  various  assumptions  have  not  yet  been 
compared  to  experimental  data  in  a  purely  predictive  mode.  This  requires  obtaining  the 
proper  thermodynamic  and  transport  relationships  that  are  needed  as  input  to  the  model 
and  obtaining  sufficient  dynamic  discharge  data  that  can  be  used  to  validate  the  various 
models  utilizing  various  simplifying  assumptions.  This  will  be  another  focus  of  the  work 
to  be  carried  out  in  HAPS  II.  It  is  anticipated  that  even  this  more  rigorous  model  or  one 
very  similar  to  it  can  be  used  in  the  VTB  simulator)  hence,  it  will  constitute  a  major  com¬ 
ponent  in  the  design  and  development  of  a  hybrid  advanced  power  source  via  the  VTB. 


System-level  Model 

A  simplified  dynamic  model  of  the  discharge  of  hydrogen  from  a  cylindrical 
metal  hydride  bed  has  been  developed  for  use  in  the  Virtual  Test  Bed  process  simulator. 
The  following  simplifying  assumptions  were  utilized  in  this  model:  no  heat  transfer  resis¬ 
tance  at  the  interface,  plug  flow,  no  internal  cooling  system,  negligible  heat  loss  by  con¬ 
vection,  negligible  thermal  conductivity,  perfectly  insulated  system,  and  linear  flux  drop 
in  the  bed.  This  model  still  accounts  for  heat  and  mass  transfer  resistances  through  over¬ 
all  (or  lumped)  heat  and  mass  transfer  coefficients.  At  present,  however,  these  two  coef- 


10 


'  On-Board  Fuel  Cell  Power  Plants 


Final  Report 


ficients  must  be  obtained  by  fitting  the  model  to  experimental  discharge  data  as  discussed 
below.  This  dynamic  model  is  currently  being  implemented  in  the  VTB. 

Preliminary  experimental  data  obtained  from  the  experimental  metal  hydride  hy¬ 
drogen  storage  test  facility  described  above  were  used  to  calibrate  this  simple  model  by 
obtaining  a  set  of  heat  and  mass  transfer  coefficients  from  one  run  and  applying  them  to 
predict  the  behavior  of  other  runs.  Note  that  the  modified  Virial  model  was  used  in  these 
simulations.  Only  three-fourths  of  the  cylinder  is  filled  with  the  metal  hydride  to  allow 
the  material  to  expand  when  taking  up  hydrogen  during  charging.  As  a  result,  it  was  an¬ 
ticipated  that  this  simplified  first  generation  VTB  model  will  not  be  able  to  accurately 
predict  the  dynamic  behavior  of  this  system,  especially  the  heat  effects.  Nevertheless, 
preliminary  modeling  results  of  the  dynamic  discharge  of  hydrogen  from  this  system  un¬ 
der  different  molar  discharge  flow  rates  were  quite  encouraging  (Figure  9).The  model 
does  very  well  at  the  low  flow  rate  demand,  but  fails  at  higher  flow  rate  demand  (not 

shown). 
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Figure  9.  Comparison  of  the  experimental  pressure  history  with  different 
model  predictions  during  discharge  at  5  SLPM  H2  demand. 

Figure  1 0  shows  data  obtained  from  VTB  simulation  of  the  simplified  model  that 
discharges  hydrogen  to  a  constant  sink  (fuel  cell).  The  graphs  show  variation  of  the  pres¬ 
sure,  temperature  and  molar  flow  rate  of  the  hydride  bed,  H2  concentration  inside  the  hy¬ 
dride  bed.  From  Fig  10,  it  is  shown  that  the  pressure  declines  rapidly  initially  from  15.0 
atms  and  then  decreases  slowly  after  the  pressure  drops  below  5.0  atms.  It  can  also  be 
seen  that  the  temperature  and  H2  concentration  decrease  linearly  with  the  time  when  the 
molar  flow  is  constant. 

Currently,  the  effort  to  refine  the  model  is  continued  in  conjunction  with  other 
projects.  It  will  be  further  validated  by  comparing  to  the  experimental  data. 
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Figure  10.  Simulation  results  of  the  hydride  bed  model  under  a  constant 
hydrogen  flow  sink  load. 

CENTRIFIGUAL  COMPRESSOR 

We  have  developed  several  compressor  models  concerning  different  applications 
and  different  research  foci.  The  simple  and  linear  model  is  based  on  empirical  data  and 
valid  near  and  within  a  limited  range  of  the  operating  point.  The  output  pressure  is  a  lin¬ 
ear  function  of  mass  flow  and  rotational  speed.  The  compressor  efficiency  is  a  constant 
that  users  can  input  based  on  their  field  experience  in  steady  state  operation. 

The  latest  development  is  a  nonlinear  centrifigual  compressor  model  derived  from 
the  first  principle,  based  on  which  energy  transfer  considerations  are  used  to  develop 
compressor  characteristics.  The  centrifugal  compressor  model  is  presented  jointly  by  the 
states  of  mass  flow,  pressure  rise  and  rotational  speed  of  the  spool.  Pressure  nse  is  de¬ 
rived  from  the  temperature  changes  in  the  compression  process.  Accordingly,  ij-essure 
rise  is  further  related  to  energy  transfer  by  considering  the  temperature  chanp.  The  ac¬ 
tual  energy  transfer  in  the  compressor  is  equivalent  to  ideal  energy  transfer  minus  energy 
losses.  There  are  two  major  losses:  incidence  losses  and  friction  losses.  Incidence  Losses 
are  the  losses  due  to  incidence  onto  the  rotor  and  vane  diffuser.  Frictional  Losses  are  the 
losses  due  to  friction  in  the  impeller  and  diffuser.  The  variable  speed  compressor  charac¬ 
teristic  is  developed  based  on  these  losses  consideration. 
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The  new  model  valid  in  a  wide  range.  Rotational  speed  should  be  over  300  rpm  to 
guarantee  convergence.  The  isentropic  efficiency  of  the  compressor  varies  at  different 
operating  points.  It  is  defined  as: 

Efficiency  =  (ideal  energy  transfer-energy  losses)  /  ideal  energy  transfer 

Furthermore,  the  efficiency  will  be  corrected  with  losses  in  the  volute  and  additional 
losses  arising  from  clearance  and  backflow. 

Figure  1 1  shows  the  simulation  results  of  the  compressor  transient  and  steady 
state  operation,  under  the  condition  of  a  constant  load.  The  model  is  under  validation 
now. 
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Figure  11.  Transient  and  steady  state  of  the  VTB  compressor  model. 

METHANOL  REFORMER 

The  reformer  is  assumed  to  be  a  CSTR  model  i.e.  the  temperature,  concentrations 
variations  inside  the  bed  are  neglected.  The  heat  input  to  the  reformer  is  assumed  to  be 
provided  by  a  proportional  controller  that  is  within  the  reformer  and  is  not  shown  explic¬ 
itly.  The  reformer  is  designed  to  operate  for  steam  to  methanol  ratio  of  0.5  to  2.0,  150  C 
to  266°  C  temperature  and  2  to  3  bar  pressure. 

The  model  allows  one  to  specify  the  nominal  power  rating  of  the.  fuel  cell  stack 
that  the  reformer  will  feed  (in  order  to  estimate  the  hydrogen  production  rate  and  hence 
internal  component  sizes),  the  target  operating  temperature  and  pressure,  and  the  design 
conversion  ratio.  The  reformer  model  does  not  yet  incorporate  explicit  closed-loop  con- 


V  -  V 


'■  . ^ 

1  LJcomO_Mass_flow j-ate 

r  -  '  ;  JnJidI 

Picked;  none  Independent:  tinne  Triggering:  disabled  Frame  rate:  1 .0  Hz 

~ 

_ ^ _ 

I 

i 

f0.1 

- d:oi 

0.03;  ■'  o:o4|  o.os 

HOB 

0,0/j  u.Ubj  uua 

_ U 

► 

-  ^ 

Picked:  non 

e  Independent:  time  Triggering:  disabled  Frame  rate:  0.7  Hz 

140000 

— 

= 

130000 

! 

120W^ 

! 

li«1^000.01 

0.02 

0.03 

0.04 

0.0^ 

0.06 

0.07 

0.08'  0,09 

-  - ^ _ u 

► 

-  : _ 

P  -loixl 

Picked:  none  Independent:  time  Triggering:  disabled  Frame  rate:  5.0  Hz 

0.8/" - 

_ 

O.f 

oi 

(j|2  0.01 

0.02 

0.03 

0.04 

0.05 

0,06 

0.07 

0.08 

0.09 

_ _ _ _ ^ 

► 

_ _ _ _ Ul 

13 


’  On-Board  Fuel  Cell  Power  Plants 


Final  Report 


trol  functions  to  control  its  operation,  but  instead  it  is  forced  to  operate  at  near  nominal 
operating  conditions.  Hence  the  hydrogen  gas  production  often  varies  from  the  consump¬ 
tion  rate  of  the  fuel  cell,  depending  on  system  load.  Excess  production  is  stored  in  a  hy¬ 
drogen  tank,  and  insufficient  production  is  made  up  for  by  drawing  from  the  tank. 

Figure  12  shows  the  simulation  results  of  the  reformer  feeded  by  a  constant 
methanol  source  and  loaded  with  a  hydrogen  storage  tank.  The  fractional  conversion  (up- 
the  per  left)  of  methanol  in  the  reformer  bed  changes  as  a  function  of  temperature.  Final 
value  is  0.95.  The  temperature  rise  (the  upper  right)  of  the  reformer  bed  from  room  tem¬ 
perature,  30  °C  to  the  operating  temperature,  266  °C  with  time.  H2  (the  lower  left)  ob¬ 
tained  from  the  reformer  exit  in  moles/s.  Steady  state  value  is  0.187  moles/s.  The  heat 
input  in  watts  is  shown  by  the  lower  right  graph.  The  steady  state  value  is  13700  watts 

(corresponding  to  a  value  of  81.69  W/°C). 
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Figure  1 2.  VTB  simulation  results  of  the  reformer  that  connects  to  a  con¬ 
stant  methanol  source  and  discharges  to  a  hydrogen  storage  tank. 


FEEDBACK  CONTROL  SYSTEM 

The  feedback  control  system  is  designed  using  MatLab/Simulink,  as  shown  in 
Figure  13.  It  represents  a  simple  PID  controller  with  a  steady-state  offset.  This  controller 
measures  the  current  that  is  produced  by  the  fuel  cell,  and  calculates  the  necessary  air- 
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flow  for  that  operating  condition.  The  “Current  Scale”  is  the  product  of  357*10  kg/amp- 
cell,  the  number  of  cells  (1600),  and  the  current  that  is  measured.  Normally,  a  fuel  cell  is 
feed  anywhere  between  50%  and  200%  excess  airflow  for  cooling  and  humidity  control. 
For  this  simulation,  100%  airflow  is  chosen.  Thus,  the  difference  between  the  excess  air¬ 
flow  and  the  necessary  airflow  is  the  error  of  the  set  point.  This  value  is  then  passed  to  a 
proportional-integral-derivative  controller,  which  is  tuned  intuitively.  There  is  also  a 
steady  state  offset  (“Rest  Duty  Cycle”)  summed  to  this  value  to  ensure  the  compressor 
motor  is  idling,  and  needs  not  start  from  a  dead  stop.  The  “Saturation”  block  is  simply  to 
limit  the  duty  cycle  to  a  reasonable  range  (0-1).  The  input  and  output  block  numbers  cor¬ 
respond  to  the  input  and  output  pins  of  the  “Simulink  Controller”  object  in  VTB. 


Figure  13.  The  PID  controller  design  using  Simulink. 


SYSTEM  MODEL 

Figure  14  illustrates  the  system  configuration  as  represented  in  the  Virtual  Test 
Bed  environment.  Salient  components  of  the  system  model  include  the  fuel  cell  stack,  the 
fuel  reformer,  a  motor-driven  compressor,  motor  drive  electronics,  and  a  controller  that 
senses  load  current  and  stack  exhaust  gas  flow.  The  other  auxiliaries  include  the  thermal 
environmental  elements  and  the  air  bypass  valve.  The  system  currently  is  not  rated  yet 
since  relative  little  information  on  the  practical  power  plant  is  available. 
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Figure  14.  Schematic  diagram  of  the  fuel  cell  power  plant  in  the  VTB  en¬ 
vironment. 


The  fuel  cell  stack  can  be  specified  for  its  number  of  series  and  parallel  cells,  and 
the  cell  area.  The  heat  generated  inside  the  stack  is  stored  in  an  external  thermal  capaci¬ 
tance  that  represents  the  lumped  heat  capacity  of  the  fuel  cell  stack.  Thermal  energy  is 
dissipated  to  the  ambient  through  an  external  thermal  conductance.  The  aerodynamics  of 
air  flow  fields  are  not  represented. 

The  compressor  model  accounts  for  the  inlet  and  outlet  gas  temperatures,  heating 
of  the  gas,  efficiency,  and  speed-dependent  pressure  rise.  It  can  be  used  in  either  direc¬ 
tion,  either  as  a  pump  or  as  an  expander.  The  latter  allows  exhaust  gases  from  the  fuel 
cell  to  be  recuperated  and,  with  reheat,  to  run  the  air  feed  pump.  That  energy  saving  fea¬ 
ture  is  not  yet  incorporated  in  the  system  model,  but  the  component  models  will  now 
support  that  activity. 

The  reformer  model  represents  steam  reformation  of  methanol  into  hydrogen.  The 
model  allows  one  to  specify  the  nominal  power  rating  of  the  fuel  cell  stack  that  the  re¬ 
former  will  feed  (in  order  to  estimate  the  hydrogen  production  rate  and  hence  internal 
component  sizes),  the  target  operating  temperature  and  pressure,  and  the  design  conver¬ 
sion  ratio.  The  reformer  model  does  not  yet  incorporate  explicit  closed-loop  control  func¬ 
tions  to  control  its  operation,  but  instead  it  is  forced  to  operate  at  near  nominal  operating 
conditions.  Hence  the  hydrogen  gas  production  often  varies  from  the  consumption  rate  of 
the  fuel  cell,  depending  on  system  load.  Excess  production  is  stored  in  a  hydrogen  tank, 
and  insufficient  production  is  made  up  for  by  drawing  from  the  tank. 
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The  Simulink  control  object  that  controls  the  power  converter  that  operates  the 
compressor  drive  motor.  The  VTB  environment  now  includes  sophisticated  links  to  for¬ 
eign  simulation  objects  (such  as  Simulink),  allowing  the  use  of  either  compiled  or  inter¬ 
active  co-simulation  models  in  the  system  definition. 

Figures  15  and  16  show  results  from  a  typical  system  study.  Here,  we  show  the 
response  of  the  system  to  a  recurring  3-step  change  of  the  load  from  20  kW  for  20  sec¬ 
onds,  increasing  to  600  kW  for  40  seconds,  then  reducing  to  300  kW  for  the  next  20  sec¬ 
onds,  and  repeating. 


Figure  15.  From  top  left  to  bottom  right,  exit  pressure  of  the  air  pump  in 
Pa,  air  consumption  rate  in  kg/sec,  compressor  speed  in  rad/s,  and  elec¬ 
tric  power  consumed  by  the  motor  in  W. 


One  can  observe  several  interesting  phenomena  in  these  graphs.  When  the  load 
demands  full  power  (600  kW),  the  stack  voltage  drops  from  approximately  1500  V  open 
circuit  to  850  volts.  This  suggests  a  need  for  extensive  power  conditioning  between  the 
stack  and  the  load  bus.  This  power  conditioning  has  not  yet  been  incorporated  into  the 
system  model.  The  stack  temperature  increases  by  approximately  40  degrees  C.  The  rate 
of  heating  is  not  constant,  but  follows  somewhat  the  air  mass  flow.  The  air  pump  motor 
consumes  a  significant  fraction  of  the  system  power,  approximately  300  kW  when  the 
load  draws  600  kW,  but  even  more  when  the  load  draws  only  300  kW  (though  in  the  ex¬ 
ample,  it  never  reaches  steady-state  at  the  300  kW  power  level.)  The  lag  time  of  the  mo¬ 
tor  speed  in  responding  to  the  fuel  cell  stack  air  demands  accounts  for  a  significant  part  of 
the  system  dynamics.  Heating  of  the  fuel  cell  stack  has  a  noticeable  affect  on  system  per- 
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formance.  The  cooling  effect  of  excess  air  blowby  is  evident  following  each  reduction  in 
power  demand,  and  the  compressor  motor  can  be  seen  to  overspeed  in  those  conditions. 


The  system  is  not  yet  optimally  designed,  nor  that  the  control  algorithms  are  fully 
developed,  so  all  of  these  results  should  be  taken  simply  as  examples,  rather  than  as  de¬ 
finitive  predictions  of  the  performance  of  a  whole  class  of  power  plants.  We  have  not  yet 
taken  pains  to  size  all  of  the  components  correctly  in  relation  to  each  other  nor  to  size 
them  carefully  with  respect  to  the  total  load  power  demand. 

In  the  future  research,  we  expect  to  more  fully  explore  and  understand  methods 
for  handling  mixed  streams  of  compressible  gases  within  the  VTB  environment,  to  fur¬ 
ther  improve  the  fuel  cell  stack  and  compressor  models,  to  augment  the  system  with 
power  conversion  equipment  capable  of  regulating  the  bus  voltage,  and  to  incorporate 
additional  features  of  the  thermal  systems.  We  will  also  study  designs  that  use  exhaust 
gas  recuperators  and  expanders  to  power  the  inlet  air  pump.  In  addition,  we  will  attempt 
to  validate  all  models  against  real  test  data,  where  such  data  are  available  for  comparable 
systems.  Note  that  there  are  no  PEM  fuel  cell  plants  operating  at  the  MW  level,  so  there 
is  a  dearth  of  test  data  to  use  for  model  validation.  However,  comparisons  can  be  made  to 
systems  operating  at  lower  power  levels. 
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DELIVERABLES 

Models 

VTB  Models 

The  following  models  were  developed.  Full  descriptions  of  each  appear  in  Ap¬ 
pendix  I  in  the  form  of  the  model  help  files. 


Model  ID 

Description 

Features 

Fuel  Cell  12.15_010716 

Fuel  cell 

Temperature  dependence, 
heat  power  generation, 
pressure  dependence, 

Nemst  equation  based  volt¬ 
age  and  current  characteris¬ 
tics,  fuel  and  air  consump¬ 
tion 

Compressor! 

General  centrifigual  com¬ 
pressor 

The  model  presents  states  of 
mass  flow,  pressure  rise  and 
rotational  speed  of  the  spool 
based  on  the  energy  transfer 
considerations.  Temperature 
change  and  efficiency  of  the 
isentropic  compression  are 
considered.  Incidence  loss 
and  frictional  loss  are  in¬ 
cluded. 

Creformer 

Steam  methanol  reformer 

This  is  a  CSTR  model  of  the 
reformer.  The  temperature  of 
the  catalyst  bed,  methanol 
conversion  and  the  pressure 
inside  the  reactor  change  with 
time  but  do  not  vary  within  the 
reformer.  External  heat  input 
to  the  reformer  is  also  in¬ 
cluded  to  control  the  reformer 
bed  temperature  at  the  desired 
value. 

Hydride  Bed 

Cylindrical  metal  hydride 
bed 

Mass  and  heat  transfer  re¬ 
sistances  ignored.  Isother¬ 
mal  local  equilibrium  as¬ 
sumed.  Heat  exchange  of 
bed  with  the  external  heat¬ 
ing/cooling  system  indued. 

It  can  be  connected  in  series 
with  a  gas  pressure  regula¬ 
tor  downstream  to  provide 
hydrogen  at  a  relatively 
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constant  pressure. _ 

5  in_5out_Matlab 

PID  controller  (Mat- 
Lab/Simulink  object) 

Simple  PID  controller  de¬ 
signed  in  MatLab/Simulink. 
The  Simulink  PID  object 
embedded  in  VTB  allows 
co-simulation  of  the  Simu¬ 
link  with  VTB  time  domain 
solver 
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PUBLICATIONS  AND  PRESENTATIONS 

Publications 

Journals 

In  print 

1.  A.  T.  Haug,  J.W.  Weidner,  R.E.  White,  and  W.  Huang  “Development  of  a  Novel 
CO  Tolerant  Anode,”  accepted,  J.  Electrochem.  Soc.,  December,  2001. 

2.  A.  T.  Haug,  J.W.  Weidner,  R.E.  White,  W.  Huang,  Steven  Shi,  Timothy  Stoner, 
and  Narender  Rana,  “Using  Sputter  Deposition  to  Increase  CO  Tolerance  in  a 
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4.  A.  T.  Haug,  R.E.  White,  J.W.  Weidner,  W.  Huang,  Steven  Shi,  Timothy  Stoner, 
and  Narender  Rana,  “Increasing  PEM  Fuel  Cell  Catalyst  Effectiveness 
Through  Sputter  Deposition,”/.  Electrochem.  Soc.,  149,  A280-A287  (2002). 


Presentations 

Presentations  are  listed  here  only  if  there  was  not  a  corresponding  publication  in 
conference  proceedings.  Otherwise,  see  above  list  of  Proceedings  . 

1.  Huag.  A..  J.  W.  Weidner,  R.  E.  White,  and  W.  Huang  “Development  of  a  Novel 
CO  Tolerant  Anode,”  The  Electrochemical  Society,  San  Francisco,  CA,  Sep¬ 
tember,  2001. 

2.  Huag.  A..  R.E.  White,  J.W.  Weidner,  W.  Huang,  Steven  Shi,  Timothy  Stoner,  and 
Narender  Rana,  “Increasing  PEM  Fuel  Cell  Catalyst  Effectiveness  Through 
Sputter  Deposition,”  The  Electrochemical  Society,  Philadelphia,  PA,  May,  2002. 

3.  Blackwelder.  M.  “Fuel  Cell  Model  Implementation  and  System  Evaluation”, 
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CONCLUSIONS 

Major  component  models  of  the  shipboard  fuel  cell  power  plant  were  developed. 
A  example  system  for  the  power  system  was  assembled  and  it  can  be  used  for  study  of 
system  first  order  transient  and  steady  state  operation.  The  future  work  will  focus  on  the 
validation  of  all  component  models  against  experimental  data,  and  protype  the  system  in 
VTB  environment  with  a  feasible  parameters,  so  that  the  system  design  optimization  and 
operation  performance  can  be  obtained. 
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APPENDICES 


Appendix  1  -  Students  Education 

Foilwing  table  lists  the  graduate  students  who  were  or  have  been  supported  by  the 
fuel  cell  shipboard  power  plant  project. 


Name 

Pursuing  Degree 

Major  Field 

Advisor 

Enrolling  Status 

Mark  Blackwelder 

Ph.  D. 

Electrical 

R.  A.  Dougal 

Current 

Ramasrishna  Gun- 
dala 

Ph.  D. 

Chemical 

J.  Weidner 

Current 

Andrew  Haug 

Ph.  D. 

Chemical 

R.  E.  White 

Current 

Ananda  Mondal 

Ph.  D. 

Graduated 

Venka  Subraminian 

Ph.  D. 

Chemical 

R.  E.  White 

Graduated 

Wei  Jiang 

Ph.  D. 

Mechanical 

J.  Khan  /  R.  A. 
Dougal 

Current 

Zhenhua  Jiang 

Ph.  D. 

Electrical 

R.  A.  Dougal 

Current 

Ramaraja  Rama- 
samy 

Ph.  D. 

Chemical 

J.  Weidner 

Current 

Vijay  Sethuraman 

Ph.  D. 

Chemical 

J.  Weidner 

Current 

Appendix  2  -  Fuel  Cell  CO  Modeling 

See  attached  paper  entitled  “A  Mathematical  Model  Predicting  the  Effects  of  Re¬ 
formate  and  Air  Bleed  on  a  PEM  Fuel  Cell”  by  Andrew  T.  Haug,  John  W.  Weidner,  and 
Ralph  E.  White. 

Appendix  3  -  Model  Help  Files 

Description  of  model  help  files 

Appendix  4  -  Publications  (copies  attached) 

1.  A.  T.  Haug,  J.W.  Weidner,  R.E.  White,  and  W.  Huang  “Development  of  a  Novel 
CO  Tolerant  Anode,”  J.  Electrochem.  Soc.,  149  (7)  A862-A867  (2002). 

2.  A.  T.  Haug,  J.W.  Weidner,  R.E.  White,  W.  Huang,  Steven  Shi,  Timothy  Stoner, 
and  Narender  Rana,  “Using  Sputter  Deposition  to  Increase  CO  Tolerance  in  a 
PEM  Fuel  Cell,”  J.  Electrochem.  Soc.,  149  (7)  A868-A872  (2002). 

3.  D.  Tianping,  J.W.  Weidner,  and  R.E.  White,  “Extension  of  Newman’s  Method 
to  Electrochemical  Reaction-Diffusion  in  a  Fuel-Cell  Catalyst  Layer,”/. 

Power  Sources,  107  (1),  24-33,  2002. 

4  A  T  Haug,  R.E.  White,  J.W.  Weidner,  W.  Huang,  Steven  Shi,  Timothy  Stoner, 
and  Narender  Rana,  “Increasing  PEM  Fuel  Cell  Catalyst  Effectiveness 
Through  Sputter  Deposition,”  /  Electrochem.  Soc.,  149,  A280-A287  (2002). 
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Abstract 

A  model  has  been  developed  to  explain  the  effects  of  introducing  an  air  bleed  to  a 
PEMFC  anode  feed  containing  hydrogen  and  carbon  monoxide.  The  basis  for  this  model 
was  the  work  by  Springer  et  al  (32)  for  a  PEMFC  anode  stream  containmg  hydrogen  and 
CO  The  appropriate  coverage  and  flux  equations  have  been  modified  to  take  into 
account  the  electrochemical  oxidation  of  CO  that  has  been  proven  m  earlier  work  to  be 
the  primary  means  by  which  oxygen  oxidizes  CO.  This  model  shows  that  the  CO 
electrooxidation  rate  is  greatly  improved  through  the  addition  of  oxygen  to  an  anode 
feedstream.  The  model  shows  a  good  quantitative  fit  to  performance  data  taken  under 
various  air-bleed  concentrations  for  a  PEMFC  containing  a  Pt  anode.  The  rnodel  then 
further  explores  the  selectivity  of  the  electrochemical  reaction  of  oxygen  with  CO  versus 
the  catalytic  oxidation  with  H2.  Finally,  qualitative  projections  are  made  on  the  required 
increases  in  certain  parameters  that  would  result  in  improving  the  selectivity  of  O2  usage 
in  the  anode  feedstream. 


Introduction 

Proton-exchange  membrane  fuel  cells  (PEMFCs)  are  gaining  popularity  due  to 
their  benefits  such  as  environmental  friendliness  and  increased  fuel  efficiency.  Because 
of  the  difficulties  inherent  to  storing  hydrogen,  liquid  fuels  such  as  propane,  natural  gas 
and  gasoline  are  used  to  produce  reformate  gas.  Dry  reformate  is  typically  composed  of 
35  -  45  %  hydrogen,  15-25  %  carbon  dioxide,  50  -  10000  ppm  carbon  monoxide  and  a 
balance  of  nitrogen.  It  has  been  shown  extensively  that  CO  poisons  the  platinum  catalyst 
used  in  PEMFC  systems  (1,2,3).  Fuel  clean-up  units  including  high  and  low  temperature 
water-gas  shift  (WGS)  reactors  and  a  preferential  oxidation  (PROX)  unit  typically  reduce 
CO  levels  to  10-100  ppm  in  the  fuel  stream.  However,  at  concentrations  as  low  as  10 
ppm,  CO  lowers  power  output  of  the  PEMFC  containing  Pt  electrodes  by  50%  (4,5). 

Springer  et  al  (32)  have  developed  a  model  that  rigorously  investigates  the  effects 
of  CO  on  the  PEMFC  anode  for  reformate  feedstreams.  This  includes  the  effects  of  H2 
dilution  and  starvation  and  the  effects  of  CO  on  the  Pt  anode  performance. 
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Alloying  Pt  with  various  metals  (8,9,10,11,12),  increasing  the  fuel  cell  operating 
temperature,  and  adding  oxygen  to  the  anode  feedstream  (6,7)  are  all  being  investigated 
as  means  of  increasing  CO  tolerance  (13,14,15,16,17,18).  The  injection  of  oxygen  into 
the  fuel  stream  has  also  been  shown  to  increase  catalyst  tolerance  to  CO  (6,7,19,20). 
Injected  as  an  ‘air  bleed’,  this  method  provides  a  greater  concentration  of  oxygen  on  the 
catalyst  surface  that  then  reacts  with  CO  to  form  CO2.  It  has  been  proven  that  the 
primary  means  by  which  oxygen  oxidizes  CO  on  the  catalyst  surface  is  by  reacting  with 
adsorbed  hydrogen  to  form  a  hydroxyl  species  that  then  reacts  electrochemically  with 
adsorbed  CO  to  form  CO2  (19,20).  It  is  these  kinetics  occurring  within  a  reformate  feed 
containing  oxygen  that  are  modeled  here. 

To  develop  a  model  predicting  the  steady-state  performance  of  a  PEMFC  in  the 
cases  where  CO  +  H2  and  CO+H2+O2  are  present  in  the  anode  feed  stream,  two  basic 
components  are  required:  the  necessary  kinetic  expressions  to  accurately  describe  the 
effects  of  CO  on  the  PEMFC  anode  and  the  necessary  governing  equations  for  the 
various  regions  of  the  PEMFC  (anode  GDL,  anode  catalyst  region,  membrane,  cathode 
catalyst  region  and  cathode  GDL). 

Those  necessary  equations  are  presented  here.  Three  separate  derivations  of 
expressions  for  CO  interactions  with  Pt  are  given  based  on  work  by  Vogel,  Dhar  and 
Springer  (21,22,23,32).  The  1-D  version  of  Springer’s  model  is  then  described  and  used 
to  explain  the  Plug  Power  steady  state  data  in  the  case  of  a  reformate  stream  containing 
CO  interacting  with  a  Pt  anode.  Springer’s  work  is  then  modified  to  describe  the  effects 
of  an  airbleed  as  shown  in  previous  chapters. 
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Model  Development 


ASSUMPTIONS 

Below  are  the  assumptions  for  the  1-D  (through  the  cell)  model  predicting  the 
steady-state  performance  of  a  PEMFC  membrane-electrode  assembly.  These  are  valid 
for  both  Springer’s  model  and  the  modified  version  that  accounts  for  the  addition  of 
oxygen  to  the  anode  feedstream. 

•  One  dimensional  flow  of  the  of  the  species  through  the  anode  (24,25). 

•  The  anode  operates  under  isothermal  conditions  with  uniform  temperature 
distribution  (24,26).  This  assumption  is  based  on  the  high  thermal  conductivities  and 
of  the  carbon  portion  of  the  electrode  and  even  mixing  that  occurs  in  flow  through 
porous  media  (24,27). 

•  Species  are  distributed  evenly  throughout  the  cross-sectional  area  of  the  anode. 

•  Ha,  CO,  and  water  vapor  will  be  treated  as  ideal  gases  that  are  well  mixed  (24,26). 

•  The  porosity  of  the  GDL  and  in  the  ACR  is  considered  uniform. 

.  The  catalyst  in  the  ACR  is  assumed  to  be  spread  uniformly  throughout  the  region. 

•  Pressure  is  considered  uniform  throughout  the  cell. 

.  The  inlet  gas  is  considered  fully  hydrated  (24,25,27,28),  however  any  liquid  water 
entering  the  PEMFC  anode  is  neglected. 

•  Reactions  only  occur  in  the  active  catalyst  regions  (ACRs)  of  the  PEMFC.  This 
reaction  is  considered  homogeneous. 

•  It  is  assumed  that  the  entire  Pt  particle  area  is  active. 

•  The  CO2  and  N2  present  in  the  anode  ACR  are  considered  inert. 

•  The  rate  of  reaction  in  the  ACR  is  mass  transfer  limiting  (29,30). 

H2  +CO  IN  THE  ANODE  FEED 

In  order  to  generate  an  expression  relating  anodic  performance  to  the  amount  of 
CO  in  the  feed  stream,  the  derivation  must  first  be  made  of  the  current  generated  by  the 
dissociation  of  H2  in  terms  of  surface  coverage  of  the  Pt  catalyst  (23,3 1): 


J  =  Jo 


F^y/RT 


^H2  ^F/7/RT 


■"H2 


(0.1) 
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where  0h  is  the  steady  state  hydrogen  coverage,  jo  is  the  exchange  current  density 
in  A  cm■^  Ch2  and  c^^  are  the  initial  and  local  hydrogen  concentrations  in  mol  cm■^  r\  is 
the  anode  overpotential.  A  more  detailed  derivation  is  presented  in  Appendix  A. 

Vogel  et  al  (23)  deduced  the  empirical  relationship  shown  in  equation  (0.3)  and 
(0.2)  for  relating  the  reaction  rates  and  current  density  of  pure  Ha  to  that  of  gas  containing 
CO  and  Ha  in  phosphoric  acid  fuel  cells. 

kf°  =  k.(l-«co)’ 


=  (0.4) 

Jh2 


where  and  are  the  rate  constants  for  Ha  adsorption  onto  Pt  for  feedstreams 

containing  pure  Ha  and  Ha+CO,  respectively.  The  performance  loss  by  CO  arises  due  to 
its  ability  to  block  reaction  sites  which  otherwise  would  have  been  accessible  to  Ha.  This 
is  distinct  from  the  Nemst  loss,  which  arises  because  of  reactant  dilution,  and  is  the  basis 
for  the  relationship  with  Oco  (21).  Thus,  a  final  expression  for  the  current  density  can  be 
found,  expressed  in  terms  of  the  steady-state  hydrogen  concentration,  Ch^  ,  equilibrium 

hydrogen  coverage,  0cO)  &nd  CO  coverage,  0co. 


where  jH2  w/co  the  current  generated  when  CO  is  present  in  the  feed  stream. 

Dhar,  Christner,  et  al  (6)  developed  an  expression  for  the  coverage  of  CO  on 
PAFCs  at  over  the  range  of  130  to  190°C: 


6 


r 


^co  ~ 


-AG“  RT,  RT, 
- 5. - ln(H)  + — In 

r  r  ^  '  r 


^co 

) 


(0.6) 


ACp  rn  AH  /n  n\ 

lnH  =  A  +  ^lnT-—  (0-7) 

where  AGo°  is  the  standard  free  energy  of  adsorption  for  CO  (kcal  mol"^),  r  is  an 
interaction  parameter  (kcal  mol'^)  and  H  is  Henry's  constant  (atm  mol'^  L'^)  for  CO 
absorption  in  water,  ACp  and  AH°  are  the  change  in  heat  capacity  and  standard  change  in 
enthalpy,  respectively  of  CO  oxidation  and  A  is  a  constant. 

The  effects  of  CO  on  fuel  cell  performance  increase  dramatically  as  temperature 
decreases  from  PAFC  conditions  to  typical  PEMFC  conditions  (33).  While  Dhar  and 
Vogel’s  work  might  be  adaptable  for  typical  PEMFC  operating  temperatures  (40  •  100 
°C),  the  effects  of  CO  on  fuel  cell  performance  increase  dramatically  as  the  temperature 
is  decreased  from  PAFC  conditions  to  PEMFC  conditions  (33).  Springer  et  al  (32)  have 
developed  solutions  specifically  to  explain  the  effects  of  CO  poisoning  on  the  PEMFC 
anode.  A  detailed  description  of  this  model  is  described  in  Appendix  C.  This  model  is 
both  rigorous  in  the  description  of  the  PEMFC  anode  catalyst  region  and  allows  for  the 
simulation  of  V-I  curves  that  may  be  compared  to  existing  data.  The  model  takes  into 
account  cell  temperature,  pressure,  hydrogen  concentration  and  CO  concentration.  The 
original  model  is  two  dimensional,  considering  effects  both  through  the  cell  and  along  the 
channel.  However,  the  1-D  (through  the  cell)  version  is  modified  to  account  for  the 
effects  of  an  air  bleed.  This  includes  the  addition  of  O2  to  the  coverage  and  flux 
equations. 
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The  model  considers  the  PEM  fuel  cell  as  three  separate  regions  as  shown  in 


Figure  20. 


MEMBRANE- 

ANODE  GDL 

ANODE  ACR 

CATHODE 

Stephan-Maxwell 

Eq’s 

Flux  Eq’s 
Coverage  Eq’s 

Numerical  Fit 
based  on 
principles 

.1  \ 

x=0  X“(1gdl  ^acr) 


Figure  1  The  Springer  model  sees  the  PEMFC  as  an  anode  GDL,  active  catalyst 
region  (ACR)  and  a  membrane-cathode  region. 

The  basic  features  of  the  PEMFC  membrane-electrode  assembly  (MEA)  are  all 
present.  Because  the  focus  of  the  model  is  on  the  anode,  each  of  the  anode  regions  by  an 
individual  set  of  equations  while  the  membrane  and  cathode  portions  of  the  cell  are 
grouped  together  and  represented  by  a  single  equation. 

The  governing  equations  within  the  anode  GDL  region  are  given  by  the  standard 
Stephan-Maxwell  equations  for  multi-component  diffusion  as  shown  in  equation  (0.8). 

dz  8  PsDsij 

where  i  and  j  represent  the  various  anode  species,  xi  is  the  mole  fraction  of  each 
species,  and  Vs,  Ps  and  Ts  are  standard  molar  volume,  temperature  and  pressure.  The 
addition  of  oxygen  to  the  anode  feedstream  increases  the  number  of  species  in  the  model 
from  five  to  six  and  they  are  defined  as  follows:  1  =  H2,  2  =  CO,  3  =  N2,  4  -  CO2,  5-02 
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and  6  =  H2O.  In  the  case  where  H2=l  and  CO-2,  O, -1^/2?,  Oj-Jco^^F  and 
=  O4  =  <1>6  =  0 .  The  flux  of  oxygen  across  the  GDL  is  represented  by  (0.9): 


O5=[(j;:“U-J?0)  +  keffX02Paeh]/2F 


(0.9) 


The  multicomponent  pressure-diffusion  coefficients,  PsDsij,  are  calculated  by  the 
following  expression  (26): 


PsP^sij  ~  ^ 


^CELL 


/t  t 

V  C,1  -^CJ 


1 


(0.10) 


where  Te,i  and  pc.i  are  the  critical  temperatures  and  pressures  of  the  various 
components,  Mc,i  is  the  molecular  weight  of  each  component,  a  and  b  are  fitting 
parameters  and  e  is  the  Bruggeman  coefficient. 

The  reaction  mechanism  becomes  significantly  more  complex  when  oxygen  is 
present  anode  feedstream  (36,37,38,39,40,41,42,43,44,45,46,47,48,49,50): 


M  +  HjF  2(M-H) 

(0.11) 

M  +  CO€  (M-CO) 

(0.12) 

02  +  2M6  2(M-0) 

(0.13) 

(M-H)  +  (M-0)€  (M-OH)  +  M 

(0.14) 

(M  -  CO)  +  (M  -  OH)  -4  CO2  +  H"  +  e-  +  2M 

(0.15) 

(M-H)  +  (M-OH)€  H2O  +  2M 

(0.16) 

9 


(M-C0)  +  (M-0)-^C02+2M 


(0.17) 


(M-H)->H^+e'+M 


(0.18) 


In  the  case  of  the  Ruthenium  filter,  M  =  Ru.  As  described  in  the  previous 
chapters,  the  addition  of  oxygen  to  the  anode  feedstream  in  the  form  of  an  air  bleed 
results  in  the  formation  of  addition  hydroxyl  (M-OH)  groups.  It  is  the  extent  to  which 
these  groups  are  used  for  the  electro-oxidation  of  CO  that  determines  the  effectiveness  of 
the  air  bleed.  For  purposes  of  modeling  the  filter  region,  the  rigorous  incorporation  of  the 
entire  reaction  mechanism  described  in  equations  (0.1 1)  -  (0.18)  is  both  too  complex  and 
unnecessary.  This  would  result  in  four  site  coverage  expressions  (0h.0co»®o>®oH)» 
flux  equations  (for  Xh,Xco,Xo2),  along  with  the  appropriate  inert  species  equations  and 
voltage  relationships.  Like  Springer  et  al  have  done  with  the  H2+CO  reaction 
mechanism,  the  H2+CO+O2  mechanism  may  be  simplified  for  purposes  of  developing  a 
solvable  set  of  coverage  expressions  without  sacrificing  the  effects  of  the  air  bleed. 
Simplifications  of  the  mechanism  have  been  proposed  based  on  the  following  criteria: 

1.  The  rate  of  CO  oxidation  through  equation  (0.17)  is  negligible.  It  has  been 
shown  in  the  previous  chapters  that  O2  oxidizes  very  little  CO  catalytically. 

2.  The  model  must  take  into  account  that  O2  increases  the  CO  electrooxidation 
rate  above  the  level  found  when  only  CO  is  present  in  the  feedstream. 

3.  The  model  must  take  into  account  that  oxygen  oxidizes  both  Hads  and  COads, 
and  it  is  the  selectivity  of  this  oxidation  between  these  two  competing 
reactions  ((0.15)  and  (0.16))  that  determine  the  effectiveness  of  the  airbleed. 

4.  All  oxygen  is  consumed  before  it  penetrates  any  significant  distance  into  the 
ACR  As  shown  earlier  in  this  work,  a  sputtered  Ru  filter  of  less  than  100  nm 
was  shown  to  oxidize  200  ppm  CO  when  2  vol.  %  O2  was  present  in  the 
feedstream. 

5.  The  model  must  simplify  to  the  Springer  model  when  no  oxygen  is  present. 
Thus  a  single  set  of  equation  will  represent  all  feeds  (pure  H2,  dilute  H2, 
H2+CO,  H2+CO+O2)  entering  the  PEMFC  anode. 
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In  accordance  with  #1,  reaction  (0.17)  can  be  neglected  in  building  the 
mechanism.  Following  #2  and  #5  lead  to  the  assumption  that  there  is  an  exponential 
relationship  between  the  addition  of  O2  and  the  rate  of  CO  electrooxidation.  Thus,  when 
the  amount  of  oxygen  present  reduces  to  zero  (Xqj  ->0),  the  CO  electrooxidation  rate 

reduces  to  that  described  in  the  Springer  model  (jcow/02  jco)-  ^he  O2  selectivity 
described  in  #3  is  one  of  the  key  features  of  the  model.  Thus,  the  developed  expressions 
for  X02,  %  or  Gqh  must  describe  the  competitive  consumption  of  O2. 

Simplifying  the  reaction  mechanism  while  taking  into  account  the  effects  of 
0  is  the  most  desirable  as  it  is  the  amount  of  OHads  present  that  most  directly  affects 
the  rate  of  CO  oxidation.  However,  the  proposed  simplified  reaction  mechanism 
(described  in  Appendix  B)  is  still  too  complex  to  solve.  Similar  difficulties  arise  when  Gq 
is  used  to  account  for  increased  CO  oxidation  and  selective  O2  consumption  incorporated 
into  the  reaction  mechanism. 

Thus,  the  following  simplified  model  was  proposed; 


CO  +  M^ 

- ^(M-CO) 

bfc(^Co)^fc 

(0.19) 

H2+2M<= 

(0.20) 

io,+2(M 

H)-J^H20  +  2M 

(0.21) 

(M-H)- 

— +e“  +M 

(0.22) 

+  (M-CO)- 

-l^M  +  C02+2H"+2e- 

(0.23) 
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(0.24) 


k  -k®  exof— 1  (^-28) 

^f02  ”  ^f02  I  ^ 

This  version  satisfies  all  the  criteria  set  out  for  this  model,  while  maintaining  the 
simplicity  needed  achieve  stable  solutions.  Instead  of  using  0oh,  the  competitive 
consumption  of  O2  is  now  directly  related  to  oxygen  concentration  by  using  X02  as  shown 
in  (0.21)  and  (0.24).  The  model  then  maintains  the  simplicity  that  comes  from  having 
only  two  coverage  expressions  (for  0h  and  0co)  provides.  As  required  by  the  mechanism, 
oxygen  is  consumed  in  the  formation  of  water  and  the  consumption  of  M-CO  in  (0.21) 
and  (0.24).  However,  the  adsorption  and  reaction  of  O2  with  adsorbed  hydrogen 
described  in  (0.13),  (0.14)  and  (0.16)  is  now  represented  solely  by  (0.21). 


(0.29) 
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(0.30) 


2FRT,D 


PsVsT 


S^h  . 


dx. 


’  ^hloDL’ 


dz 


=  0 


filt/ACR 


^  ^CO  -ii  CO  w/02  )^ACR 
8Z^  Ir 


^DCO 


2FP,T,D 


a  ^S'-^CO  . 


9x 


PsVsT 


’  ^colfiit  ^COIgdL’  02 


CO 


=  0 


filt/ACR 


d\ 


02  _ 


dz^ 


(jcOw/02  Jco) 


^^fD2^02^a^h 


V  ^D02 


^^ACR 


_  2FPJsDo2 

PsVgT 


I  -  I  •  ^02 

>  ^02|fiu  “  ^02|gdl’  02 


=  0 


filt/ACR 


_  (Jh  jcow/02) 

dz^  crljiit 


(0.31) 

(0.32) 

(0.33) 

(0.34) 

(0.35) 


5V, 


dz 


—  Ya  filt/ACR  '^filler 


GDL/filt 


(0.36) 


where  the  local  and  overall  current  densities  for  H2  and  CO  are  given  by 
equations  (0.37)  -  (0.41): 


rv.+v,, 

] 

j02=ke02X02PA^C0Ce 

)  (0.37) 

Jco  “  J^ec^CO® 

(0.38) 

jh=2keAsmhj^  ""j 

(0.39) 

0 

11 

K) 

11 

! 

h— < 

0 

zr 

?  ^ 

'•  zr 

(0.40) 
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(0.41) 


JS=2F«.j=-I^^ 

where  jo2  is  the  CO  electrooxidation  current  resulting  from  O2  breaking  down 
into  OHads  and  then  reacting  electrochemically  with  COads-  Note  that  because  this  CO 
electrooxidation  pathway  is  distinct  from  the  pathway  involving  H2O,  it  is  represented  by 
a  different  rate  equation.  The  surface-area-to-volume  ratio,  ap, ,  is  set  at  1.0  x  10®  cm  ^ 
The  CO  that  is  oxidized  by  hydroxyl  group  formation  from  water  is  represented  by  the 
traditional  jco  described  in  Springer’s  model.  The  bulk  of  the  oxygen  is  consumed 
catalytically  with  adsorbed  hydrogen  and  is  represented  by  second  term  (k(Q2Xo2ii*a^h  )  i^ 
the  O2  and  H2  flux  equations.  However,  in  order  to  accurately  determine  the  amount  of 
O2  consumed  along  the  ACR  path  via  electrooxidation  and  hence  determine  the  total  rate 
of  O2  consumption,  jco  must  be  subtracted  from  j^o  w/02  • 

The  membrane-cathode  portion  of  the  model  is  represented  by  the  experimental 
relationship  developed  by  Kim  et  al  (35),  shown  in  equation  (0.42).  This  is  used  with 
equation  (0.43)  to  complete  the  model  in  1-D  form  by  relating  the  anode  cell  potential, 

Vanode,  to  the  cell  voltage,  Vceii- 

V.j=K,-K,J,„-K,ln{J,,)-K,exp(K5J„)  (0.42) 

V.„oDB=V.,(J„)-V„„  (0.43) 

where  the  constants  K1  -  K5  account  for  the  open  circuit  potential,  tafel  effects, 
IR  (current-resistance)  losses  and  gas  diffusion  limitations  occurring  in  the  membrane 
and  cathode  portions  of  the  fuel  cell.  Equations  (0.42)  and  (0.43)  reduce  the  entire 
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membrane-cathode  region  to  a  boundary  condition  which  is  used  in  conjunction  with  the 
ACR  flux  equations. 

The  model  equations  were  solved  using  MBAND(J).  Due  to  the  high  non¬ 
linearity  associated  with  the  effects  of  an  airbleed  addition,  4001  nodepomts  were 
required  for  the  ACR  region. 


Results  and  Discussion 

The  membrane-electrode  assemblies  tested  had  a  cross-sectional  area  of  50  cm  . 
The  cell  was  assembled  and  incubated  for  4-to-8  hours  at  ambient  pressure,  cell 
temperature  of  70°C,  an  anode  feed  of  hydrogen,  a  cathode  feed  of  air,  a  stoichiometnc 
ratio  ([actual  flow]/[stoichiometric  flow]  required  for  a  1.0  A/cm^  current)  of  1.5  at  the 
anode  and  2.0  at  the  cathode.  Single-cell  performance  curves  were  obtained  under  the 
conditions  set  in  Table  1. 


Table  1  Fuel  Cell  Test  Conditions. 


Pressure 

1  atm  (Anode  and  Cathode) 

Cell  Temoerature 

70°C _ _ _ 

Stiocheometric  Ratio  (at  1  A/cm'^) 

1.5  Anode  (Hydrogen) 

2.0  Cathode  (Air) 

Feedstreams 

Anode:  Hydrogen,  Reformate,  Reformate 
+  Air  Bleed,  H2  +  CO  +  Air  Bleed 

Cathode:  Air 

Dry  Reformate  Composition 

40%  H2 

20%  CO2 

50  ppm  CO 
balance  N2 
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Pt  Anode:  H2+CO  (Reformate)  Feed 

Figure  2  compares  the  model’s  predictions  to  data  obtained  at  Plug  Power  for  a 
PEMFC  containing  0.4  mg  Pt/cm^  in  the  anode  for  pure  H2  and  reformate  feeds.  First 
order  kinetics  are  used  here  to  describe  hydrogen  adsorption  and  oxidation. 

Several  parameters  were  modified  from  the  Springer  model  in  order  to  attain  the 
performance  curves  shown  in  Figure  2.  Differences  in  materials  used  (Nafion  type, 
catalyst  loading,  wt  %  catalyst  on  carbon,  wt  %  Nafion  in  the  catalyst  regions)  give  rise 
to  different  membrane-cathode  performance  characteristics  that  require  different 
constants  to  result  in  an  accurate  representation  of  the  cells  tested.  Of  more  importance  is 
that  the  values  for  kec,  and  ^(AGeo)  (1  x  10'^^  A  cm-^  1  x  lO'^^  atm  and  15.0  kJ 
mol'^  respectively)  used  in  this  work  are  significantly  lower  than  in  Springer  s  model  (1 
X  10‘®  A  cm■^  1.51  X  10'®  atm  and  20.0  kJ  mol'^  for  the  same  parameters).  PEMFCs 
tested  in  this  work  demonstrated  a  lower  CO  tolerance  than  those  examined  at  LANL. 
kec,  bfc  and  <5(AGco)  all  relate  to  rates  of  removal  of  CO  and  lowering  any  one  of  these 

three  parameters  results  in  the  lowering  of  the  CO  removal  rate.  Therefore  these  values 
were  adjusted  to  reflect  the  less  CO  tolerant  PEMFC  Pt  anode  catalyst  used  here. 

Figure  3  shows  the  increase  in  anode  potential  with  increasing  current  density. 
The  eell  potential  required  to  oxidize  pure  H2  ranges  from  0  to  0.05  V  as  the  eurrent 
density  increases  from  0  to  1.3  A/cml  As  the  amount  of  CO  in  the  anode  increases  from 
0  to  100  ppm,  the  anode  potential  required  to  achieve  similar  current  densities  increases 
(from  0.05  to  0.60  V  for  0.1  A/cm^).  Figure  3  and 
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Figure  4  show  the  correlation  between  the  CO  coverage  and  current  density.  As  the 
anode  potential  is  reached  at  which  the  CO  coverage  begins  to  decrease  (roughly  0.5  V 
for  50  ppm  CO)  the  cell  current  density  increases  significantly.  The  anode  potential  at 
which  this  occurs  shows  a  steady  increase  as  the  amount  of  CO  in  the  feed  increases. 
These  results  are  consistent  with  Springer’s  model  for  a  similar  system,  showing  that  the 
model  effectively  reduces  to  Springer’s  version  when  no  oxygen  is  present. 

Pt  Anode:  H2+CO+O2  (Reformate)  Feed 

Figure  5  compares  the  cell  performance  predicted  by  the  model  for  an  anode  feed 
containing  H2  +  CO  +  O2  to  actual  performance  data  of  a  PEMFC  MEA  containing  a  Pt 
anode.  Comparing  Figure  5  to  Figure  2  of  H2  +  CO  only  in  the  anode  feed  (no  oxygen  is 
present),  the  model  becomes  a  1-D  version  of  Springer’s  model.  This  is  as  it  should  be. 
Figure  5  also  shows  that  the  model  also  successfully  demonstrates  the  effect  of  adding  an 
airbleed  to  the  reformate  stream.  The  addition  of  0.5  vol.  %  O2  (relative  to  the  hydrogen 
concentration)  results  in  a  significant  increase  in  the  oxidation  of  all  CO  in  the  front 
portion  of  the  Pt  anode  (as  shown  in  Figure  6).  As  the  amount  of  oxygen  in  the  anode 
feed  is  increased  to  2.0%,  even  more  CO  is  oxidized.  Increasing  amounts  of  oxygen  in 
the  feedstream  results  in  increasing  local  CO  current  densities,  jco .  as  shown  in  Figure  7. 
All  CO  oxidized  through  an  electrochemical  reaction  with  O2  occurs  in  the  front  portion 
of  the  anode.  This  is  consistent  with  experimental  data  where  extremely  thin  Ru  filters  ( 
<  1.0  X  10'^  cm)  were  able  to  oxidize  up  to  200  ppm  CO  when  a  2%  air  bleed  is  added  to 
the  filter.  Figure  8  shows  the  effect  of  adding  O2  on  hydrogen  electrooxidation  rates.  In 
the  region  of  the  anode  where  oxygen  is  present  (shown  in  Figure  10),  H2  is  oxidized  at  a 
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higher  rate  due  to  the  reduction  of  CO  coverage.  For  higher  concentrations  of  02,  more 
CO  is  oxidized,  leaving  the  hydrogen  oxidation  rate  higher  even  after  all  02  is  consumed. 
This  is  because  lower  CO  concentrations  reduce  CO  coverage,  allowing  more  free  sites 
for  Ha  to  be  oxidized.  In  all  cases  the  bulk  of  the  O2  in  the  anode  feed  is  consumed 
through  catalytic  reaction  with  Ha  resulting  in  the  oxygen  concentration  reducing  to  zero 
well  before  the  end  of  the  ACR.  When  0.5  and  2.0  %  Oa  is  added  to  the  anode 
feedstream,  initial  CO  electrooxidation  rates  increase  by  a  factor  of  500  and  1000 
respectively.  However,  Oa  concentration  quickly  reduces  to  zero  along  the  length  of  the 
ACR,  returning  CO  oxidation  rates  to  the  same  order  of  magnitude  as  a  normal  reformate 
stream.  In  all  cases  where  Pt  anodes  were  used,  Oa  was  unable  to  oxidize  all  CO  present 
in  the  anode  feed  and  the  ACR  effectively  remained  poisoned. 

Figure  11  and  Figure  12  show  the  effects  of  voltage  on  the  amount  of  Ha  and  CO 
consumed  along  the  ACR.  As  cell  voltage  increases  (and  anode  voltage,  Vanode, 
decreases),  the  amount  of  CO  consumed  by  electrooxidation  with  Oa  and  HaO  decreases. 
As  the  cell  voltage  increases  from  0.25  to  0.80  V  (and  Vanode  decreases  from  0.37  to 
0.029  V),  the  amount  of  CO  oxidized  along  the  ACR  is  roughly  cut  in  half.  Both  means 
of  CO  electrooxidation  are  reduced  as  both  reactions  are  voltage  dependent.  Ha 
consumption  behaves  differently  from  that  of  CO  because  the  Ha  is  oxidized  catalytically 
by  Oa.  Figure  1 1  shows  that  the  amount  of  Ha  consumed  in  the  first  0.5  pm  of  the  ACR 
changes  very  little.  However,  the  amount  of  Ha  consumed  by  Oa  is  roughly  1%  in  all 
cases,  and  thus,  the  amount  of  Ha  consumed  overall  is  still  largely  dependent  upon  the 
operating  voltage  and  amount  of  CO  present. 
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The  catalytic  and  electrochemical  consumption  of  oxygen  is  represented  by 

and  ^  respectively,  where  kfo2  and  keo2  are  the 

preexponential  factors  for  those  expressions  given  in  A  cm  ^  atm  \  Decreasing  kf02 
results  in  an  increase  in  the  effect  of  oxygen  on  PEMFC  performance  to  the  point  where 
O2  is  able  oxidize  all  CO  present  in  the  feedstream  as  shown  Figure  13.  As  the  CO  is 
removed  near  the  front  of  the  ACR,  the  anode  operates  mostly  “CO  free”  resulting 
PEMFC  performance  nearly  equivalent  to  that  of  pure  H2  operation  as  shown  in  Figure 
14.  Figure  15  shows  that  the  H2  electrooxidation  rate  increases  dramatically  in  the  case 
where  CO  has  been  eliminated  from  the  feedstream.  As  the  H2  is  then  depleted  along  the 
length  of  the  electrode,  the  electrooxidation  rate  appropriately  reaches  a  maximum  and 
begins  to  decrease  to  the  Nemstian  effects  arising  from  hydrogen  dilution. 

Increasing  the  rate  of  O2  electrochemically  consumed,  keo2,  effects  PEMFC 
performance  similarly  to  decreasing  kfo2.  Figure  16  shows  that  as  keo2  is  increased, 
PEMFC  performance  increases. 

The  natural  Tafel  slope  of  the  electrochemical  oxidation  of  CO  by  O2,  bo2,  affects 
the  dependence  of  the  oxidation  rate  of  CO  on  anode  voltage,  Vg.  The  Tafel  slope  for  the 
electrochemical  oxidation  of  CO  by  O2  is  high  relative  to  bn  and  be-  This  makes  sense,  as 
the  electrooxidation  of  H2  and  CO  electrooxidation  by  H2O  are  far  more  voltage 
dependent.  Literature  has  shown  that  water  breaks  down  into  OHads  at  a  significant  rate 
only  at  high  anode  voltages  (-0.7  V).  H2  kinetics  are  well  documented  as  being  sensitive 
to  applied  voltage.  When  O2  is  added  to  the  feedstream,  the  performance  curves 
presented  in  Figure  5  show  that  significant  performance  increase  even  at  high  voltages  at 
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both  low  and  high  voltages  cell  voltages.  This  corresponds  to  significant  CO 
electrooxidation  occurring  at  both  high  and  low  anode  potentials.  As  shown  in  Figure  18, 
for  lower  values  of  bo2,  CO  oxidation  primarily  occurs  at  higher  anode  voltages 
(corresponding  to  lower  cell  voltages). 

It  is  believed  that  catalytic  oxidation  rate  between  hydrogen  and  oxygen  is  slower 
on  Ru  than  on  Pt.  This  gives  rise  to  the  increased  amount  of  CO  desorbed  and  resultant 
increase  in  PEMFC  performance  relative  to  pure  Pt  and  Pt:Ru  alloys. 

In  comparing  the  model  for  the  case  when  bo2  =  0.48  V,  kfo2  =  7  x  10  ^  A  cm  ^ 
atm'^  and  keo2  =  3  x  10“^  A  cm'^  atm"^  to  the  data  for  the  Ru  Filter  (shown  in  Figure 
19),  one  can  see  that  when  all  CO  is  oxidized  in  the  filter  region,  performance  does  not 
match  that  of  a  pure  H2  feed.  This  is  because  CO  is  present  as  part  of  a  reformate  feed 
that  contains  only  40  vol.  %  hydrogen  before  humidification.  The  model  successfully 
incorporates  the  effects  of  hydrogen  dilution  into  the  resulting  performance  curves. 
Increased  voltage  losses  through  the  dilution  effects  of  hydrogen  are  consistent  with 
literature  (32,51).  Some  performance  loss  can  also  be  contributed  to  hydrogen 
consumption  by  oxygen.  However,  this  amounts  to  a  maximum  of  4%  (when  2  vol  % 
oxygen  relative  to  hydrogen  is  present  in  the  feedstream).  This  performance  reduction  is 
far  less  than  that  resulting  from  CO  poisoning. 

Conclusions 

A  model  demonstrating  the  effects  of  adding  an  air  bleed  to  a  reformate  feed 
containing  CO  has  been  presented  here.  The  model  compares  well  with  data  gathered  for 
PEMFC  containing  a  Pt  anode  when  0.5  and  2.0  vol  %  oxygen  with  respect  to  hydrogen 
is  present  in  the  anode  feed.  All  oxygen  is  consumed  before  it  has  a  chance  to  diffuse 
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very  far  into  the  anode  filter,  which  is  also  consistent  with  existing  data, 
oxygen  is  present  in  the  anode  feed,  the  model  successfully  reverts  to  the  1 
Springer’s  model  for  an  anode  reformate  feed. 


When  no 
-D  form  of 
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Anode  overpotential  versus  current  density  for  a  variety  of  CO 
concentrations.  P  =  latm,  T  =  70  C. 
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CO  coverage  at  the  ACR  -  membrane  interface  versus  anode  overpotential 
for  a  variety  of  CO  concentrations.  P  =  latm,  T  =  70°C. 
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Cell  Voltage  (V) 


Figure  5  Comparison  of  single  cell  performance  data  to  the  developed  model  for 
the  Pt  +  Ru  filter  for  anode  feeds  of  pure  H2  (before  humidification),  40  % 
H2  +  50  ppm  CO  (before  humidification),  and  H2  +  50  ppm  CO  +  1%  O2 
(before  humidification).  P  =  latm,  T  =  70°C. 
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Local  CO  Current  Density,  A  cm' 


Distance  Along  ACR  (cm  x  1 0'*) 

igure  7  The  calculation  of  local  CO  current  densities  along  the  length  of  the  anode 
catalyst  region  for  an  anode  feeds  containing  CO  and  CO  +  O2.  P  =  latm, 

T-70°C,Vcell  =  0.4V. 
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The  calculation  of  H2  concentration  along  the  length  of  the  anode  catalyst 
region  when  CO  and  O2  are  added  to  the  feedstream.  P  =  latm,  T  -  70  C, 
=0.4  V. 


02  Mole  Fraction 


Distance  Along  ACR  (cm  x  10“*) 

Figure  10  The  calculation  of  O2  concentration  along  the  length  of  the  anode  catalyst 
region  for  an  anode  feed.  P  =  latm,  T  =  70  C,  Vcell  “  0.4  V. 
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ieure  1 1  The  calculation  of  CO  concentration  along  the  length  of  the  anode  catalyst 
region  for  a  fully  humidified  anode  stream.  P  =  latm,  T  =  70  C.  Initial 
feedstream  composition  (dry):  Xh2  =  0.40,  Xco  50  ppm,  and  l/o  O2 
(relative  to  H2).  bo2  =  0.48  V,  kfo2  =  7  x  lO"'  A  cm'^  atm'  and  keo2  =  3 
X  10"^  A  cm’^  atm”^. 


Mole  fraction  of  CO  in  ppm 


Figure  12  The  caleulation  of  CO  concentration  along  the  length  of  the  anode  catalyst 
region  for  a  fully  humidified  anode  stream.  P  =  latm,  T  =  70  C.  Initial 
feedstream  composition  (dry);  Xh2  ~  0.40,  Xco  ~  50  ppm,  and  1  ^  O2 
(relative  to  H2).  bo2  “  0.48  V,  kfo2  =  7x10  ^  A  cm  atm  and  keo2  ~  3 
X  10"^  A  cm"^  atm”^. 
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Figure  13  The  effect  of  varying  the  rate  constant,  kfo2,  for  the  catalytic  reaction 
between  H2  and  O2  on  the  amount  of  CO  oxidized  along  the  ACR.  P  = 
latm,  T  =  70°C.  Initial  feedstream  composition  (dry):  Xh2  =  0.40,  Xco  = 
50  ppm,  and  1%  O2  (relative  to  H2).  bo2  =  0.48  V,  keo2  =  3  x  10  ^  A  cm 


atm" 
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The  effect  of  varying  the  rate  constant,  kfo2,  for  the  catalytic  reachon 
between  Ha  and  O2  on  PEMFC  performance.  P  =  latm,  T  =  70°C.  Initial 
feedstream  composition  (dry):  Xh2  ”  0.40,  Xco  ~  50  ppm,  and  1  ^  O2 
(relative  to  Ha).  bo2  =  0.48  V,  keo2  =  3  x  10  ^  A  cm  atm  . 
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The  effect  of  varying  the  rate  constant,  kf02,  for  the  catalytic  reaction 
between  H2  and  O2  on  local  hydrogen  current  density,  jh,  along  the  length 
of  the  ACR.  P  =  latm,  T  =  70°C.  Initial  feedstream  composition  (dry): 
Xh2  =  0.40,  xco  =  50  ppm,  and  1%  O2  (relative  to  H2).  bo2  =  0.48  V,  keo2 
=  3  X  10"^  A  cm'^  atm-^ 
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The  effect  of  varying  the  rate  constant,  keo2,  for  the  electrochermcal 
oxidation  of  CO  by  O2  on  PEMFC  performance.  P  =  latm,  T  -  70  C 

Initial  feedstream  composition  (dry):  Xh2  -  0.40,  Xco  -  50  ppm,  and  0 

O2  (relative  to  H2).  bo2  =  0.48  V,  kf02  =  7  x  10  ^  A  cm  atm  . 
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The  effect  of  varying  the  rate  constant,  bo2)  for  the  electrochemical 
oxidation  of  CO  by  O2  on  PEMFC  performance.  P  =  latm,  T  =  70°C. 
Initial  feedstream  composition  (dry):  Xh2  =  0.40,  Xco  =  50  ppm,  and  1  /o 
O2  (relative  to  H2).  kfo2  =  7  x  10'^  A  cm'^  atm'^ ,  keo2  =  3x10'  A  cm 
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Figure  19  Comparison  of  the  performance  of  the  Pt  +  Ru  filter  to  the  model  for  a 
case  where  all  CO  is  electrooxidized  within  the  ACR.  P  =  latm,  T  - 
70°C  Initial  feedstream  composition  (dry):  Xh2  =  0-40,  Xco  =  50  ppm,  and 
P/o  O2  (relative  to  H2).  bo2  =  0.48  V,  kfo2  =  7  x  10'^  A  cm'^  atm"'  and 
ke02  =  3  X  10’^  A  cm"^  atm'\ 
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List  of  Symbols 

Table  below  lists  the  parameters,  their  definitions,  values  and  units.  (32,36,52) 


Symbol 

Definition 

Units 

Value  (Pt  ACR) 

Iacr 

ACR  thickness 

cm 

5x  10“' 

inrw 

GDL  thickness 

cm 

0.03 

keh 

pre-exponential  hydrogen  Tafel  electro¬ 
oxidation  rate 

A  cm"^ 

4.0 

kec 

pre-exponential  CO  Tafel  electro¬ 
oxidation  rate  times  2F  _  . 

A  cm‘^ 

3  X  10"^^ 

ke02 

pre-exponential  02  Tafel  electro¬ 
oxidation  rate  times  2F 

A  cm  atm 

3x  10^ 

kfh 

forward  hydrogen  adsorption  rate  times 
2F 

A  cm’'*  atm"' 

4000 

bfh 

back-to-forward  hydrogen  adsorption 
ratio  (1/Kea) 

atm 

0.5 

kfc 

forward  CO  adsorption  rate  times  2F 

A  cm''‘  atm'‘ 

10 

bfc 

back-to-forward  CO  adsorption  ratio 

(1/K.a) 

atm 

b£c° 

bfr  at  0CO  ”  0 

atm 

2.5  X  10‘^° 

kf02 

forward  oxygen  adsorption  rate  (convert 
to  mol  cm'^  atm"'  by  dividing  by  RT) 

s-' 

k£02° 

forward  oxygen  adsorption  rate  pre- 
exDonential  factor 

s-' 

7x  10' 

be 

Natural  Tafel  slope  for  CO  electro¬ 
oxidation  via  H2O 

V 

0.06 

K 

Natural  Tafel  slope  for  hydrogen  electro¬ 
oxidation 

V 

0.032 

bo2 

Natural  Tafel  slope  for  CO  electro¬ 
oxidation  via  H2O 

V 

0.48 

Dh 

Hydrogen  diffusion  coefficient  in 
catalyst  layer 

2  -1 
cm  s 

2  x  10-' 

Deo 

CO  diffusion  coefficient  in  catalyst  layer 

cm^ 

5.34  X  10'^ 

Do2 

O2  diffusion  coefficient  in  catalyst  layer 

cm  s 

Do,.D;.exp(|5:] 

Do  =  3.38x  lO'' 

Ed  =  16040  J 
mol'^  K'^ 

PsDsij 

pressure-diffusion  coefficient  of  species 
i  in  species  i  in  the  GDL 

atm  cm^  s""^ 

- - - - 

asvR 

surface  to  volume  ratio  of  Pt  in  the  ACR 

cm“^ 

IxlO" 

0h 

hydrogen  surface  coverage 

CO  surface  coverage 

veu 

0o 

oxygen  surface  coverage 

0nH 

OH  surface  coverage 

Xi 

mole  fraction  of  species  i 

Oi 

molar  flux  of  species  i 

mol  cm"^ 

Tc  i 

Critical  temperature  of  species  i 

K 

Pc,l 

Critical  pressure  of  species  i 

atm 
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Mi 

Molecular  weight  of  species  i 

g  moV^ 

Va 

Local  anode  overpotential 

V 

Vanode 

Anode  overpotential  at  the  ACR- 
membrane  interface 

V 

Vn 

Nemstian  potential:  =5iin(xhP3) 

nF 

V 

Vs 

Standard  molar  volume 

cm^  moL^ 

22414 

Ps 

Standard  pressure 

atm 

1.0 

Ts 

Standard  temperature 

K 

273.15 

Hvdroeen  diffusion  current  density 

A  cm'^ 

CO  diffusion  current  density 

A  cm’^ 

Oxveen  diffusion  current  density 

A  cm‘“ 

i  h 

_ JLS. - — — — - ; - — ^ 

local  hvdroeen  current  density 

A  cm'^ 

local  CO  current  density 

A  cm"^ 

TtOt 

total  hydrogen  current  density 

A  cm‘^ 

TtOt 

**co 

total  hydrogen  current  density 

A  cm‘^ 

J.O. 

total  current  density 

A  cm‘^ 

F 

Faradays  constant 

C  mol'^ 

96487 

R 

Universal  gas  constant 

J  mof'  K-^ 

8.315 

cm^  atm  mol'^ 

K-' 

82.06 

5(AGc„) 

Variation  of  free  energy  of  adsorption  of 
CO  between  zero  and  full  coverage  of 
CO 

J  mol'^ 

20,000 

Eo2 

The  activation  energy  for  the  reaction  of 
oxygen  with  adsorbed  hydrogen 

cal  moL^ 

11,400 

1  n 

exchange  current  density 

A  cm"^ 

_Ji2 - 

Z 

^ - - - - - * -  ' 

distance  through  the  anode  _ _ 

cm 

Xi 

mole  fraction  of  species  i 

S; 

stoichiometric  coefficient  of  species  i 

*^1 

Oi 

surface  coverages  of  species  I  (CO,H2) 

0h 

equilibrium  surface  coverage 

n 

electrochemical  overpotential 

V 

a 

conductivity  of  the  electrically 

conductive  phase 

mho  cm‘‘ 

ref 

Cl44- 

H2  reference  concentration 

mol/cm*^ 

T 

Temperature 

c 

eff 

effective  (accounting  for  a  pourous 
medium) 

ACR 

. — z -  - - - 

Active  Catalyst  Region 

GDL 

Gas  Diffusion  Layer 

MEM 

membrane 

BEST  AVAILABLE  COPY 
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Appendix  A:  Derivation  of  Vogel’s  anode  current  density 
expression 

The  derivation  of  the  current  generated  by  the  dissociation  of  H2  in  terms  of 
surface  coverage  of  the  Pt  catalyst  done  by  Vogel  is  presented  here.  The  surface 
coverage  of  the  adsorbed  H2  is  governed  by  the  electrode  potential  given  by  the  Nemst 

equation  (23,31)- 

£  =  In^  (0.44) 

F  ^MH 


Ev°  is  the  open  circuit  electrode  potential,  s  is  the  actual  electrode  potential  and 
6m,  6mh  are  the  empty  sites  and  sites  covered  by  adsorbed  hydrogen,  respectively.  At 
equilibrium,  this  expression  becomes: 


(0.45) 


Where  6h  is  the  number  of  sites  occupied  by  adsorbed  hydrogen  at  steady  state. 


Using  the  relationships  shown  in  (0.46)  and  (0.47) 

?]  =  s-E 

(0.46) 

^MH  ”1 

(0.47) 

0mh  can  be  solved  for  using  equations  (0.44)  and  (0.45)  to  give  (0.48): 


^MH  ~ 


Ou 


(0.48) 
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This  expression  can  then  be  substituted  into  the  rate  expression  for  H2 
+2e.  The  current  density  generated  is  related  to  the  rate  of  reaction  (23)  which  is  derived 
from  the  Tafel-Volmer  reaction  in  equation  (0.49): 

Jhi  =r„„nF=  nF(k„c„;«^ -K^m) 

kat  is  the  rate  of  adsorption  of  H2  onto  the  catalyst  and  kct  is  that  constant  for 
hydrogen  desorption.  The  number  of  electrons  generated,  n,  is  2  and  ch2  is  the 
concentration  of  H2  along  the  ACR.  At  equilibrium  no  current  is  generated  0=0)  and  the 
equation  (0.50)  is  obtained  for  the  exchange  current  density  (23): 

jo=2Fk,CH,^^=2FkA  (0.50) 


Noting  that  0m  =1-0mh,  the  current  density,  j,  can  be  solved  for  using  (0.49)  aand 
(0.50),  resulting  in  equation  (0.51): 


j  =  jo 


_ 1_ 


eF,/RT 


^  F;7/RT 
C« 

V'^H2 


A 

1 

) 


(0.51) 
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Appendix  B:  Alternative  H2/CO/O2  Mechanism 

Incorporating  6ok 

The  following  is  the  simplified  mechanism  that  uses  0oh  to  determine  the 
selectivity  and  effectiveness  of  O2  usage. 


CO  +  M:t 


bfc(^Co)^fc 


H2+2M4 


^^lh(^co) 


02+2(M-H)-^ 


i(M-CO) 

±2(M-H) 

->2(M-OH) 


( M  -  H)  +  (M  -  OH)  — H^O  +  2M 
(M  -  H)— H'' +  e"  +  M 
(M  -  OH)  +  (M  -  CO)  >2M  +  CO2  +  2H"  +  e' 


(0.52) 

(0.53) 

(0.54) 

(0.55) 

(0.56) 

(0.57) 


(V.+Vn 


p(9co=kf,XcoPA(l-^co-^h-^oH)-bfcMco-kece^“’'°''^e'^  =0 


p^h  =kn,x,,PA(l-^co  “^h  "^oh)  k^H^oH^h 


-  2kj.,,^h  sinh 


V  +V, 


V  7 


(0.59) 


=  0 


P  0OH  =  kf^XgjPA  (1  ~  ^CO  ~^h~  ^OH  )^h  kjQH^h^OH  kfQH^OH^h 

X  »h  J^o 


bfc  =b?eexp 


0, 


<?(AGco) 


CO 


RT 


(0.60) 


(0.61) 
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(0.62) 


kro2=kro2exp[^^^ 

bfDH  =  broH  exp  (a,^oH ) 


Like  in  the  case  where  only  CO  and  H2  are  present  in  the  feedstream,  first  order 
hydrogen  adsorption  kinetics  is  used  (n  =  1)  for  simplification  of  the  model. 

The  critical  element  of  this  reaction  mechanism  is  the  usage  of  the  hydroxyl 
group,  M-OH.  The  production  of  M-OH  beyond  what  occurs  in  a  feedstream  where  only 
H2  and  CO  are  present  is  the  result  of  the  reaction  shown  in  (0.53)  -  which  is  a 
simplification  of  O2  adsorption  and  reaction  with  M-H.  As  O2  adsorption  is  the  rate- 
limiting  step,  that  reaction  coefficient  is  used  for  kf02  as  shown  in  (0.62).  Equations 
(0.13)  and  (0.14)  have  been  combined,  and  since  (0.14)  occurs  very  quickly,  the  reaction 
rate  is  that  of  (0.13).  In  order  for  the  model  to  simplify  to  the  Springer  model,  the  reverse 
reaction  through  which  water  dissociates  in  (0.16)  is  neglected.  This  is  acceptable, 
because  this  phenomenon  is  taken  into  account  by  the  base  Springer  model.  CO 
oxidation  through  equation  (0.17)  has  also  been  neglected  in  this  reaction  mechanism. 

As  it  is  the  additional  M-OH  that  gives  rise  to  the  added  electrooxidation  of  CO, 
this  has  been  incorporated  into  j^o  w/02  shown  in  (0.64). 


Jco  w/02 


(“l^OH  ) 


(0.64) 


This  exponential  term,  — » las  ^  1  •  This  model  assumes  that  0oh~O.O 

when  no  oxygen  is  present  in  the  feedstream. 
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Appendix  C:  Springer’s  Model 

They  have  developed  separate  expressions  for  hydrogen  and  CO  electrooxidation 
occurring  on  from  the  surface  of  the  Pt  anode: 


jh  =2k,,^hSinh 


(0.65) 


Jco  —  2^ec^CO® 


(0.66) 


where  jh  and  jco  are  the  local  electro-oxidation  current  densities  for  H2  and  CO,  0h 
and  0CO  are  the  surface  coverages  of  H2  and  CO,  respectively,  Va  and  Vn  are  the  anode 
and  Nemst  potentials,  keh  and  kec  are  the  electrooxidation  rate  constants  for  H2  and  CO, 
and  be  and  bh  are  the  Tafel  slopes  for  H2  and  CO  electro-oxidation,  respectively.  While 
the  expression  for  jh  doesn’t  directly  contain  0co,  they  are  related  through  the  coverage 
expressions  and  are  based  on  the  simplified  reaction  mechanism  developed  by  Springer 
for  an  anode  feed  containing  H2  and  CO  as  shown  below  (32): 


CO  +  M,  ‘‘^(M-CO) 

Ofcl^CO 

(0.67) 

(0.68) 

(M  H)— ^H^+e"  +  M 

(0.69) 

H^O  +  (M  -  CO)— ^ M  -h  CO2  +  2H"  +  e- 

(0.70) 

p0co=kf,XcoPA(l-^co-^h)-bfckfAo-kece^  -0 

(0.71) 
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pk=k^x,?^{l-e^o-^J-  KKK  -  ^KA  sinh 


V  bh  y 


=  0 


(0.72) 


bfc=bfcoexp 


0. 


S[AGqq) 


CO 


RT 


kf,=k?,exp 


^(AEJ 


RT 


1-e’ 


P^co  ^ 
1-^co. 


(0.73) 


(0.74) 


where  n  =1  or  2  (1*'  or  2"‘‘  order),  kfc,  k^,  bfckfc,  and  bfhkfh  are  the  forward  and 
reverse  rate  constants  and  are  functions  of  ^co  shown  in  equations  (0.73)  and 

(0.74).  In  this  work  kfc  is  kept  constant.  PA  is  cell  pressure,  Xco  and  Xh  are  mole  fraction 
of  H2  and  CO  respectively,  bh  and  be  are  the  natural  Tafel  slopes  for  H2  and  CO  electro¬ 
oxidation  respectively. 

This  model  is  both  rigorous  in  the  description  of  the  PEMFC  anode  catalyst 
region  and  allows  for  the  simulation  of  V-I  curves  that  may  be  compared  to  existing  data. 
The  model  takes  into  account  cell  temperature,  pressure,  hydrogen  concentration  and  CO 
concentration.  The  original  model  is  two  dimensional,  considering  effects  both  through 
the  cell  and  along  the  channel.  However,  only  the  1-D  (through  the  cell)  version  is  used 
described  here.  Thus,  the  Springer  model  chosen  as  a  basis  to  develop  a  model  for 
describing  the  Ru  Filter. 

These  kinetic  expressions  are  integrated  into  the  overall  model.  The  model 
considers  the  PEM  fuel  cell  as  three  separate  regions  as  shown  in  Figure  20. 
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Figure  20  The  Springer  model  sees  the  PEMFC  as  an  anode  GDL,  active  catalyst 
region  (ACR)  and  a  membrane-cathode  region. 

The  flux  equations  for  the  anode  active  catalyst  region  (ACR)  are  given  in 
equations  (0.75)  -  (0.80): 


=  (0.75) 

loH 


^  2FP3TSD, .  I  _  I  . 
^Dh  PVT  ’  ^hlcDL’ 


ax. 


dz 


=  0 


ACR /MEM 


(0.76) 


^  ^co  _  Jco 

dz^  I 


DCO 


(0.77) 


T  2FPJsDco  .  I  ^1  .  ^ 

■‘dco  pvt  ’  coIgdl’ 


CO 


’S' 


=  0 


ACR/MEM 


(0.78) 


dX  (Jh+Jco) 


dz^  <j\ 


ACR 


(0.79) 


av. 


dz 


aCR/MEM  X\N0DE 


GDL /ACR 


(0.80) 


The  membrane-cathode  portion  of  the  model  is  represented  by  the  experimental 
relationship  developed  by  Kim  et  al  (35)  as  shown  in  equation  (0.81).  The  membrane- 
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cathode  potential,  Vgtd,  is  calculated  based  on  overall  current  density,  Jtot,  and  then  used 


to  relate  the  anode  potential,  Vanode,  to  the  cell  potential,  Vceii,  in  (0.82). 


V.,  =  K,  -  K,  In  (J„, )  -  K.  exp(K,J„ ) 


(0.81) 

(0.82) 


The  constants  used  in  this  work  to  calculate  Vstd  are  given  in  Table  2  below: 
Table  2:  Values  for  the  Membrane-Cathode  Relationship 


Symbol 

Value 

units 

K1 

0.795 

V 

K2 

0.1895 

V  cm^  A”^ 

K3 

0.0176 

V 

K4 

4.54  X  10''‘ 

V 

K5 

4.59 

cm^ 

Thus,  cell  performance  curves  can  then  be  calculated  that  take  into  account  the 
losses  from  the  various  regions  as  shown  in  Figure  21. 
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Theoretical 

Maximum 

Potential 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0  1.1  1.2 


Current  Density,  A/cm^ 


Figure  21  The  various  PEMFC  performance  losses  taken  into  consideration  by  the 
Springer  model  are  represented  above.  This  model  focuses  primarily  on 
the  anode  activation  losses  and  losses  due  to  the  addition  of  CO  to  the 
anode  feed  stream. 
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Appendix  D:  Model  Source  Code 


On-Board  Fuel  Cell  Power  Plants 


Final  Report 


Appendix  3  -  Model  Help  Files 

A  Methanol  Reformer  VTB  Model  dated  1 1/15/00 
“Simulink”  Interface  VTB  Models  dated  6/23/01 
“Fuel  Cell  12.15_010716”  VTB  model  dated  7/16/01 
Metal  Hydride  Bed  VTB  Model  dated  3/18/02 
i  VanePumpli  VTB  Model  dated  6/10/02 
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A  Methanol  Reformer  VTB  Model 

Author:  Ramakrishna  Gundala 
Date:  11/15/00 
Model  name:  Reformer 
DLL  name:  CReformer 
Version  number:  1.0 

Report  errors  or  changes  to:  vtb@engr.sc.edu 

Pictorial  Representation 


Brief  Description  of  Model 


This  is  a  CSTR  model  of  the  reformer.  The  temperature  of  the  catalyst  bed,  methanol  conversion 
and  the  pressure  inside  the  reactor  change  with  time  but  do  not  vary  within  the  reformer.  External 
heat  input  to  the  reformer  is  also  included  to  control  the  reformer  bed  temperature  at  the  desired 

value. 


Validity  Range  and  Limitations 


The  reformer  is  designed  to  operate  for  steam  to  methanol  ratio  of  0.5  to  2.0, 150®C  to  266®  C 
temperature  and  2  to  3  bar  pressure. 


Terminals  and  Connectivity  Information 


This  model  has  the  tollowmg  terminals. 

Terminal  Designation 

Description 

A 

Input  for  steam  and  methanol  mixture 

B 

Output  for  hydrogen,  CO  and  C02  and  unused 
methanol 

Parameters  and  Output  Variables 

This  is  a  complete  list  of  user  accessible  parameters.  SI  units  are  used  throughout. 


Parameter 

Description 

Default 

Value 

Units 

Pow^ 

Required  Power 

1000 

KW 

Td 

Operating  Temperature 

266 

°C 

Pd 

Operating  Pressure 

3 

atm. 

SIM 

Steam  to  Methanol  Ratio 

1.1 

- 

Xd 

Design  Conversion  of  methanol 

0.95 

To 

Initial  temperature  of  reformer 
bed 

266 

®c 

- m - 

UA, 

Heat  transfer  conductance 

3.99 

Watts/ "C 

Fo 

Methanol  feed  rate  for  a  IMW 
system 

3.30 

moles/s 

L - - - - - - - 

The  following  variables  are  eomputed  within  the  model  and 

I  are  user  accessible. 

Variable  Name 

Units 

Operating  Temperature 

oc 

Ky  produced 

moles 

Weight  of  the  bed 

Kg _ _ _ 

Fraction  of  CO  in  the  reformer 
exit  stream 

Fraction  of  methanol  converted 

- 

Power  delivered 

Watts 

Heat  requirement  of  the  reformer 

KJ/s 

Pressure  at  the  reformer  outlet 

atm. 

Total  methanol  into  the  reformer 

moles/s 

Assumptions  in  Model  Derivation 

The  reformer  is  assumed  to  be  a  CSTR  model  i.e.  the  temperature,  concentrations  variations 
inside  the  bed  are  neglected.  The  heat  input  to  the  reformer  is  assumed  to  be  provided  by  a 
proportional  controller  that  is  within  the  reformer  and  is  not  shown  explicitly. 


Mathematical  Description  of  Model 


Methanol  (methyl  alcohol)  is  catalytically  steam  reformed  via  the  overall  reaction: 

CH^OH  +  H2O - >  CO2  +  Si/j  (reforming)  [  1  ] 


Some  portion  of  methanol  also  decomposes  directly  to  CO  via  the  reaction: 

CH^OH  ^CO  +  lHj  (decomposition) 

The  rates  r,  and  of  the  above  reactions  can  be  described  adequately  as 
^1  = 

ri=-K 


[2] 

[3] 

[4] 


where  stands  for  the  methanol  concentration  and  and  k2  are  the  rate  constants  which  are  functions 
of  temperature  and  pressure  as  given  below. 


h  = 


(4  +  5ilog(— )) 

- ^^^exp(- 


i?(r^  +  273) 


) 


[5] 


A2  exp(- 


k^  =■ 


R{Td+21l) 


) 


At  constant  temperature  and  pressure,  the  concentration  of  methanol  as  a 
X  can  be  expressed  as 


function  of  the 


[6] 


conversion  ratio 


Cat  W  =  0”^) 

1  +  — +  2x 
M 

Combining  equations  (3)  and  (7)  reaction  rate  f  |  can  be  written  as 


[7] 


1  +  - 


'*iW  =  0-^) 


M 


^  ■' 
1  +  — +  2x 
M 


^.(0) 


[8] 


The  weight  of  the  reformer  can  be  easily  obtained  as  a  function  of  the  methanol  conversion  (X )  after  doing 
some  simple  algebra  and  is  given  as 


U  2Y  Y 


[9] 


Where,  for  simplicity,  we  use  17  =  1  +  5  /  M,  r^  —  r,(0)andr2  —r2(0).Y  U(r^  +r2  )  and 

Z  =  f/r,''-2r2®. 

The  total  output  flow  equals  {l  +  S  /  M  +  2x)F  and  the  H2  flow  3xF  —  ^2  fV(x) 

Rate  of  heat  generation  within  the  reformer  due  to  the  reaction  is  given  by 

Q,e.=-i^F-k2W{x))m,,„, -k2W(x)AH„„2  [10] 

Where  AH,.„,|  is  the  enthalpy  change  of  the  reaction  1  and  AH^j„2 *1^®  enthalpy  change  of  reaction  2. 
The  heat  balance  gives 

(T-T,)-±F,._,Cp  (T-T^-Q,.  tm 

1=1  out_i 

Where  out  _  i  =  (CH^OH,  H2O,  H2 ,  CO2 ,  CO) ,  in  _  i  =  (CH2OH,  H2O) ,  Q?,  is  the  heat  capacity 

of  component  land  is  the  reference  temperature. 

A  proportional  controller  is  used  to  regulate  the  reformer  temperature  according  to  the  law, 

Q^=va,(t,-t)  [12] 

Example  of  Model  Use 

In  figure  2,  the  reformer  is  connected  to  a  constant  methanol  source  at  the  input  and  a  H2 
tank  at  the  output. 
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Figure  3.  The  fractional  conversion  of  methanol  in  the  reformer  bed  changes  as  a  function  of 
temperature.  Final  value  is  0.95 


temperature,  266  C  with  time. 
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Figure  5.  H2  obtained  from  the  reformer  exit  in  moles/s.  Steady  state  value  is  0.187  moles  /  s 
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Figure  6.  The  heat  input  into  the  reformer.  The  steady  state  value  is  81 .69 
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Model  Validation 


To  validate  the  model,  simulations  were  performed  under  different  conditions  and  the 
results  obtained  conformed  to  the  behavior  of  the  reformer  model  along  the  expected 
lines  For  example,  the  starting  temperature  is  assumed  to  be  30  0  C  and  the  reformer 
takes  about  10  minutes  to  reach  the  steady  state  value.  The  simulations  gave  good  results 
for  a  range  of  design  temperatures,  power  requirements,  and  pressures. 
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“Simulink”  Interface  VTB  Models 

Author:  Teems  Lovett,  Antonello  Monti 
Date:  June  23, 2001 
Model  name:  5in_5out  Simulink 
DLL  name:  simulink55.dll 
Version  number:  2.0 

Report  errors  or  changes  to:  lovett@engr.sc.edu  or  monti@engr.sc.edu 

Installation  of  Simulink  Models 

In  order  to  install  and  use  the  Matlab-Simulink  interface  models  there  is  a  set  of  instructions  that  must  be 
followed.  These  instructions  include  the  setup  to  use  Matlab  remotely  from  VTB.  Once  these  steps  ^e 
accomplished  then  adding  more  Matlab-based  interface  models  is  simple  and  does  not  require  repeating 
these  steps. 

Step  1 .  Copy  the  following  DLL’s  from  the  Matlab  installation  folder  on  the  host  computer  into 
the  installation  folder  ofthe  Virtual  Test  Bed.  The  DLL’s  must  be  copied  into  the 
directory  with  the  VTB  and  Model  Manager  executables.  The  DLL’s  will  be  located  in 
the  bin  folder  ofthe  Matlab  installation  path. 
libeng.dll 
libmat.dll 
Iibmi.dll 
libmx.dil 
libut.dll 
libmatlbmx.dll 

Step  2.  Download  the  calcsim.p  from  the  model  library  web  page  and  place  this  file  into  the  same 
directory  as  step  #1. 

Step  3.  Start  Matlab  and  set  the  same  directory  as  used  in  steps  #1  and  #2  as  a  permanent  path  tor 
Matlab. 

Pictorial  Representation  of  Model 
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Brief  Description  of  Model 

VTB-Simulink  interface  uses  and  requires  a  registered  copy  of  Matlab  to  be  loaded  onto  the  host  computer. 
The  input  ports  are  modeled  as  infinite  impedance  loads.  The  output  ports  are  modeled  as  ideal  voltage 
sources  connected  to  ground.  The  VTB>Simulink  interface  requires  a  Simulink  .mdl  file  to  be  located  in  an 
active  or  permanent  path  of  Matlab.  For  this  specific  model  the  file  is  Test55.mdl.  The  voltages  at  the 
inputs  ofthe  VTB  model  are  sent  into  the  input  ports  of  the  Simulink  file.  The  magnitude  ofthe  ideal 
output  voltage  sources  is  with  respect  to  the  outputs  received  from  the  Simulink  file.  Every  time-step  the 


Simulink  file  is  executed  from  Matlab  and  solved  for  the  time  difference  between  the  previous  time  value 
and  the  present  time  value. 


Model  Validity  Range  and  Limitations 

The  VTB-Simulink  models  are  restricted  to  the  necessity  of  certain  Matlab  DLLs  being  located  in  the 
executable  directory  of  the  schematic  editor.  Also  the  loading  of  the  VTB-Simulink  models  the 

Matlab  DLLs  being  either  in  the  working  directory  of  Model  Manager  or  the  location  of  the  DLL  to  e 
added  to  the  database.  A  Simulink  .mdl  with  the  correct  name  and  number  of  inputs  and  output  is  required 
to  be  in  a  working  directory  of  Matlab.  .  ,  .  .  j  .u  I 

The  validity  range  of  the  VTB-Simulink  models  depends  on  the  meeting  of  requirements  and  the  validity  ot 
the  Simulink  file  being  interfaced  into. 


List  of  Model  Pins  with  Connectivity  Information 


Pin  Designation 

Description 

N-input 

Infinite  impedance  input  nodes 

N-output  — 

Ideal  voltage  sources 

List  of  Parameters  and  Output  Variables 

This  is  a  complete  list  of  all  parameters  of  the  model.  All  models  use  SI  units. 


Parameter  Name 

Description 

Default  Value 

Units 

Iterations-skipped 

Number  of  VTB  iterations  to  be 
preformed  before  executing  one 
step  of  the  Simulink  model. 

1 

NA 

This  is  a  list  of  output  variables. 


Variable  Name 

Description 

Units 

None 

None 

NA 

Assumptions  in  Model  Derivation 

Infinite  impedance  inputs  and  ideal  voltage  source  outputs  connected  to  ground. 


Mathematical  Description  of  Model 

The  heart  of  the  communication  protocol  with  the  Simulink  solver  is  the  Matlab  engine.  The  Matlab 
allows  data  exchange  between  Matlab  and  C/C++  custom  software.  The  interface  is  realized  by  a  dedicated 
class  inside  the  VTB  architecture.  This  class  manages  the  following  procedure: 

•  At  simulation  start  the  VTB-class  calls  the  Matlab  engine 

•  The  VTB  calls  a  null-time  execution  of  the  target  Simulink  model.  The  result  of  this  step  is  the  calculation  of  the 
Simulink  model  size  in  terms  of  state  variables,  input  and  output  number. 

•  At  every  step  the  VTB  calls  a  step  execution  of  the  Simulink  model  giving  as  input  the  required  data  and  receiving 
back  the  data  for  the  next  step 

•  The  output  is  inserted  in  the  input  vector  of  the  RCM 

This  data  exchange  is  implemented  by  the  authors  via  a  dynamically  linked  library  object  and  new  Matlab 
commands  that  run  through  the  Matlab  engine  a  simulation  step  of  the  target  Simulink  model.  A  great 
advantage  of  this  solution  from  the  user  viewpoint  is  the  possibility  to  interact  (hrectly  during  the 
simulation  both  with  the  VTB  Schematic  Editor  and  the  Simulink  user  interface.  Parameters  can  be 
changed  on  the  fly  in  both  the  environments  thereby  speeding  up  any  testing  activity. 


Example  of  Model  Use 


A  Vector  Control  of  permanent  magnet  linear  motor 

This  example  is  a  vector  control  for  a  permanent  magnet  linear  motor.  As  per  the  theory  described 
previously,  the  power  plant  is  described  using  VTB  native  models  while  the  control  is  described  by  means 
of  a  Simulink  model.  In  Figure  1  the  VTB  model  is  introduced. 


Figure  1 :  The  VTB  model  used  as  example 

The  “Matlab  Model”  block  represents  the  n-port  (5-port)  interconnection  previously  introduced. 


Figure  The  correspondent  Simulink  file  (test55.mdl) 

Note  that  if  you  are  not  using  all  the  5  inputs  or  outputs  as  real  variables  for  the  Simulink  program,  you  still 
have  to  declare  5  input  and  5  output  ports  in  the  Simulink  schema. 


Model  Validation 

Figure  3  shows  a  3D  visualization  obtained  at  simulation  time.  In  this  case  the  speed  is  controlled  to  follow 
reLence  signal  that  is  determined  by  means  of  the  programmable  source.  The  transient  dynamics  are 


apparent. 


References 

MATLAB  online  help  documents 


“Fuel  Cell  12.15.01071 6”  VTB  Model 

Author:  Mark  J.  Blackwelder 

Date:07/ 16/01 

Model  name:  FC_12 

DLL  name:  Fuel_Cell_12_010716.dll 

Version  number:  12.15 

Report  errors  or  changes  to:  blackwel@engr.sc.edu 


Pictorial  Representation 


Brief  Description 

The  fuel  cell  model  imitates  the  operation  of  a  polymer  electrolyte  membrane  fuel  cell.  This  includes  the 
electrochemical  voltage-current  relationship,  heat  production,  fuel  and  air  consumption,  and  temperature 
and  pressure  dependence. 

Validity  Range  and  Limitations 

The  model  is  valid  only  for  current  sourced.  If  current  is  forced  against  the  electrochemical  potential,  then 
the  results  are  invalid.  For  example,  this  model  does  not  electrolyze  water.  The  results  are  also  invalid  for 
pressure  and  temperature  ranges  outside  of  the  experimental  data  used  to  fit  the  empirical  parameters.  For 
the  default  case,  pressure  outside  the  7.4-30PSIG  range,  and  temperature  outside  of  the  -10  to  60  degree 
Celsius  range,  are  invalid.  Also,  note  that  air  is  consumed,  not  just  oxygen.  When  mixed  material  streams 
are  integrated’ into  VTB,  this  will  be  modified  such  that  oxygen  is  consumed,  and  the  characteristics  are 
dependent  on  the  partial  pressure  of  oxygen.  The  thermal  transfer  by  mass  streams  is  considered 
negligible. 


Terminals  and  Connectivity  Information 


inis  moaei  nas  inc  icmunaio. 

Terminal  Designation 

Description 

+ 

Positive  electrochemical  potential  terminal 

Negative  electrochemical  potential  terminal 

H2 

Hydrogen  fuel  connection 

Through  Variable 

Moles/Second 

Across  Variable 

Pascals 

Air 

Air  mass  stream  connection 

Through  Variable 

Kilograms/Second 

Thermal  generation  and  transfer  terminal 

Through  Variable 

Watts 

Across  Variable 

Degrees  Celsius 

Parameters  and  Output  Variables 

This  is  a  complete  list  of  user  accessible  parameters.  SI  units  are  used  throughout. 


Parameter  Name 

Description 

Default  Value 

Units 

Open  circuit  potential 
npr  cell 

Depends  on  fuel  and 
oxidant  used 

1000 

millivolts 

Exponential  Squared 
pressure  dependent 
term  (no) 

The  portion  of  the 
exponent  that  is 
dependent  on  the  square 
of  the  pressure 

.2.843678952E-13 

cm^milliamp/pascal" 

Exponential  Linear 
pressure  dependent 
term  (nO 

The  portion  of  the 
exponent  that  is  linearly 
dependent  on  the 
pressure 

1.162035969E-7 

cmVmilliamp/pascal 

Exponential  constant 

The  constant  portion  of 
the  exponent 

-6.390737001E-3 

cm  /milliamp 

Logarithmic  Squared 
temperature  dependent 
coefficient  (bo) 

The  portion  of  the 
logarithm  that  is 
dependent  on  the  square 
of  the  pressure 

3.091E-3 

millivolts/cell/degree 

Kelvin^ 

Logarithmic  Linear 
temperature  dependent 
coefficient  (bi) 

The  portion  of  the 
logarithm  that  is 
linearly  dependent  on 
the  pressure 

-2.0524 

millivolts/cell/degree 

Kelvin 

Logarithmic  Constant 
coefficient  (b?) 

The  constant  portion  of 
the  logarithm 

398.57 

millivolts/cell 

Resistive  Squared 
temperature  dependent 
coefficient  (ro) 

The  portion  of  the 
resistance  that  is 
dependent  on  the  square 
of  the  temperature 

1.5168E-5 

Ohms/degree 

Kelvin^*cmVcell 

Resistive  Linear 
temperature  dependent 
coefficient  (rO 

The  portion  of  the 
resistance  that  is 
linearly  dependent  on 
the  temperature 

.1.0024E-2 

Ohms/degree 

Kelvin*cmVcell 

Resistive  Constant 
coefficient  (r2) 

The  constant  portion  of 
the  resistance 

1.9202 

Ohms*cmVcell 

Exponential  Squared 
pressure  dependent 
coefficient  (mo) 

The  portion  of  the 
exponential  coefficient 
that  is  dependent  on  the 
square  of  the 
temperature 

1.314167823E-11 

millivolts/cell/Pascals'^ 

Exponential  linear 
pressure  dependent 
coefficient  (mO 

The  portion  of  the 
exponential  coefficient 
that  is  linearly 
dependent  on  the 
temperature 

-7.331599169E-6 

millivolts/cell/Pascals 

Exponential  constant 
coefficient  (m2) 

The  constant  portion  of 

the’exponential 

coefficient 

1.046788384 

millivolts/cell 
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Cell  Area  (Aceu) 

The  active  area  of  the 
membrane  electrode 
assembly 

292 

cm^ 

Number  of  cells  (Nceii) 

The  number  of  cells 
stacked  in  series 

60 

No  Units 

Heat  Creation  (Uo) 

Higher  Heating  value  of 
the  H2-O2  chemical 
reaction,  1.23  if  water 
vapor  is  the  end 
product,  1.482  if  liquid 
water 

1.482 

Joule-second/mole 

Newton  Steps 

Maximum  number  of 
iterations  to  attempt  to 
find  a  voltage-current 
operating  point 

25  1 

No  Units 

The  following  variables  are  computed  within  the  model  and  are  user  accessible. 


Variable  Name 

Description _ 

Units 

Current  Supplied 

Amount  of  current  sourced  from 
fuel  cell 

Amps 

Voltage 

Fuel  cell  electrochemical 
potential  of  positive  terminal 
with  respect  to  negative  terminal 

Volts 

Air  Consumed 

Amount  of  “air”  consumed  by 
chemical  reaction 

Kg/second 

H2  Consumed 

Amount  of  H2  consumed  by 
chemical  reaction 

Moles/second 

Electrical  Power  out 

Power  supplied  by  fuel  cell 
electrical  terminals 

Watts 

Heat  Power  out 

Power  expelled  by  thermal 
terminal 

Watts 

Electrical  Conductance 

Nonlinear  electrical  conductance 
of  fuel  cell 

Mhos 

Assumptions  in  Model  Derivation 

The  thermal  transfer  by  material  streams  is  negligible.  Water  management  effects  on  performance  are 
negligible.  This  is  not  always  the  case.  Air  is  consumed.  This  will  be  replaced  by  consumption  of  oxygen 
and  performance  dependence  on  partial  pressure  of  oxygen,  when  VTB  allows  mixed  material  stream 
connections.  This  will  also  allow  for  thermal  balance  of  material  streams. 

Mathematical  Description  of  Model 

For  each  cell  of  a  fuel  cell  stack, 

Eo  is  the  standard  potential  of  reaction  for  H2-O2  reaction,  b  and  r  are  functions  of  temperature.  And,  m 
and  n  are  functions  of  pressure.  Assuming  that  each  cell  has  equivalent  potential,  the  current  for  the 
complete  stack  can  be  found  using  iterative  methods.  Thus  the  overall  equation  for  the  electrical  terminals 

is: 
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Where  I  is  the  stack  current,  V  is  the  stack  voltage,  t  is  the  stack  temperature,  and  p  is  the  fliel/air  pressure. 
All  other  variables  are  defined  as  model  parameters  above. 

The  mass  flow  is  simply  proportional  to  current  and  number  of  cells. 

M  =  K,,rK-i 

The  proportionality  constant  K  is  in  units  of  Kg/Coulomb  for  air,  and  Moles/Coulomb  for  hydrogen,  and  is 
fixed  internally. 

The  total  power  contribution  to  the  system  from  the  reaction  is: 

\n-FJ 

Where  Uq  is  the  higher  heating  value  of  the  reaction.  The  difference  between  this  and  the  electrical  power 
release  is  the  thermal  power  that  must  be  expelled.  Thus,  the  thermal  power  expelled  from  the  fuel  cell 
stack  is: 

Q  =  Ncens-Uo-I-V-I 

Where  I  is  the  stack  current,  V  the  stack  voltage,  and  the  other  variables  are  defined  above  in  the  model 
parameters. 

Example  of  Model  Use 

The  following  Figures  show  a  test  schematic  that  can  be  used  to  evaluate  parameter  settings,  and  sample 
results. 


a  'a  B 


Model  Validation 

The  model  was  compared  with  experimental  data.  An  exact  match  is  difficult  to  achieve  and  will  change 
with  water  management,  construction  materials,  and  construction  methods.  Therefore,  the  model  only 
approximates  actual  performance  characteristics  of  the  desired  fuel  cell. 


Metal  Hydride  Bed  VTB  Model 


Author:  Zhenhua  Jiang 

Date:  03/18/2002 

Model  name:  HydrideBed 

DLL  name:  HydrideBed  _020318.vtm 

Version  number:  1.0 


Pictorial  Representation  of  Model 


Fig.l  Pictorial  representation  of  constant  flow  source  model 

Brief  Description  of  Model 

This  model  represents  the  dynamics  of  the  discharge  of  hydrogen  from  a 
cylindrical  metal  hydride  bed.  This  model  takes  into  account  the  heat  exchange  of  metal 
hydride  bed  with  the  external  heating/cooling  system  (Not  included  in  the  current 
version). 

Model  Validity  Range  and  Limitations 

This  model  is  valid  to  be  used  in  a  system  level  simulation.  It  can  be  connected  in 
series  with  a  gas  pressure  regulator  downstream  to  provide  hydrogen  at  a  relatively 
constant  pressure,  for  example,  to  the  fuel  cell.  This  model  can  be  used  for  actual  process 
design  and  hypothetical  process  feasibility,  respectively.  It  can  also  be  used  to  judge 
complexity  of  a  particular  hydrogen  discharge  process  by  comparing  with  experimental 
data. 

List  of  Terminals  with  Connectivity  Information 


Terminal  Designation 

Description 

H2 

Hydrogen  mass  flow  terminal 

List  of  Model  Parameters 


Parameter  Name 

Description 

Default  Value 

Units 

Initial  Pressure 

Initial  Pressure 

1500000.0 

Pa 

Initial  Temperature 

Initial  Temperature 

298.0 

K 

Initial  H2  concentration 

Initial  H2  concentration 

1.059 

Mol/mol 

Radius 

Radius  of  hydride  bed 

0.1 

m 

Length 

Length  of  hydride  bed 

0.25 

m 

CP  . . . . . . 

Porosity 

Hydride  bed  porosity 

0.4 

- 

Heat  capacity  of  metal 
hydride 

Heat  capacity  of  metal 
hydride 

10000.0 

J/(mol.K) 

List  of  Output  Variables 


Parameter  Name 

Description 

Units 

Pressure 

Pressure  across  the  source 

Pascal 

Molar  flow 

Molar  flow  rate  out  of  the  terminal 

Mol/second 

Temperature 

Temperature  inside  the  hydride  bed 

K 

H2  Concentration 

H2  Concentration 

Mol/mol 

Assumptions  in  Model  Derivation 

Ignored  are  various  mass  and  heat  transfer  resistances  in  both  the  axial  and  radial 
directions  and  internal  heat  exchange  with  a  circulating  fluid.  In  this  model,  isothemal 
local  equilibrium  conditions  prevail  and  the  isotherm  relationship  is  represented  by  a 
composite  Langmuir  model. 


Mathematical  Description  of  Model 


The  set  of  equations  below  represent  a  simplified  dynamic  model  of  the  discharge  of 
hydrogen  from  a  metal  hydride  bed  [1-5]. 
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where  both  andq*  are  functions  of  P  and  T  which  are  given  by 
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and  the  initial  conditions  are  as  follows:  T  =  T^,  P  =  -Po >  9  -  ^ *o > t  -  0 .  Note  that 

the  explanation  of  the  variables  used  here  is  given  in  Appendix  1. 

Model  Validation 

This  is  a  simplified  model  that  comes  from  a  set  of  complex  equations  describing 
the  dynamics  of  the  hydrogen  discharge  from  hydride  bed.  The  model  will  be  validated 
with  the  experimental  data. 

Model  Verification 

The  test  bench  consists  of  a  constant  flow  source  and  hydride  bed.  The  flow  r^e 
of  the  source  is  1.0  mol/second.  All  the  parameters  are  selected  as  the  default  values.  The 
schematic  diagram  of  the  simulation  system  is  shown  in  Fig. 2. 
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Fig  3.  Simulation  results  of  the  test  system 


Example  of  Model  Use 

See  Model  Verification  Section. 
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Appendix  1.  Variables  used  in  the  equations 

Hydride  bed  cross  sectional  area, 

Henry’s  law  constant  defined  in  Eq.Q,  Pa"* 

Temperature  independent  coefficient  of  the  Henry’s  law  constant,  Pa 
Heat  capacity  of  H2,  J  l{Kg  ■  K) 

Heat  capacity  of  the  metal  hydride,  J  l{Kg  ■  K) 

Diameter  of  metal  hydride  particles,  m 

Armenians  activation  energy  of  the  diffusion  coefficient,  J !  mol 
Mass  transfer  coefficient,  s  ' 

Hydride  bed  length,  m 

Molecular  weight  of  H2,  kg  /  mol 

Average  molecular  weight  of  metal  hydride,  kg  /  mol 

H2  molar  flow,  mol  /  s 

Pressure  inside  the  hydride  bed,  Pascal 

Pressure  inside  the  hydride  bed  at  the  initial  state,  Pascal 

Equilibrium  pressure  of  H2  in  metal  hydride  at  the  envelop,  Pascal 

Concentration  of  hydrogen  on  the  metal  hydride  bed,  mol  H2/mol  metal 

hydride 

Equilibrium  concentration  of  hydrogen  on  the  metal  hydride  bed,  mol 
H2/mol  metal  hydride 

Maximum  equilibrium  concentration  of  hydrogen  on  the  metal  hydride 
bed,  mol  H2/mol  metal  hydride 

Equilibrium  concentration  of  hydrogen  on  the  metal  hydride  bed  at  the 
initial  state,  mol  H2/mol  metal  hydride 

Saturation  concentration  of  hydrogen  on  the  metal  hydride  bed  at  the  low 

concentraion  phase,  mol  H2/mol  metal  hydride 
Ideal  gas  constant  (8.314  J / mol  •  K) 

Temperature  inside  the  hydride  bed,  K 

Temperature  inside  the  hydride  bed  at  the  initial  state,  K 

Isosteric  heat  of  adsorption  of  the  hydrogen  in  the  metal  hydride,  J/mol 

Hydride  bed  porosity,  dimensionless 

Density  of  the  metal  hydride,  kg  Im^ 


I  VanePumpli  VTB  Model 

Author;  Wei  Jiang 

Date;  06/10/02 

Model  name;  VanePumpl 

DLL  name;  VanePumpi_020610.dll 

Version  number; 

Report  errors  or  changes  to;  jiangw@ener.sc.edu 


Pictorial  Representation 


;  i 


Description 

This  model  represents  a  rotary  vane  pump. 


Validity  Range  and  Limitations 

1.  The  vane  pump  model  is  only  valid  for  small  (Mass  flow  Rate  is  around  6.0E-05  kg/s)  one-way  vane  pump. 

2.  The  equations  (1)(2)  for  pressure  and  pump  efficiency  are  derived  from  experimental  data  for  a  specific 
pump  (G12/04-E,  see  referencel).  So  this  model  is  only  valid  for  the  Pumps  similar  to  G12/04-E  type. 

3.  For  better  performance,  Rotational  speed  of  shaft  should  be  (2806  730)  rad/s 


Terminals  and  Connectivity  Information 


1  ms  mouei  nas  inc  luiiuw 

Variable  Name 

Description 

Type  of  variable 

Units 

Mass  Flow  Rate(  M  ) 

Mass  Flow  Rate  of 
fluid 

Through  Variable 

kgis 

Outlet  Pressure  ( p ) 

Outlet  Pressure  of  the 
pump 

Across  Variable 

Pa 

Torque(T) 

Input  torque 

Through  Variable 

N  ‘tn 

Rotational  speed  {CO) 

The  rotational  speed  of 
shaft 

Across  Variable 

radians!  s 

Parameters  and  Output  Variables 

This  is  a  complete  list  of  user  accessible  parameters.  SI  units  are  used  throughout. 


Parameter  Name 

Description 

Default  Value 

Units 

Density(  p  )) 

Density  of  the  gas 

1.205 

kRind 

Gas  constant  (  ) 

Gas  specific  heat 

1.005 

kJ/kgK 

Gas  constant  {k) 

Gas  specific  heat  ratio 

1.4 

No  Units 

Pressure 

atmosphere  (Pq) 

Pressure  atmosphere 

101000 

pa 

Temperature 
atmosphere  (Tq) 

Temperature 

atmosphere 

25 

"C 

The  following  variables  are  computed  within  the  model  and  are  user  accessible. 


Name 

Description  — - - 

Units 

Tnrniift 

The  torque  applied  to  the  Pump  _ _ _ 

N*m 

Rotational  Soeed 

Rotational  Speed  of  Pump 

rad/s 

Mass  Flow  Rate 

Mass  Flow  Rate  of  fluid 

kg/s 

Pressure  difference 

Pressure  difference  of  Pump 

pu 

Efficiency 

The  Efficiency  of  the  Pump 

No  Units 

Assumptions  in  Model  Derivation 

1.  The  process  in  the  vane  pump  is  adiabatic. 

2.  Air  is  considered  ideal  gas. 

3.  Pressure  difference  is  the  function  of  mass  flow  rate  and  rotational  speed. 

4.  Efficiency  is  the  function  of  mass  flow  rate. 


Mathematical  Description 

For  a  vane  pump  the  following  equations  encapsulate  the  useful  pump  characteristics: 

Ap  =  a-M^  +b-M  +  c- CO 
r}  =  d-M^  +e-M  +  f 


From  the  following  data  [1]  we  can  obtain  constant  a,b,C  and  d,e,f  when  CO  is  471.23  rad/s. 


t^rom  me  loiiowiug  udia  [j. 

Prp<?siire  ( rnhar^ 

0.0 

68.9 

103.4 

Efficiency  (NU) 

0.02% 

9.5% 

14 .07  % 

Mass  flow  (kg/s) 

10.27  E-05 

6.032  E-05 

3.95  E-05 

We  can  get 

a  =  -4.5782£+ll 
f>  =  -9.8635£  +  07 
c  =  31.746 


d  =  -7.5025£+06 
e  =  -1.16£;+03 
/  =  1.983£-01 


From  energy  conversion  law  and  ideal  gas  performance,  we  can  get  following  equations 


TCO  = 


Mc-(T,-T,) 


V 


T  =T 

■‘■2  -^0 


kPo  j 


^P  =  P2-Po 


Pi 

Po 

M 


-Pressure  difference 
—Outlet  pressure 
— Pressure  atmosphere 
— Mass  flow  rate 


(1) 

(2) 


(3) 


- Rotational  speed 

Yj  _ Vane  pump  efficiency 

j  . Outlet  temperature 

rr  _ Temperature  atmosphere 

^  0 

J  - Torque 

The  following  graphic  is  to  compare  experimental  data  with  results  from  equations  (1)  (2). 


Vane  Pump(G12/04-E) 


Example  of  Model  Use 

This  example  is  to  test  the  characteristics  of  air  vane  pump. 


Model  Validation  „ 

In  order  to  verify  this  Pump  model,  we  compare  VTB  output  results  with  Equations  (1)  (2)  just  as  follows  w  en 
is  471  rad/s. _ _ _ 


Mass  flow  rate  (kg/s) 


Output  pressure(pa) 


Efficiency 


VTB 


Equa. 


VTB 


4.6018E-05 


110444 


110450 


0.129 


0.1290 


5.4601E-05 

109202 


109210 


0.113 


0.1126 


6.705  lE-05 


107281 


107290 


0.087 


0.0868 


8.6642E-05 


103970 

103980 


0.041 


0.0415 
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On-Board  Fuel  Cell  Power  Plants 
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Tvnicallv  Pt  is  alloved  with  metals  such  as  Ru.  Sn.  or  Mo  to  provide  a  more  CO-tolerant.  high-perfomanCe  proton  exchange 
m^mbranc^fucl  celHPEMFC)  anode.  In  this  work,  a  layer  of  carbon-supported  Ru  is  placed  between  the  Pt  catalyst  and  the  anode 
flosv  irto  form  a  fiL.  When  oxygen  is  added  to  the  fuel  stream,  it  was  predicted  that  the  slow  H,  kmet.cs  of  Ru  m  th.s  fiber 
lould  become  an  advantage  compared  to  Pt  and  Pt:Ru  alloy  anodes,  allowing  a  greater  of  O;  to  oxtd,^ 

to  COj.  With  an  anode  feed  ofH,,  2%  O,.  and  up  to  100  ppm  CO.  the  Pt  +  Ru  filter 

Pt'Rii  allov  The  oxvceii  in  the  anode  feed  stream  was  found  to  form  a  hydroxyl  species  within  the  filter.  Tlie  reaction  of  these 
hvdroxy  groujs  tXadsorbed  CO  was  the  primary  means  of  CO  oxidation  within  the  filter.  Because  of  the  resultmg  proton 
foiation  the  Ru  filter  must  be  placed  in  front  of  and  adjacent  to  the  Pt  anode  and  must  contain  Nafion  in  order  to  prcide 
ionic  pathways  for  proton  conduction,  and  hence  achieve  the  maximum  benefit  of  the  filter. 
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Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  gaining 
popularity  due  to  their  benefits  such  as  environmental  friendliness 
and  increased  fuel  efficiency.  Because  of  the  difficulties  inherent  to 
storing  hydrogen,  liquid  fuels  such  as  propane,  natural  gas,  and 
gasoline  are  used  to  produce  reformate  gas.  Dry  reformate  is  typi¬ 
cally  composed  of  35-45%  hydrogen,  15-25%  carbon  dioxide,  50- 
10,000  ppm  carbon  monoxide,  and  a  balance  of  nitrogen.  It  has  been 
shown  extensively  that  CO  poisons  the  platinum  catalyst  used  in 
PEMFC  systems.*'^  Carbon  monoxide  chemically  adsorbs  onto 
available  Pt  catalyst  sites  as  shown  in  Eq.  1 

C0+Pt->C03d3  [1] 

At  concentrations  as  low  as  10  ppm,  CO  lowers  power  output  of 
the  PEMFC-containing  Pt  electrodes  by  50%."^'  The  addition  of  a 
preferential  oxidation  (PROX)  unit  to  the  fuel  processing  system  can 
reduce  CO  concentrations  in  the  reformate  gas  stream  to  approxi¬ 
mately  50  ppm.  Attempts  to  find  catalysts  both  tolerant  to  CO  and 
equivalent  in  performance  to  Pt  have  led  to  the  alloying  of  Pt  with 
Ru,  Mo,  W,  Co,  Ir,  Ni,  and  Sn.^‘‘”  Used  by  themselves  as  anode 
catalysts,  these  metals  do  not  provide  the  high  rate  of  hydrogen 
oxidation  necessary  to  achieve  the  current  densities  that  make  PEM¬ 
FCs  competitive  in  the  marketplace.’'*'*'^  The  most  commonly  used 
alloy  is  Pt:Ru.  The  Pt:Ru  alloys  combine  the  high  catalyst  activity  of 
Pt  with  the  increased  CO  tolerance  of  Ru."*'^ 

The  oxidation  of  CO^ds  the  platinum  catalyst  surface  in  the 
anode  shown  in  Ea  2  has  been  found  to  follow  Langmuir- 
Hinshelwood  kinetics 

CO.a,  +  CO2  +  +  2M  +  e"  [2] 

where  M  represents  Pt  or  Ru.  The  reactions  by  which  OHgds  ts 
formed  on  Pt  and  Ru  are  shown  in  Eq. 

M  +  HjO  -►  OH.d5  -b  H*  +  e"  [3] 

The  formation  of  OHgjst  shown  in  Eq.  3,  is  the  rate-determining 
step  of  this  reaction  and  occurs  on  platinum  at  potentials  of  0.7 
V/RHE  and  above. Ruthenium  has  the  ability  to  form 


*  Electrochemical  Society  Student  Member 
♦*  Electrochemical  Society  Fellow, 

**♦  Electrochemical  Society  Active  Member. 

*  E-maib  weidncr@engr.sc.edu 


from  water  at  significantly  lower  potentials  than  Pt,  0.35  V  for  50 
atom  %  Ru,  and  0,2  V  for  90  atom  %  Ru."*'^’'^  This  allows  the 
catalytic  desorption  of  CO  as  CO2  to  commence  at  lower  potentials. 
There  is  a  linear  relationship  between  the  onset  of  CO  oxidation  and 
Ru  composition  (/.e.,  the  shift  in  potential  is  linear  with  respect  to 
atomic  fraction  of  Ru  in  the  alloy).  However,  the  benefit  of  alloying 
Pt  with  Ru  has  only  been  shown  to  provide  near  equivalent  perfor¬ 
mance  to  pure  H,  on  Pt  for  CO  concentrations  up  to  100  ppm^in  the 
feed  stream  for  low-temperature  fuel  cell  PEMFC  operation.  ' 

The  injection  of  oxygen  into  the  fuel  stream  has  also  been  shown 
to  increase  catalyst  tolerance  to  CO.'^*'^  This  “air  bleeding  pro¬ 
vides  a  greater  concentration  of  active  oxygen  on  the  catalyst  that 
will  then  react  with  CO  to  form  CO2 .  The  following  reactions  are 
assumed  to  occur  at  the  Pt  catalyst  surface  (in  addition  to  Eq.  1) 

O2  +  2Pt  — >  20ads  W 

C03ds+  Ogds-^  CO2  +  2Pt  [5] 

2H3d..+  03ds->H20+3Pt  [6] 

The  goal  of  the  addition  of  an  air  bleed  to  the  anode  feed  is  to 
promote  the  reaction  shown  in  Eq.  5,  the  oxidation  of  COg^s »  lo 
CO2.  Note  the  difference  between  Eq.  2  and  5.  The  reaction  of 
COgas  with  OHgds  on  Pt  requires  an  overpotential,  and  hence  ionic 
contact  to  the  cathode  region  (in  the  form  of  the  Nafion  electrolyte), 
whereas  the  reaction  shown  in  Eq.  5  does  not.  One  of  the  problems 
of  adding  O2  to  the  gas  inlet  stream  is  that  the  reaction  sho\yn  in  Eq, 
6  occurs  at  a  much  higher  rate  than  the  reaction  shown  in  Eq,  5. 
Only  about  one  of  400  O2  molecules  oxidizes  COgds  to  CO2 .  This, 
and  the  fact  that  the  mixture  becomes  combustible  for  concentra¬ 
tions  of  oxygen  in  hydrogen  above  4  vol  %,  limits  the  amount  of 
oxygen  in  the  feed  stream  and  limits  the  effectiveness  of  the  air 
bleeding  technique. 

The  first  objective  of  this  study  was  to  determine  if  depositing  a 
layer  of  carbon-supported  Ru  on  top  of  a  typical  carbon-supported 
Pt  anode  (see  Fig.  1)  increases  the  effectiveness  of  the  air-bleeding 
technique  in  preventing  CO  poisoning.  By  first  coming  in  contact 
with  Ru  instead  of  Pt,  it  was  speculated  that  a  larger  percentage  of 
the  O2  in  the  H2/CO/O2  feed  will  react  wth  COgds  than  in  either  a 
Pt  or  a  Pt:Ru  alloy  electrode.  The  selectivity  of  Eq.  5  and  6  would 
then  shift  toward  the  oxidation  of  COgds .  Studies  have  sho^vn  that 
the  rate  of  hydrogen  oxidation  on  pure  Ru  at  62°C  (roughly  the 
operating  temperature  used  in  this  experiment)  was  two  orders  of 
magnitude  lower  than  that  of  Pt  measured  under  similar 
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Figure  1.  Diagram  of  the  Ru-filtcred  anode.  The  feed  gases  through  the 

GDL  and  come  in  contact  with  the  Ru  filter  before  reaching  the  Pt  anode 
The  Rii  filter  acts  as  a  chemical  barrier  on  which  Oj  oxidizes  the  CO  present 
in  the  feed  stream  to  COj .  preventing  a  loss  in  the  Pt  anode  performance  due 
to  CO  poisoning. 


conditions.'^  CO  has  been  shown  to  have  similar  adsorption 
strengths  onto  Ru  as  Pt.'^'”’'’  Thus,  the  slow  H,  kinetics  of  pure  Ru 
would  become  an  advantage,  allowing  a  greater  percentage  of  Oj  to 
oxidize  CO,*  to  COj  rather  reacting  with  adsorbed  hydrogen  to 
form  water.  Pt  catalyst  is  then  placed  after  Ru  to  oxidize  and 

maintain  the  high  PEMFC  performance.  .  ■ 

The  second  objective  was  to  characterize  the  mechanism  of  CO 
oxidation  occurring  within  the  Ru  filter  by  varying  the  Nafion  con¬ 
tent  (10  or  35  wt  %)  within  it  and  its  placement  along  the  anode.  The 
Ru  filter  can  be  placed  between  the  gas  diffusion  layer  (GDL)  and  P 
anode  as  shown  in  Fig.  1,  or  between  the  flow  '^annels  and  the 
GDL.  By  varying  the  placement  of  the  Ru  filter,  the  role  of  Nafion 
electrolyte  in  CO  oxidation  could  be  determined. 


Experimental 

Development  of  Pt  and  Pt.Ru  catalyst  The  method  de_ 

scribed  in  Patent  no.  5,211,984  provides  an  outline  for  NCI  and 
membrane  electrode  assembly  (MEA)  manufacture  'p 
The  Pt  and  Pt:Ru  (1:1  atomic  ratio)  catalysts  were  bonded  to  the 
Nafion  112  proton  exchange  membrane  (PEM)  in  several  steps. 
These  steps  included  the  formulation  of  the  catalyst  inks,  application 
of  the  inks  to  decals,  and  transference  of  the  dried  catalyst  •nk  frofn 
the  decal  to  the  membrane.  Inks  were  prepared 
adding  the  E-TEK  catalyst  to  a  solution  of  5  wt  %  Nafion  (DuPont). 

The  components  of  the  ink  were  added  to  an  appropriate  size 
bottle  (evacuated  with  helium  to  avoid  sparking)  and  then  stirred  for 
a  minimum  of  8  h  to  ensure  uniformity  of  the  ink.  Three-ply  Teflon 
decals  were  weighed  prior  to  applying  the  application  of  the 
ink.  The  ink  was  drawn  across  the  surface  of  the  Teflon  decals 
using  a  Meyer  rod.  The  loadings  of  0.40  mg  Pt/cm^  (anode), 
0  61  mg  Pt:Ru/cm’,  and  0.50  mg  Pt/cm’ (cathode)  were  achieved. 
The  coated  decals  were  then  dried  in  an  oven  at  105“C  and  ambient 
pressure  for  10  min  to  remove  any  remaining  alcohols  from  the  ink 

To  form  an  MEA  with  an  active  area  of  50  cm  decals  coated 
with  an  anode  and  cathode  were  placed  on  either  side  of  a  Nafion 
112  membrane  (proton  form).  The  assembly  was  then  hot-pressed 


Table  I.  Types  of  MEUs  tested. 

MEU  name _  Dcscnpt\on _ 

Pt  (baseline)  0.4  mg  Pt/cm^  anode 

0.5  mg  Pt/cm^  cathode 

Uncatalyzed  Zoltek  GDLs  on  anode  and  cathode 

pt:Ru  0.61  mg  Pt:Ru/cm-  anode 

0.5  mg  Pt/cm^  cathode 

Uncatalyzed  Zoltek  GDLs  on  anode  and  cathode 
The  atomic  ratio  of  Pt  to  Ru 
in  the  Pt:Ru  alloy  is  1:1,  ^ 

resulting  in  an  anode  Pi  loading  of  0.4  mg/cm  . 

Pi  +  Ru  (Ru  filter)  0.4  mg  Pt/cm-  anode 
0.5  mg  Pt/cm^  cathode 
0.21  mg  Ru/cm^  coaled  Zoltek  GDL 
on  the  anode 

Uncatalyzed  Zoltek  GDL  on  the  cathode  side 


for  2  min  at  205°C.  The  assembly  was  cooled  to  room  temperature 
before  the  decals  were  carefully  pealed  from  the  assembly,  leaving 
an  MEA. 

Development  of  the  Ru  filter.— Ru  filters  on  both  the  catalyst  and 
flow  channel  sides  of  the  GDL  were  studied.  Ru  inks  containing  35 
and  10%  Nafion  solids  by  weight,  respectively,  were  prepared  in  a 
similar  method  to  Pt  and  Pt:Ru  catalyst  inks.  The  ruthenium  ink  was 
applied  to  the  microlayer  side  of  an  uncatalyzed  Zoltek  GDL  with  a 
cross-sectional  area  of  50  cm^  A  Meyer  rod  was  used  to  achieve  a 
target  loading  of  0.21  mg  Ru/cm^  After  the  Ri^C  catalyst  was  ap^ 
plied  the  GDLs  were  dried  at  105°C  and  ambient  pressure  for  10 
min.  Ru  filters  placed  on  the  flow  channel  side  of  the  GDL  were 
prepared  in  a  similar  manner  and  with  identical  loadings. 

Cell  assembly  and  testing, — The  three  types  of  50  cm 
units  (MEUs)  prepared  and  tested  are  shown  m 


Table  I.  -i  ,  „  .4 

The  MEUs  were  placed  in  a  50  cm^  cell.  The  cell  was  assembled 
and  incubated  for  4-8  h  at  ambient  pressure,  cell  ternperature  ot 
70°C  an  anode  feed  of  hydrogen,  a  cathode  feed  of  a>r.  h™  ® 
stoichiometric  ratio  (actual  flow/stoichiometric  flow  required  for  a 
1.0  A/cm^  current)  of  1.5  at  the  anode  and  2.0  at  the  cathode. 
Single-cell  performance  curves  were  obtained  under  the  conditions 

set  in  Table  11.  ,  •  i  • 

An  alternate  test  station  was  used  to  conduct  trials  using  LU 

concentrations  greater  than  50  ppm  (Hj,  H2  +  50  ppm  CO, 

+  100  ppm  CO,  and  Hj  +  200  ppm  CO).  The  data  attained  on 
this  station  were  consistent  and  reproducible,  however,  prob  ems 
arising  from  excess  water  entering  the  test  cell  on  these  tnals  limit 
comparisons  of  this  data  to  a  qualitative  nature. 


Results  and  Discussion 


Performance  of  the  Ru  filter.-Figutes  2  and  3  show  the  perfor¬ 
mance  of  a  Pt;Ru  and  Pt  +  Ru  filter  anode  under  five 
anode  feeds:  (/)  hydrogen,  (//)  reformate;  (i/i)  reformate  +  2/« 
O,  (/u)  reformate  +  1%  O^.  and  (u)  reformate  +  0.5/.,  O,.  In 
the' case  of  Pt:Ru,  as  the  amount  of  oxygen  in  the  feed  is  reduced, 
the  performance  decreases  in  a  continuous  fashion.  For  a  Pt  +  Ru 
anode,  there  is  very  little  difference  in  performance  between 
reformate  -F  2%  O,  and  reformate  -f  1%  Oj,  but  as  the  concen¬ 
tration  of  the  oxygen  in  the  reformate  is  reduced  to  0.5%,  there  is  a 
sudden  drop  in  performance  such  that  the  performance  cuwe  re¬ 
sembles  that  of  pure  reformate.  For  conditions  of  1-2%  Oj  (by  vol¬ 
ume  of  H,)  addition  to  the  reformate  feed  stream,  the  Pt  +  Ru  filter 
performed  better  than  the  Pf.Ru  alloy  anode. 

Figure  4  compares  the  cell  performance  of  Pt,  Pt:Ru,  and  Pt 
+  Ru  filter  anodes  under  conditions  of  reformate  +  1.0%  O2 
bleed.  Under  these  conditions,  the  Ru  filter  outperforms  a  Pt:Ru 
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Table  II.  Fuel-cell  test  conditions. 

Pressure  _ 

Cell  temperature 

Stoichiometric  ratio  (at  I  A/cm  ) 


1  atm  (anode  and  cathode) 
70®C 

1.5  Anode  (hydrogen) 

2.0  Cathode  (air) 


Fcedstreams 


Anode:  hydrogen,  reformate, 
reformate  +  air  bleed, 

H2  +  CO  +  air  bleed 
Cathode:  air 


Dry  reformate  composition 


40%  Hz 
20%  CO2 
50  ppm  CO 
balance  Nz 


Humidification 


Complete  humidification 
of  anode  and 

cathode  gas  streams  for  all  trials. 


Figure  3.  Single-cell  performance  comparison  of  the  Pi  +  Ru  filter  for  van- 
ous  anode  feeds.  P  =  I  atm,  T  =  70°  C. 


Air  bleed 


CO  amounts 


0.5,  1.0,  2.0%  O2  (in  the  form 
of  an  air  bleed  relative 
to  the  volumetric  flow  of 
hydrogen  in  slm). 

50  ppm  in  reformate 
50,  100,  200  ppm  in  Hz 


alloy  anode  of  an  identical  loading  by  the  greatest  margin.  At  0.6  V, 
the  current  density  produced  from  the  cell  containing  the  Pt  +  Ru 
filter  is  almost  double  that  of  an  identical  cell  using  the  Pt:Ru  alloy. 
The  Pt;Ru  shows  performance  improvement  over  the  plain  Pt  anode 

only  at  voltages  below  0.5  V.  n 

inversely  in  Fig.  5.  the  cell  performance  of  the  Pt:Ru  alloy  for 
reformate  +  0.5%  O,  shows  the  greatest  tolerance  to  CO  of  the 
three  anodes.  At  0.6  V,  the  Pt:Ru  cell  provides  50%  more  current 
than  the  Pt  +  Ru  filter.  Both  Ru-containing  anodes  do  perform  bet¬ 
ter  than  the  pure  Pt  anode,  but  there  is  little  benefit  gained  from 
using  the  Pt  +  Ru  filter  instead  of  a  pure  Pt  anode.  Similarly,  m  Fig. 
6,  Whh  no  air  bleed  in  the  reformate  stream,  the  Ru  filter  shows 
performance  resembling  that  of  a  pure  Pt  anode,  while 
alloy  anode  exhibits  some  CO  tolerance.  In  all  cases  these  results 
were  reproducible  over  several  trials  using  different  MEAs  of  iden- 

'"'^Thusfwfthout  sufficient  oxygen  in  the  anode  feed  to  oxidize  CO 
to  CO.,  the  CO  penetrates  the  Ru  filter  and  poisons  the  Pt  region  of 
the  anode.  A  representation  of  this  is  shown  in  Fig.  7  where  the 


Figure  2.  Single-cell  performance  comparison  of  PtiRu  for  various  anode 
feeds.  P  ~  I  atm,  T  =  70°  C. 


Fieure  4.  Single-cell  performance  comparison  of  Pt,  PtiRu,  and  Pt  +  Ru 
filter  for  an  anode  feed  of  reformate  +  1%  O,  bleed.  P  -  atm,  T 


Figure  5.  Single-cell  performance  comparison  of  Pt,  PtiRu  and  Pt  +  Ru 
filter  for  an  anode  feed  of  reformate  +  0.5%  O,  bleed.  /■  -  1  a  m. 


shading  indicates  the  relative  concentration  of  CO  in  the  ^ 

CO  is  oxidized,  the  concentration  of  CO  decreases  as  it  p 
through  the  filter  region.  Figure  7a  is  indicative  of  adding  l-2  /»  ^ 
to  the  reformate  stream.  All  CO  is  oxidized  from  the  feed  stream 
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Figure  6.  Single-cell  performance  comparison  of  Pt,  Pt:Ru.  and  Pt  +  Ru 
filter  for  an  anode  feed  of  reformate.  P  =  I  atm,  T  -  10 


before  it  reaches  the  Pt  portion  of  the  electrode  and  * 

cell  operates  with  minimal  performance  losses.  F'gure  '  h ‘mmh 
that  when  the  Ru  filter  fails  to  oxidize  CO  before  it  diffuses  through 
the  filter  region,  the  remaining  CO  then  encounters  the  pure  ^  e 
trode  where  it  poisons  the  Pt  through  the  reaction  shown  in  Eq.  1. 

A  decrease  of  the  Oj  concentration  in  the  anode  feed  (Fig.  3) 
or  an  increase  in  the  CO  concentration  in  the  anode  feed  stream 
results  in  the  eventual  failure  of  the  Pt  -1-  Ru  filter  anode.  Figure  8 
shows  the  performance  of  the  Ru  filter  anode  when  using  an  anode 
feed  of  hydrogen,  2%  Oj ,  and  various  levels  of  CO.  When  50  ppm 
•CO  is  present,  there  is  almost  no  reduction  m  performance 
(0.43-0.42  A/cm^  at  0.6  V).  This  is  inconsistent  with  the  same 
amount  of  CO  contained  in  a  reformate  stream  (Fig.  3)  and  can  be 
attributed  to  the  dilution  of  hydrogen  in  the  reformate  stream.  (Com¬ 
pletely  humidified  at  70°C,  reformate  -I-  2%  Oj  m  the  form  of  an 
air  bleed  contains  28%  by  volume  compared  to  a  70  vol  /.  Hj 
for  an  H,  +  2%  O^  feed.)  Because  hydrogen  must  diffuse  ‘bro^h 
the  Ru  filter  in  order  to  react  on  the  Pt  electrode,  additional  diffu- 
sional  resistance  accounts  for  the  remaining  performance  losse 
compared  to  the  perfomance  of  a  baseline  MEA  with  an  anode  feed 
of  pure  hydrogen.  As  the  CO  level  increases  to  100  ppm,  the  current 
density  decreases  from  0.42  to  0.38  A/cm  at  0  _6  Y’ 
the  limit  of  the  Ru  filter  may  have  been  reached  At  200  Ppm  CO 
the  drop  in  performance  is  much  larger  (0  38-0.^  Mem  at  0  6  V). 
In  this  case,  not  all  the  CO  is  oxidized  m  Ru 
resulting  in  the  poisoning  of  the  Pt  portion  of  the  electrode.  Th  s 
trend  is^  qualitatively  similar  to  the  reduction  of  the  air  bleed  as 
shown  in  Fig,  3. 


Fiuurc  7.  A  schematic  representation  of  the  species  present  m  the  various 
regions  of  the  Pt  +  Ru  filler  anode.  The  shading 

centration  of  the  CO  in  each  region  (the  ''ghter  the  shade,  the  lower  the 
concentration).  In  (a)  all  CO  is  oxidized  within 

catalyst  receives  CO-free  gas.  In  (b)  the  concentration  of  CO  m  the  inlet  gas 
is  too  great  to  be  completely  oxidized  within  the  filter  and  poisons  the  anode 
catalyst. 


Mechanism  of  the  Ru  fiUer.-Equttlxons  2  and  5  present  two  dif- 
irent  methods  for  CO  oxidation.  In  the  mechanism  shown  in  bq.  2, 
le  oxidation  of  CO  generates  H^,  whereas  the  CO  oxidation  de- 
'ribed  in  Eq.  5  produces  no  ions.  Thus,  by  removing  the  electrolyte 
lat  allows  H*  to  diffuse  away  from  the  reaction  site,  CO  oxidation 
irough  the  mechanism  in  Eq.  2  can  be  effectively  prevented  This 
nables  the  CO  oxidation  through  Eq.  5  to  be  isolated  and  analyzed. 

The  isolation  of  the  Ru  filter  was  achieved  by  applying  o  me 
lowfield  side  of  the  electrolyte,  thus  separating  it  from  the  Pt  e  ec- 
rode  by  the  GDL.  The  Toray  GDL  contained  no  Nation  electrolyte 
nd  was  sufficiently  thick  (175  ptm)  to  provide  adequate  separation 
rom  the  Pt  electrode.  To  further  characterize  these  two  mechanisms, 
tu  filters  containing  10  wt  %  Nafion  were  tested  on  both  the  mem- 

irane  and  flowfield  sides  of  the  GDL. 

Figure  9  shows  the  effect  of  varying  the  Nafion  content  within 
he  Ru  filter  and  the  effect  of  the  placement  of  the  Ru  filter  for  an 
mode  feed  of  reformate  +2%  O,.  Of  the  four  types  tested,  the  Ru 
liter  placed  on  the  membrane  side  of  the  GDL  and  containing  35  wt 

r  .  .  .  .  _ C _ DrtfVi  Pn  fi  tPT  fvne.'? 


Figure  9.  Single-cell  performance  comparison  of  the  Pt  +  Ru  fiber  for  an 
Slfeed  of  mformate  -P  2%  O,.  P  =  1  atm  F  = 
figurations  are  defined  as  follows:  filter  ^ 

N^on  and  is  placed  on  the  membrane  side 

contains  35  wt  %  Nafion  and  is  placed  on  the  flowfield  side  of  the  GDL. 
filter  no  3  (— C3— )  contains  10  wt  %  Nafion  and  is  placed  on  the  membrane 
^de  of  ie  GDL;  filter  no.  4  (-A-)  contains  10  wt  %  Nafion  and  is  placed 
on  the  flowfield  side  of  the  GDL. 


A866 

placed  on  the  flowfield  side  of  the  GDL 

over  the  pure  Pt  anode,  providing  0.14  and  0.16  ^cm  . 

A/cm^  XheRu  filter  on  the  membrane  side  of  the 

K  nta  0  wt  %  performed  sirnilar  to  the  35  v. 

Ster  for  0  6  V  and  above,  but  the  lack  of  electrolyte  present  most 
likelv  induced  H*  diffusion  limitations  at  a  lower  cuircnt  density. 

tLs  it  is  believed  that  the  following  mechanism  is  occurring  m 
the^ruthenium  filter  for  a  fully  humidified^  foel  stream  containing 


Journal  of  The  Electrochemical  Society,  149  (7)  A862-A867  (2002) 


Ru  +  Ht  2Ha<is 

[7] 

Ru  +  CO  ^  COajs 

[8] 

O2  +  2Ru  2  Oads 

[9] 

H,,,  +  0,„  -  OH,„  +  Ru 

[10] 

-F  OH,„  -b  CO2  +H^  +  2Ru  +  e- 

[11] 

+  e"  -h  OH,d,  •r-b  H2O  +  Ru 

[12] 

CO,,,  -F  0,*  -*  CO2  -F  2Ru 

[13] 

H,„  H^  -F  e-  +  Ru 

[14] 

Reactions  7-9  represent  the  adsorption  of  species  onto  the  Ru 
cataWst  Reaction  10  represents  an  intermediate  reaction  on  Ru  re- 
sultii  in  the  formation  of  Like  on  Pt,  oxygen  dissociatively 

peting  desorption  reactions  on  the  Ru  catalyst.  >3 's  we 

LuLnted  as  following  Langmuir-Hmshelwood  k-ne.  s  on  Pt  and 
n„  34.35  Of  the  desorption  reactions,  Eq.  11  and  13  are  oesirea  as 
Stey  resuU  in  die  oxidation  of  CO,,-  above-descnbed  mecha- 
Ssm  is  not  a  comprehensive  list  of  all  reactions 
anode  region.  It  focuses  on  those  reactions  that  lead  to  CO  oxi 
rn  Fm  example,  evidence  has  been  found  that  ts  formed  a 

“”nSi..u  Ririns  ox«..  red.c.ion  ..  ••  .s 

0.5  V.“  However,  breaks  back  down  to  ® 

compounds  at  or  even  before  reaching  the  catalyst  surface. 

The  mtionale  for  the  placement  of  the  filter  on  the  Aowfie^d  side 

of  the  GDL  is  that  if  the  oxidation  of  CO,,,  were  ‘o  Pjoce^  ^a  Eq^ 
11  a  three-phase  interface  (anode  gas,  Nafion,  and 
needed  in  the  catalyst  mixture  because  there  are  no  tons  or  electrons 
o  be  tr^sported  aLy  from  the  reaction  Z 

is  necessary  in  the  Ru  filter  to  act  as  a  binder,  Ae  Ru  filter  was 
merely  placed  on  the  flowfield  side  of  the  GDL.  There  7"°^^ 
electrolyte  in  the  GDL  to  allow  proton  transport  and  the  GDL 
sufficiently  thick  (~5  mil).  Without  the  Nafion  electrolyte  as  a  con¬ 
duit  and  because  of  the  long  distance  that  the  protons  must  trave  , 
the  formation  of  CO,  through  Eq.  11  does  ^ 

cant  extent  and  ‘“h^  ^gTl^e  B  is 

amount  ox.dued  m W 

Ftg°  s'and  8.  This  is  consistent  with  literature,  where  it  is  also  shown 
Sft  Ru  l  one  of  the  least  active  metals  for  this  method  of  CO 
SdaL  at^^the  oxygen  partial  pressures  used  in  these 

experiments.^*’ 


‘’^uSevidence  of  the  role  of  Nafion  content  in  ‘hf  ^mbrnne 
is  shown  for  the  cases  in  Fig.  9.  Significant  CO  oxtdationonly 
occurs  when  the  Ru  filter  is  placed  adjacent  to  the  Pt  anode.  TTus  is 
ScTthat  protons  are  formed  as  product  of  CO  oxidation  Jhe 
Nafion  electrolyte  is  the  necessary  conduit  for  the  protons  to  move 
away  frl  the  reaction  site  and  when  this  conduit  is  Pre^.  cell 

perfLance  increases  when  all  Sn 

The  increased  performance  under  conditions  of  increased  Nano 
tonf  stinwn  in  Fie  9  is  further  proof  that  CO,,,  oxidation  is 

a  small  ’extent  as  evidenced  by  the  slight  increase  in  performance 


when  the  filter  is  placed  on  the  flowfield  *7  f  ‘J  heTthe  SteHs 
rate  of  Reaction  11  is  far  greater  than  Reaction  13  when  the  filter 

"mVdroxTgrru^in  >1  T 

oxygen  through  Ructions  9  and  10  or  from  water  through  Reaction 
12  fn  fact,  both  are  occurring  within  an  Ru  filter  when  « 
to  the  Pt  anode.  The  mechanism  primarily  responsible  tor  the  CO 
tolerance  of  the  Pt:Ru  alloy  is  comprised  of  Reactions  11  and  12. 
ns5“  staw,  ,h..  i„  te  .b.c.=  or ..  .i,  b'f. 
small  performance  increase  relative  to  a  pure  Pt  anode.  Thus,  the 
formation  of  OH,,,  through  the  reverse  of  Reaction  12  ts  no^  rte 
primary  source  of  the  hydroxyl  species  needed  for  Eq.  1 1^  i  77 
that  the  Pt’Ru  alloy  performs  better  under  reformate  con^ions. 
Gasteieer"’'-’''*’’'’  performed  extensive  research  on  this  mechanism 
and  foLd  that  while  OH,,,  is  indeed  formed  at  ‘ower  potentials  on 
Ru  than  Pt,  the  hydroxyl  group  also  forms  a  stronger  bon  o  P 
Ru  than  a  1:1  atomic  ratio  Pt;Ru  alloy,  resulting  in 
Sion  potentials  on  pure  Ru  than  on  Pt:R«  ‘ 

^inpciilated  that  the  rate  of  formation  of  OHads 
equivalent  to  that  in  the  Pt:Ru  alloy.’o  However, 

6  that  without  the  presence  of  oxygen  in  the  anode  feed  that  only 
Action  of  the  CO  is  oxidized  in  the  Ru  filter  and  that  the  remaining 

^^Tto^thVjrimaS'meAod  by  which  CO  tolerance  is  achieved  by 
the  Ru  filter  is  by  the  fomiation  of  OH,„  from  oxygen  in  the  air 
ll  eSThe  explanation  of  near  complete  oxidation  of  CO  at  concen- 

trations  up  to  100  ppm  is  the  oVi  0  to  he 

CO  free  oerformance  is  achieved  with  the  addition  oH/o  O2  to  the 
anode  feed  while  the  filter  provides  almost  no  benefit  u^er  refo^ 
^ate  cSditions  is  evidence  that  Reactions  9-11  are  the  primary 
means  of  CO  oxidation  within  the  filter  region.  ^ 

The  key  is  that  O,,,  forms  OH,„  before  7* 

proton  to  form  water.  This  is  where  the  benefit  of  the  Ru  filter  1  es^ 
Se  oxidation  of  H,„  occurs  within  the  filter  according 
4  forming  H^.  Because  the  kinetics  of  hydrogen  adsorption  on  Ru 
is  several  orders  of  magnitude  slower  than  on  Pt,  ^ 

concentration  of  H,i,  in  the  filter  region  and  more  ^ 

fnr  m  to  react  with  OH,,, ,  resulting  in  the  formation  of  CO2 . 
The  bulk  of  the  hydrogen  oxidation  occurs  in  the  Pt  region.  How- 
Iver  omtons  must  be  present  within  the  Ru  filter  for  Reaction  10  to 
occur,  and  because  these  protons  will  also  sh^  Ru°fiiter 

the  formation  of  H2O,  there  is  a  limit  to  the  benefit  of  the  Ru  filter. 

Conclusions 

For  an  anode  feed  stream  consisting  of  reformate  (containing  50 
rni  -mrf  1  2%  oxygen  the  Pt  -h  Ru  filter  electrode  shows 

rJ»d’co «  ^ r,ss“  1 

amounts  of  Pt  and  Ru.  For  CO  concentrations  up  to  100  ppm  and  2 
vnl  %  O.  the  Pt  +  Ru  filter  anode  also  shows  supenor  P^fPof' 
mtce.  U  h  Sely  that  the  oxidation  of  the  CO  within  the  Ru  filler  is 
primarily  due  to  oxygen  reacting  to  form  OH,*,  whic  en 
Lchemically  reacts  with  CO,,,  to  form  CO2  and  protons. 

However,  with  insufficient  oxygen  (<1  vol  %)  or  too  much 
CO  (>100  ppm),  not  all  CO  is  oxidized  in  the  Xr 

CO  reaches  and  then  poisons  the  Pt  region  of  the  Pt  -1-  Ru  filter 
Sde  M  a  result,  cell  performance  under  those  conditions  is  worse 

*''Blcauw'ben”fits  of  the  Ru  filter  occur  at  high  levels  f' f  ^ 
nv  O  )  and  the  Pt’Ru  alloy  provides  CO  tolerance  even  with^  air 
b  e  d  it  is  sSgesS  that  the  anode  configuration  that  wo^  pro- 

.Idl  opdll  CO  ..did  con*,  of  .0  R.  61.«r  pte"  ■« 

front  of  and  adjacent  to  a  Pt:Ru  alloy. 
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Using  Sputter  Deposition  to  Increase  CO  Tolerance 
in  a  Proton-Exchange  Membrane  Fuel  Cell 

Andrew  T.  Hang,'  *  Ralph  E.  White, John  W.  Weidner,*’***  Wayne  Huang," 

Steven  Shi,"’***  Narender  Rana,*  Stephan  Grunow,'  Timothy  C.  Stoner, 
and  Alain  E.  Kaloyeros'' 

^Center  for  Electrochemical  Engineering.  Department  of  Chemical  Engineering.  University  of  South 
Carolina.  Columbia.  South  Carolina  2920S.  USA 

^Phtp  Power  Incorporated.  Latham.  New  York  12110.  USA  .  .  .  t 

^Ener^’  and  Environment  Applications  Center.  Institute  for  Materials,  and  School  ^  and 

Engineering.  State  University  of  New  York  at  Albany.  Albany.  New  york  12205.  USA _ _ 

Placinc  a  layer  of  Ru  atop  a  Pt  anode  increases  the  carbon  monoxide  tolerance  of  proton-exchange  membrane  fuel  ce'ls  when 

oLfoiroancc  of  (he  Pt  +  single-layer  sputtered  Ru  filter  was  double  that  of  the  Pt:Ru  alloy  (0.205  vs.  0.103  /Vcm  at  0.6  V).  The 
lltZZl  wal  also  sVnificantly  greLr  than  that  of  the  ink-based  Ru  filter  (0.149  A/cm=  at  0.6  V).  W.th.n  the  filter  region  of 
fhe^anode,  it  is  likely  that  the  decreased  hydrogen  kinetics  of  the  Ru  (compared  to  Pt)  allow  for  more  of  the  OH,j,  formed  fro 

oxvecn  in  the  fuel  stream  to  oxidize  adsorbed  CO  to  CO2.  ,  .  . 

©  2002  The  Electrochemical  Society.  [DOI:  10.1 149/1.1479727]  All  nghts  reserve  . 
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Proto n-e.x change  membrane-based  fuel  cells  (PEMFCs)  are  gain¬ 
ing  popularity  due  to  their  high  operating  efficiency  and  environ¬ 
mental  friendliness.  Because  of  the  difficulties  inherent  to  storing 
hydrogen,  hydrocarbon  fuels  such  as  propane,  natural  gas,  and  gaso¬ 
line  are  used  to  produce  reformate  gas.  Dry  reformate  is  typically 
composed  of  40-75%  hydrogen,  15-25%  carbon  dioxide,  10-10,000 
ppm  carbon  monoxide,  and  a  balance  of  nitrogen,  depending  on  the 
fuel  processing  system  used.*-^  It  has  been  shown  extensively  that 
CO  poisons  the  platinum  catalyst  used  in  PEMFC  systems.  ^  Car¬ 
bon  monoxide  chemically  adsorbs  onto  available  Pt  catalyst  sites  as 
shown  in 

CO  +  Pt  CO,ds 

As  little  as  10  ppm  CO  in  the  fuel  stream  can  lower  the  power 
output  of  the  PEMFC  by  50%.^’’  For  the  reforming  process  to  be 
effective  in  the  fuel  cell  system,  this  problem  must  be  solved.  At¬ 
tempts  to  find  catalysts  both  tolerant  to  CO  and  equivalent  in  per¬ 
formance  to  Pt  have  led  to  the  alloying  of  Pt  with  Ru,  Mo,  W,  Co, 
Os,  Ir,  Ni,  and  Sn.’'^‘  Used  by  themselves  as  catalysts,  these  metals 
do* not  provide  the  high  hydrogen  activity  necessary  to  achieve 
the  current  densities  that  make  PEMFCs  competitive  in  the 
marketplace.®-’^’*^  The  most  commonly  used  alloy  is  Pt:Ru.  The 
Pt:Ru  alloy  combines  the  high  catalyst  activit;jj^^of  pure  Pt  with  the 
increased  CO  tolerance  of  pure  Ru  catalyst. 

The  oxidation  of  CO^d^  from  the  Pt:Ru  catalyst  surface  in^the 
anode  shown  in  Eq.  2  follows  Langmuir-Hinshelwood  kinetics 

CO3J,  +  OH,j,  CO2  +  H*  +  2M  +  e-  [2] 

where  M  represents  Pt  or  Ru.  The  reactions  by  '^bich  OH,d5  is 
formed  on  Pt  and  Ru  are  shown  in  Eq.  3'  ’ 

M  +  HjO  —  OH,*  -k  H*  +  e’  [3] 

The  formation  of  OH,*,  shown  in  Eq.  3  is  the  rate-determining 
step  of  this  reaction  and  occurs  on  platinum  at  potentials  of  0.7  V  vs. 
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the  reversible  hydrogen  electrode  (RHE)  and  above.  ’  ’  Ruthe- 
nium  has  the  ability  to  form  OH^a,  from  water  at  significantly  lower 
potentials  than  Pt,  0.35  V  for  50  atom  %  Ru  and  0.2  V  for  90  atom 
%  Ru,*^’*^’*^  allowing  for  a  certain  amount  of  CO  tolerance.  At  low 
temperatures  (70-85“C),  the  Pt:Ru  (1:1  atomic  ratio)  alloy  catalyst 
provides  nearly  equivalent  performance  to  pure^Hj  on  Pt  for  CO 
concentrations  up  to  100  ppm  in  the  feedstream. 

The  injection  of  oxygen  into  the  fuel  stream  increases  catalyst 
tolerance  to  CO. However,  even  the  addition  of  high  levels  of 
oxygen  to  the  feedstream  (2-4%  by  volume  of  hydrogen)  provides 
approximately  100  ppm  CO  tolerance.  Roughly  one  out  of  every 
400  O2  molecules  oxidizes  an  adsorbed  CO  molecule,  with  the  bal¬ 
ance  reacting  with  hydrogen.'"'  The  placement  of  a  layer  of  Ru  cata¬ 
lyst  before  the  Pt  electrode  to  act  as  a  filter  has  been  shown  to 
increase  the  effectiveness  of  oxygen  addition  over  conventional 
Pt:Ru  alloy  catalysts.'®  This  method  also  eliminates  the  process  of 
alloying  the  Pt  and  Ru  metals,  and  by  using  a  filter,  the  Pt  loading  in 
the  anode  is  free  to  be  varied.  . 

It  is  believed  that  the  following  mechanism  is  occurring  in  the 
Ru  filter  for  a  fiilly  humidified  fuel  stream  containing  oxygen,  hy¬ 
drogen,  and  carbon  monoxide'"'’'®'^’ 

Ru  +  H2  ^  2Hads 

Ru  +  CO  COads 

O2  “b  2Ru  *-*  2  Ogds 

Ha*  +  0,*  -  OH,*  +  Ru 

COad,  +  OH,*  -  CO2  +  H*  +  2Ru  +  e’ 

Had,  +  OH,*  -HI  HjO  +  Ru 


[5] 

[6] 

[7] 

[8] 
[9] 


CO,*+ O,*-!  CO2  +  2Ru  [10] 

H,* -*  H*  +  e"+ Ru  [11] 


Reactions  4,  5,  and  6  represent  the  adsorption  of  species  onto  the 
Ru  catalyst.  Reaction  7  represents  an  intermediate  reaction  on  Ru 
resulting  in  the  formation  of  OH,*.  Reactions  8-11  represent  com¬ 
peting  desorption  reactions  on  the  Ru  catalyst.  Equation  10  is  well 
documented  as  following  Langmuir-Hinshelwood  kinetics  on  Pt  and 
Ru.^"*'^’  Of  the  desorption  reactions,  Reactions  8  and  10  are  desired 
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as  they  result  in  the  oxidation  of  CO., 3-  Earlier  work'*  showed  that 
the  reaction  shown  in  Eq.  8  is  the  primary  means  by  's 

oxidized  within  the  filter,  resulting  m  the  formation  of  "  • 
like  the  Pt  region  following  the  filter,  a  three-phase  interface  of 
catalyst  (Ru  here),  carbon  and  Nafion  are  required  for  the  filter  to 
operate.  The  above  mechanism  is  not  a  comprehensive  I'st  of  al 
reactions  occurring  within  the  anode  region.  The  mechanism  focuses 
on  those  reactions  that  lead  to  CO  oxidation.  For  example,  evidence 
has  been  found  that  HjOj  is  formed  as  an  intermediate  during  oxy¬ 
gen  reduction  at  overpotentials  as  low  as  0.5  V.  however.  HjOj 
breaks  down  to  oxygen-containmg  compounds  at  or  even  before 

reachinc  the  catalyst  surface.  ’  ’  ^ 

The  goal  of  this  work  is  to  increase  the  effectiveness  of  the  Ru 
filter  by  applying  it  through  sputter  deposition  as  opposed  to  con¬ 
ventional  catalyst  ink-based  application  methods.  Sputter  deposition 
is  widely  used  for  integrated  circuit 
investigated  for  us  in  fuel  cells  for  more  than  a  decade. 
ern/.'"’  sputter  deposited  platinum  on  uncatalyzed  gas-diffusion 
layers  (GDLs)  resulting  in  cell  performances  at  loadings  o 
0  10  mg  Pt/cm^  equivalent  to  those  of  standard  methods  at  loadings 
of  0.40mgPt/cm*.  Witham  etnC"  achieved  direct  methanol  fuel 
cell  (DMFC)  anode  catalyst  activities  one  to  two  orders  of  magni¬ 
tude  higher  than  those  of  conventional  ink-based  catalysts,  sugges  - 
ing  that  DMFC  anodes  may  be  manufactured  containing  less  Aan 
one-tenth  the  amounts  presently  used  (2.5-4  mg  P^cm  )  without 
loss  in  performance.  Holleck  el  al*^  have  sputtered  catalyst  mix¬ 
tures  at  the  front  surface  of  the  anode 

included  Pt:Ru;X  where  X  represents  Ni.  Pd,  Co,  Rh,  Ir,  M";  W. 
Nb,  and  Zr.  For  low  levels  of  CO  (10  ppm),  specific  Pf.W  and 
Pt  Ru;W  alloys  performed  better  than  Pt:Ru. 

In  this  work,  the  Ru  filter  is  sputter-deposited  as  a  single  layer 
and  as  a  series  of  three  layers  (separated  by  Nafion-carbon  ink)  and 
compared  to  a  conventional  ink-based  filter.  By  manufactunrig 
membrane-electrode  assemblies  (MEAs)  composed  of  multiple 
sputter-deposited  Pt  layers,  Cha  and  Lee  were  able  '''®^®“®  *® 
Pt  catalyst  activity  significantly  over  conventional  ink-based  MEAs 
by  increasing  the  amount  of  Pt  in  contact  with  Nafion  and  carbon_ 
This  process  was  further  improved  upon  by  reducing  the  amount  of 
Nafion-carbon  ink  (NCI)  separating  the  Pt  layere,  resulting  in  thin¬ 
ner  more  effective  electrodes.''^  For  the  same  loading  of  catalyst, 
spu’tter-depositing  catalyst  between  layers  of  NCI  increases^*e  ac¬ 
tive  area  of  the  catalyst  vs.  a  single  sputter-deposited  layer  For 
Pt  catalyst,  this  results  in  greater  performance  of  the  fuel  ce 
electrodes  It  is  predicted  that  for  the  Ru  filter,  this  process  will 
create  more  sites  upon  which  CO  can  be  oxidized,  resulting  in  a 
more  CO-tolerant  PEMFC  MEA. 

Experimental 

Development  of  ink-based  MEAs.-Tht  method  described  in 
U.S.  Patent  5,21 1,984  provides  an  outline  for  the  catalyst  mk  prepa¬ 
ration  and  MEA  fabrication  performed  in  this  project.  The  follow¬ 
ing  catalyst  inks  were  prepared;  (0  Nafion  -F  carbon  only  (Nafion- 
carbon  Ink  or  NCI),  (fi)  Nafion  -F  20%  Pt  on  carbon  (11.) 
Nafion  +  20%  Pt:Ru  on  carbon,  and  (iv)  Nafion  +  20/o  Ru  on 

‘'“'The  inks  were  prepared  for  Pt,  PtiRu,  Ru  by  adding  ‘he  E-TEK 
catalyst  (20%  catalyst  on  XC-72  Carbon)  to  a  solution  of  5  wt  /o 
Nafion  (DuPont).  For  NCI,  XC-72  Carbon  was  added  to  a  solution 

of5wl%  Nafion.  ...  *  r 

Decals  (Teflon,  10  cm^  three-ply)  were  weighed  prior  to  appli- 
cation  of  catalyst  ink.  The  ink  was  drawn  across  the  surface  of  the 
decals  using  a  Meyer  rod.  The  coated  decals  were  dned  in  an  oven 
at  105°C  under  ambient  pressure  for  10  mm. 

Target  loadings  were  0.15mgPt/cm^  (anode),  0.230  mg 
Pt;Ru/cm*  (anode),  and  0.15mgPt/cm*  (cathode).  Target  loadings 
of  0.230  mg  PtiRu/cm*  and  of  0.080  mg  Ru/cm^  for  the  Ru  filter 
were  chosen  so  that  a  1 : 1  atomic  ratio  of  PtiRu  is  mamtamed,  while 
Pt  catalyst  is  maintained  at  a  loading  of  0.15  mg/cm  . 


Table  1.  Types  of  ink-based  MEUs  tested. 

MEU  name  Description 


Pt 


Pt:Ru 


0.15  mg  Pt/cm^  anode 

0.15  mg  Pt/cm^  cathode 

Uncatalyzed  Toray  GDLs  on  anode  and  cathode 

0.23  mg  Pt:Ru/cm^  anode 
0.15  mg  Pt/cm^  cathode 

Uncaialyzed  Toray  GDLs  on  anode  and  cathode 


Pt  +  Ru  (Ru  filter)  0.15  mg  Pt'em^  anode 
0.15  mg  Pt/enr  cathode 

0.080  mg  Ru/cm“  coated  Toray  GDL  on  the  anode 
Uncatalyzed  Toray  GDL  on  the  cathode  side 


To  form  a  MEA  with  ink-coated  decals,  appropriate  decals  were 
placed  on  either  side  of  the  PEM  (Nafion  117,  protonated  fonn). 
This  assembly  was  hot-pressed  to  ensure  bonding.  It  was  then 
cooled  to  room  temperature,  before  the  decals  were  carefully  peal 

from  the  assembly.  ,  .  ^ 

To  form  a  membrane-electrode  unit  (MEU,  equivalent  to 
GDL  +  MEA  +  GDL)  with  ink-coated  GDLs,  appropriate  GDLs 
were  placed  on  either  side  of  a  PEM  (Nafion  117,  protonated  form 
soaked  in  deionized  water  for  I  h).  This  assembly  was  hot-pressed 
to  ensure  a  well-bonded  MEU.  The  types  of  ink-based  MEUs  made 
are  shown  in  Table  I. 

Development  of  sputter-deposition-based  A/£'.4i.--Plasma  modi¬ 
fications  and  sputter-deposition  augmentations/additions  were  both 
completed  using  an  Anatech  Hummer  10.2  sputter-coating  tool.  A 
modified  sample  stage  was  used  to  support  Nafion  1 17  substrates  up 
to  6  X  6  cm  while  masking  1.5  cm  about  the  membrane’s  penmeter. 

It  has  been  shown  previously  that  alternating  current  (ac)  plasma 
should  be  used  for  all  plasma  modifications,  as  Nafion  17  and 
MEAs  from  Nafion  117  suffer  no  ill  effects  from  this  ‘reatmem. 

An  aluminum  target  was  used  for  ac  plasma  modifications,  while  the 
Ru  (Kurt  J.  Lesker)  target  and  carbon  evaporation  system  (Anatech) 
were  used  for  sputter  augmentations/additions. 

All  MEAs  subject  to  sputter  deposition  were  first  ac  Plasma 
cleaned  for  5  min  at  5  mA  and  1.2  kV  to  remove  ^u.  dup 

from  the  target  as  well  as  roughen  the  substrate  surface.  All  tre 
ments  were  completed  at  ~62  mTorr.  A  separate  vacuum  chamber 
was  used  to  evacuate  each  substrate  to  ~45  mTorr,  before  it  was 
placed  in  the  sputter-coating  tool  to  minimize  contaminant  outgas- 
Lg  in  the  depLition  chamber.  A  potential  of  1.8  kV  ^"<1  a  currem 
of  8  mA  were  maintained  to  control  the  deposition  rate  for  Ru.  A 
SiO,  sample  was  sputter-deposited  in  situ  with  each  MEA.  The  re¬ 
sultant  metal/SiO,  stack  was  subjected  to  "r 

ning  electron  microscopy  (SEM)  imaging  to  verify  the  thickness  of 
the  sputter-deposited  film  and  top-view  SEM  imaging  to  determine 

surface  characteristics  of  the  film.  ,  .  •  -i 

Multilayered  Ru  filters  were  manufactured  by  a  method  similar 
to  that  described  by  Cha  and  Lee.«  On  the  MEA  anode,  a  layer  of 
NCI  was  sprayed  on  the  surface  of  the  sputter-deposited  Ru  layer. 
The  MEA  was  then  dried  at  80‘’C  under  vacuum  for  10  mm.  Addi¬ 
tional  layers  of  sputter-deposited  Ru  and  NCI  were  added  until  th 
desired  Ru  loading  and  number  of  Ru  layers  were  attained. 

MEAs  containing  sputter-deposited  catal^yst  "'®.'’®  . 

factored  into  MEUs  through  similar  methods  as  the  ink-based 

MEAs. 

Cell  assembly  and  lesting.-The  MEUs  were  placed  in  a  10  cm* 
cell  assembly  and  incubated  for  4  to  8  h  at  ambient  pressure 
cell  temperature  of  10°C,  stoichiometric  ratio  ([actual  flow]/ 
[stoichiometric  flow]  required  for  a  1.0  A/cm  current)  a‘  « 

anode  and  2.0  at  the  cathode.  Fuel  cell  performance  curves  were 
obtained  under  the  conditions  given  in  Table  11. 


A870 


Journal  of  The  Electrochemical  Societ}\  149  (7)  A868-A872  (2002) 


Table  11.  Fuel  cell  test  conditions. 

Pressure  1  ^tm 

Cell  temperature  70°C 

Stoichiometric  ratio  (at  1  A/cm^)  1-5  hydrogen 
2.0  air 

Anode:  hydrogen, 

H2  +  CO.  H2  +  CO  +  air  bleed 
Cathode:  air 

Complete  humidification  of  anode 
and  cathode  gas  streams  for  all  trials 

2.0%  O2  (in  the  form  of  an  air  bleed 
relative  to  the  volumetric  flow  of 
hydrogen  in  slm) 

50,  200  ppm 


Results 

Determination  of  sputter-deposition  rates  and  catalyst  load- 
sputter-deposition  rates  for  Ru  were  determined  with 
loadings  calculated  from  cross-sectional  and  top-view  SEM  images 
of  sputter-deposited  Ru  films  on  Si02  substrates.  The  top-view  SEM 
images  in  Fig.  I  were  analyzed  to  determine  the  surface  coverage  of 
the  sputter-deposited  film  (Pt  65%,  Ru  --  56%).  Ru  did  not  form 
a  continuous  film  on  the  Si02  substrate,  but  rather  agglomerated. 
This  is  consistent  with  literature The  surface  coverages  were  used 
in  conjunction  with  the  bulk  density  of  Ru  (12.2  g/mL)  and  the  film 
thicknesses  from  the  cross-sectional  SEM  images  to  calculate  the 
subsequent  Ru  loadings.  The  Ru  sputter-deposition  rate  was  constant 
with  time  at  3.3  p.g  Ru/cm^/min. 

The  cross-sectional  view  in  Fig.  1  shows  that  the  thickness  of  30 
and  45  min  sputter-deposited  Ru  is  roughly  two  and  three  times  the 
thickness  of  a  15  min  deposition  of  Ru,  respectively.  Semi  quantita¬ 
tive  analyses  via  energy  dispersive  X-ray  spectroscopy  (EDXS)  and 
Rutherford  backscattering  spectrometry  (RBS)  confirm  these  Pt  and 
Ru  deposition  rates. 

CO  testing.^Jhc  three  types  of  Pt  +  Ru  filter  anodes  prepared 
are  shown  in  Fig.  2.  They  were  made  to  the  following  specifications: 
filter  1:  0.08  mg/cm^  ink-based  20%  Ru/C.  filter  2:  NCI  +  25  min 
(0.08  mg/cm^)  of  sputter-deposited  Ru,  and  filter  3:  NCI  +  3 
X  (8,33  min  of  sputter-deposited  Ru  +  NCI),  The  total  Ru  loading 
was  0.08  mg/cm^. 

Filter  1  is  an  ink-based  Ru  filter,  while  filter  2  and  3  are  com¬ 
prised  of  one  and  three  layers  of  a  sputter-deposited  Ru,  respec- 


Fccdstrcams 

Humidification 
Air  bleed 

CO  amounts 


Ru  Filter  #  I  Ru  Filter  #2  Ru  Filter  #3 


Figure  2.  Diagram  of  the  Ru  filtered  anode.  The  filter  is  a  layer  of  Ru 
catalyst  placed  between  the  Pt  catalyst  and  the  gas-diffiision  layer  for  o.xi- 
dizing  CO  present  within  the  anode.  Filter  1  is  the  standard  Ru  filter  prepared 
from^a  catalyst  ink.  Filter  2  is  a  single  sputter-deposited  layer  of  Ru  sepa¬ 
rated  from  the  Pt  electrode  by  a  layer  of  NCI.  Filter  3  is  a  three-layer 
sputter-deposited  Ru  filter  separated  from  the  Pt  electrode  by  a  layer  of  NCI. 
All  filters  have  a  loading  of  0.08  mg  Ru,/cm^  and  were  placed  atop  a  Pt 
electrode  with  a  loading  of  0.15  mg  Pt/cm^ 


lively.  As  shown  in  Fig.  2,  NCI  is  used  to  separate  the  layers  of 
sputter-deposited  Ru  layered  between  Nafion-carbon  ink  (NCI). 

Figure  3  compares  various  Ru  filters  to  the  Pt:Ru  alloy  and  the 
baseline  MEA  at  70°C  for  an  anode  feed  of  hydrogen  +  50  ppm 
CO.  All  MEAs  prepared  with  Pt,  Pt:Ru  alloy,  and  Pt  +  Ru  filters 
1-3  e.xhibited  a  dramatic  loss  of  performance  when  50  ppm  CO  was 
added  to  the  hydrogen  feedstream.  In  all  cases,  fuel  cell  performance 
dropped  to  less  than  40%  of  the  baseline  MEA  running  on  pure 
hydrogen  fuel.  It  is  clear  that  the  filters  do  not  completely  oxidize 
the  CO  in  the  feedstream  to  CO2 .  The  remaining  CO  passes  through 
the  filter  and  poisons  the  Pt  portion  of  the  anode,  resulting  in  per¬ 
formance  resembling  that  of  a  typical  Pt  anode  as  evidenced  in  Fig. 
3.  The  Pt:Ru  alloy  demonstrates  CO  tolerance  over  a  Pt  anode  con¬ 
sistent  with  literature,  more  than  doubling  the  current  density  of  any 
other  anode  configuration  at  0.6  V. 

When  2%  O2  (in  the  form  of  an  air  bleed)  is  added  to  the  anode 
feed,  all  Ru-containing  anodes  show  significant  performance  irn- 
provement  over  the  Pt  baseline  as  shown  in  Fig.  4.  For  the  25  min 
Ru  (0.08  mg  Ru/cm^)  sputter-deposited  filter,  there  is  almost  no  per¬ 
formance  loss  except  at  lower  voltages.  The  Pt:Ru  alloy  and  stan¬ 
dard  Ru  filter  show  almost  identical  performance,  roughly  20%  less 
than  the  baseline  (0.227  V5.  0.275  A/cm^  at  0.6  V).  The  performance 
of  the  Ru  filter  cannot  be  explained  by  the  fonmation  of  a  Pt:Ru 
alloy  as  there  is  no  interface  between  Pt  and  Ru.  A  layer  of  NCI 
separates  the  sputter-deposited  Ru  from  the  Pt  electrode  for  the 
sputter-deposited  Ru  filters.  The  multilayered  filter  performed  the 


45  minutes  Ru  30  minutes  Ru  15  minutes  Ru 


Figure  1.  Top  and  cross-sectional  SEM 
images  of  sputter-deposited  Ru  on 
Si/SiOj  substrates. 
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Currant  Density  (A/cm*) 

Figure  3.  Performance  comparison  of  MEUs  with  sputter-deposited  Ru  filter 
under  H2  +  50  ppm  CO  conditions.  P  =  1  atm,  T  =  70°C.  The  filter  con¬ 
figurations  are  defined  as  follows:  filter  I  (—A—)  is  the  standard  Ru  filter 
prepared  from  an  ink;  filter  2  (— O--)  is  a  single  25  min  sputter- deposited 
layer  of  Ru  separated  from  the  Pt  electrode  by  a  layer  of  NCI;  filter  3 

( _ P — )  is  a  three-layer  sputter-deposited  Ru  filter  (8.33  min  per  layer,  see 

Fig.  2  for  design)  separated  from  the  Pt  electrode  by  a  layer  of  NCI. 


worst  of  the  Ru-containing  anodes.  This  is  likely  due  to  the  in¬ 
creased  diffiisional  resistances  caused  by  the  thickness  of  the  Ru 
filter. 

As  the  amount  of  CO  is  increased  from  50  to  200  ppm,  the 
benefit  of  the  Ru  filters  is  seen  more  clearly.  Figure  5  shows  that  all 
three  Ru  filter  types  outperform  the  Pt:Ru  alloy.  For  the  anode  con¬ 
sisting  of  Pt  +  Ru  filter  2,  the  performance  is  double  that  of  the 
Pt:Ru  alloy  (0.205  V5.  0.103  A/cm^  at  0.6  V).  The  MEA  containing 
filter  2  loses  only  20%  of  its  performance  (0.255  vs.  0.205  A/cm^  at 
0.6  V)  as  the  CO  concentration  is  increased  from  50  to  200  ppm  CO, 
and  loses  only  25%  V5.  the  hydrogen  baseline.  There  have  been  other 
reports  of  200  ppm  CO  tolerance  with  an  air  bleed,  but  the  condi¬ 
tions  involved  higher  temperatures  and  higher  loadings.  The  MEA 
containing  filter  2  significantly  outperformed  ink-based  filter  1  by 
more  than  35%  (0.205  V5.  0.149  A/cm^  at  0.6  V).  These  results  also 
constitute  a  significant  improvement  over  the  previous  work  done  on 
the  Ru  filter. ^^Under  conditions  identical  to  those  shown  in  Table  II, 
the  Ru  filter  of  0,21  mg/cm^  developed  in  Ref.  18  showed  an 
equivalent  loss  in  performance  to  filter  2  of  this  work  when  only  100 
ppm  CO  was  added  to  the  anode  feed  containing  2%  O2  in  the  form 


Currant  Danalty  (A/cm*) 

Figure  4.  Performance  comparison  of  MEUs  with  sputter-deposited  Ru  filter 
for  an  anode  feed  of  hydrogen  +  50  ppm  CO  +  2%  O2 .  P  -  I  atm, 
T  =  70®C.  The  filter  configurations  are  defined  in  Fig.  3. 


Figure  5.  Performance  comparison  of  MEUs  with  sputter-deposited  Ru  filter 
for  an  anode  feed  of  hydrogen  +  200  ppm  CO  +2%  O2 .  F  =  I  atm, 
T  =  70°C.  The  filter  configurations  are  defined  in  Fig.  3. 


of  an  air  bleed.  The  Ru  filter  in  Ref.  18  showed  a  50%  loss  in 
performance  when  the  CO  concentration  was  doubled  to  200  ppm. 

Of  the  three  Ru-filtered  MEAs  tested,  filter  3  performed  the 
worst  at  200  ppm  CO,  but  suffered  only  a  12%  loss  compared  to  Us 
performance  at  50  ppm  CO  (0.138  V5.  0.157  A/cm^  at  0.6  V).  That 
there  is  so  little  drop-off  from  this  three-layer  filter  configuration 
suggests  that  the  lower  performance  compared  to  the  single-layer 
spuUer-deposited  filter  2  is  due  to  diffxisional  resistances  caused  by 
the  increased  thickness  of  this  filter.  The  layers  of  NCI  are  very  thick 
relative  to  each  Ru  deposition  (-^40  nm  for  each  spuUer-deposited 
Ru  layer  compared  to  -12  jxm  for  each  NCI  layer).^  accounting 
‘for  more  than  99.5%  of  the  Ru  filter  layer  thickness.  The  increased 
thickness  of  this  three-layer  filter  requires  hydrogen  to  travel  further 
to  reach  the  Pt  vs.  a  single-layer  filter,  and  vs.  the  baseline  MEA 
giving  rise  to  greater  diffusional  resistances.  By  reducing  the 
amount  of  NCI  used  in  each  sputter-deposited  layer  (by  diluting 
the  NCI),  a  thinner,  more  effective  multilayered  Ru  filter  may  be 
developed. 

The  addition  of  a  layer  of  NCI  between  the  Pt  electrode  arid  the 
sputter-deposited  Ru  filter  is  evidence  that  the  CO  oxidation  is  not 
the  result  of  a  layer  of  Pt:Ru  alloy  being  formed  at  the  interface  of 
the  filter  and  the  Pt  electrode.  Thus,  it  can  be  hypothesized  that  the 
CO  oxidation  is  occurring  very  fast  in  relation  to  the  dilTusion  of  the 
gases  through  the  filter  region.  Nearly  all  the  CO  is  oxidized  within 
a  region  -100  to  120  nm  thick  (see  Fig.  1).  Therefore,  it  is  assumed 
that  a  filter  of  identical  loading  (0.08  mg  Ru/cm^)  would  perform 
similarly  in  front  of  a  Pt  electrode  of  any  loading. 

It  is  assumed  that  all  oxygen  reacts  almost  immediately  within 
the  filter  region  and  that  the  hydroxyl  groups  produced  then  react 
with  CO  to  produce  CO2 ,  Even  with  200  ppm  CO  tolerance,  only 
one  out  200  O2  molecules  is  being  used  to  oxidize  CO.  Further 
research  into  catalysts  with  lower  activity  than  Ru  with  respect  to 
hydrogen  may  provide  even  further  CO  oxidation  by  further  shifting 
the  selectivity  of  the  oxygen  reactions  toward  the  oxidation  of  CO  to 

CO2. 

Conclusions 

The  sputter-deposited  Ru  filters  developed  here,  containing  less 
than  40%  the  amount  ofRu  (0.080  vi.  0.21  mg/cm^)  of  the  previous 
work,'*  achieved  twice  the  CO  tolerance  (200  ppm  CO  vs.  100  ppm 
CO)  under  similar  operating  conditions  (Hi  +  2%  0,  anode  feed,  a 
cathode  feed  of  air,  70‘>C  operating  temperature). 

For  an  anode  feedstream  consisting  of  hydrogen,  200  ppm  CO, 
and  2%  oxygen,  all  MEAs  containing  Pt  +  Ru  filter  anodes  showed 
increased  CO  tolerance  compared  to  a  Pt:Ru  alloy  containing  similar 
amounts  of  Pt  and  Ru  (0.150  and  0.080  mgW,  respectively).  The 
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MEA  containing  a  single  sputter-deposited  Ru  filter  exhibited 
greater  CO  tolerance  than  that  of  the  MEA  containing  the  ink-based 
Ru  filter  (0.205  vs.  0.149  A/cm^  at  0.6  V). 

Attaining  200  ppm  CO  tolerance  by  using  2%  oxygen  is  equiva¬ 
lent  to  one  out  200  O2  molecules  being  used  to  oxidize  COjjj  mol¬ 
ecule,  while  the  rest  likely  react  with  to  form  water.  Catalysts 
with  characteristics  similar  to  Ru  but  with  lower  hydrogen  activity 
would  likely  allow  for  a  higher  percentage  of  oxygen  in  the  feed- 
stream  to  react  with  adsorbed  CO,  resulting  in  a  more  effective  filter. 

While  the  three-layered  sputter-deposited  Ru  filter  (filter  3)  per¬ 
formed  the  worst  of  the  three  filters  its  drop  in  Performance  (12  /^ 
as  the  CO  in  the  anode  feed  was  increased  from  50  to  200  ppm 
(balance  H,  +  an  air  bleed  containing  2%  Oj)  was  the  smallest  of 
all  the  anodes  tested.  NCI  accounts  for  99.7%  of  the  multilayered 
Ru  filter  thickness.  Further  research  can  be  done  to  optimize  this 
three-phase  interface  area  of  the  Ru  filter  and  eliminate  the  unused 
portion  of  the  electrode,  resulting  in  thinner,  more  effective  filtere. 
However,  to  generate  a  filter  with  a  high  number  of  Ru  layers  is 
neither  the  most  economical  (due  to  the  time  required  to  generate 
such  a  multilayered  MEA)  nor  most  effective  approach.  The  genera¬ 
tion  of  a  continuous  three-phase  interface  is  the  ultimate  goal  of  this 
method  and  this  is  what  should  be  pursued  using  the  method  of 
sputter  deposition.  Simultaneously  sputter-depositing  Ru  and  spray 
depositing  NCI  could  produce  a  continuous  three-phase  interface 
region  This  would  result  in  an  extremely  thin  (~1  p.m)  filter  that 
may  provide  high  Ru  activity  with  minimal  diffusional  resistances. 
And  by  applying  this  Ru/C/Nafion  filter  in  a  single  application  to  the 
PEM  the  process  is  less  time-consuming  and  thus  more  economical. 

Because  benefits  of  the  Ru  filter  occur  at  high  level  air  bleed  (2^ 
0,)  and  the  Pt:Ru  alloy  provides  CO  tolerance  even  without  air 
bleed  it  is  suggested  that  the  anode  configuration  that  would  pro¬ 
vide  optimal  CO  tolerance  would  consist  of  a  sputter-deposited  Ru 
filter  placed  in  front  of  and  adjacent  to  a  Pt:Ru  alloy. 
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1.  Introduction 

Industrial  chemical  reactions  are  usually  accompanied  ^ 
with  mass  and  energy  transfer,  either  homogeneously  or 
heterogeneously.  Mathematical  modeling  for  th^e  pro-  r 
cesses  is  based  on  material  and  energy  balance.  One  can  b 
generate  a  set  of  differential  equations  known  as  the  reac- 
rion-<iiffusion  problem.  Owing  to  the  strong  nonlinearity  of  u 
the  reaction  rate,  mainly  from  the  effect  of  temperature, 
reaction-diffusion  equations  are  paid  more  auention  in 
analyzing  and  designing  chemical  and  catalytic  reactors 
and  are  the  major  role  in  analyzing  the  nonlinear  dynamic  e 
behaviors  in  reactor  engineering.  The  same  phenomena  exis 
in  electrochemical  processes,  with  the  add  complexity  of  a  r 
varying  potential  field,  and  considerable  research  has  been 
reviewed  for  electrochemical  reactions  occumng  in  the  c 

porous  electrode  [1].  r  v  j  \ 

Newman  discussed  the  numerical  solution  of  boundai7- 
value  problems  consisting  coupled  ordinary  differentia  ( 
equations  which  one  can  often  met  in  chemical  engineenng 
science,  and  developed  a  unique  technique  to  solve  coupled, 
linear  differential  equations  [2-4].  In  his  procedure,  a  larg^  ■ 
sparse  matrix  was  collapsed  to  a  tridiagonal  matnx  with 
elements  of  block  matrixes,  which  make  it  easy  to  be 
inversed  in  solving  the  algebraic  equations  transformed  from 

•Coaesponding  author.  Tel.:  +1-803-777-3207;  fax;  +1-803-777-8265. 
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the  original  differential  equations.  Accordingly,  subrou¬ 
tine  in  his  numerical  method  is  called  Newman  s  Band(/). 
White  discussed  Newman’s  method,  and  promoted  this 
technique  for  wide  application  [5-7],  However,  for  solving 
reaction-diffusion  equations— the  nonlinear  two-poin 
boundary  value  problem,  respective 
rithms  are  needed.  Shooting  methods,  such  as  Runge  Kutta 
integration  scheme,  is  a  marching  technique  that  must  be 
backwards  processed  from  end  point  to  beginning  point  to 
avoid  inherent  instability.  This  method  is  coinmonly  used 
but  it  depends  on  a  proper  choice  of  the  initial  value  at  the 
ending  point.  Linearizing  the  nonlinear  kinetics  term  that 
convert  the  nonlinear  differential  equations  into  a  line^  one 
is  called  the  linearization  trial-and-error  technique,  but  i 
does  not  guarantee  that  all  steady  state  solutions  can  be 
determined.  Compared  with  these  two  methods  transform¬ 
ing  the  nonlinear  differential  equations  into  nonlinear  alge¬ 
braic  equations  that  are  then  solved  by  the  Newton  method 
did  not  receive  the  appropriate  attention  because  the  ca  cu¬ 
tting  time  and  storage  is  greater.  ^ 

In  this  paper,  we  extend  Newmans’s  mamx  method  to 
solve  coupled  electrochemical  reaction-diffusion  equa- 
tions-a  coupled  high  nonlinear  two-point  boundary  value 
problem  that  is  encountered  in  modeling  of  a  fue 
catalyst  layer.  Based  on  Taylor  expansion,  a  tnal  and  etro 
iteraction  algorithm  is  developed  for  solving  the  nonline^ 
algebraic  equations  transformed  from  the  electrochemical 
reaction-diffusion  equations.  Owing  to  the  adoption  o 
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Nomenclature 

specific  interfacial  area  (cm  ) 
dimensionless  form  of  Ch+  (co  = 
dimensionless  form  of  Chi  (cj  = 
dimensionless  electric  potential  in  the  solid 
matrix  phase  {c2—nf4>\)  .  . 

dimensionless  electric  potential  in  the  pore  ionic 

phase  (C3  =  nf  <(>2)  . 

effective  concentration  of  proton  per  unit  volume 

of  solution  (mol  cm"^) 

effective  concentration  of  hydrogen  per  unit 
volume  of  solution  (mol  cm"  ) 
initial  effective  concentration  of  hydrogen 


1  Co 
Cl 
C2 

C3 


^2 

\Dh: 

/ 

F 


(mol  cm"^) 

effective  diffusion  coefficient  of  hydrogen  in  pore 
phase  (cm^  s"‘) 
formal  potential  (V) 
constant  (/  =  FIRT  (V  )) 

Faraday’s  constant  (96485  C  eq.  ) 

total  current  density  leaving  the  matrix  phase 


'2 


k2 

I 


Nh, 


(A  cm"  )  .  . 

superficial  electronic  current  density  in  the 

matrix  (A  cm"^) 

superficial  ionic  current  density  in  pore  phase 
(A  cm"^) 

current  density  of  fuel  cell,  defined  as  positive 

(Acm"^)  ..  . 

k°  heterogeneous  rate  constant  for  hydrogen  oxidation 

(cms"') 

dimensionless  parameter  {ki  =  nfllcew/c) 
dimensionless  parameter  (k2  =  nfllcew/K) 

,  thickness  of  anode  catalyst  layer  (cm) 
n  =  2  number  of  electrons  transferred  in  electrode 
reaction  (=2) 

superficial  flux  density  of  hydrogen 
(molcm"^s"')  -1  ir-K 

1  i?  universal  gas  constant  (8.314  J  mol  K  ) 

electrode  reaction  rate  of  hydrogen,  defined  as 
positive  (mol  cm  s  ) 

T  absolute  temperature  (K) 

U  dimensionless  formal  potential  {U  =  nfEr) 

X  dimensionless  form  of  distance  through  porous 
electrode  (x  =  X/l) 

X  distance  through  porous  electrode  (cm) 

Greek  symbols 
a  transfer  coefficient 

(j)^  electric  potential  in  the  solid  matrix  phase  (V) 

(^2  electric  potential  in  the  pore  ionic  phase  (V) 

0  Thiele  modulus  (<P^  =  a/^^V£>H2) 

$1  dimensionless  parameter  ($,  —  anfFk  I  C^jc) 
02  dimensionless  parameter  (<1>|  =  at^^k^i^C^jK) 
K  effective  conductivity  of  the  pore  ionic  (proton) 

phase  (n"‘  cm"*)  , 

effective  conductivity  ofthe  solid  matrix  (il  cm _ ^ 


Newman’s  unique  method,  in  each  iteration  loop  a  compact 
forward  substitution  can  be  processed.  In  each  forward  or 
backward  substitution,  we  use  Gauss-Jordan  full  pivoting 
technique  to  guarantee  the  numerical  stability.  The  who  e 
computing  is  very  fast,  even  though  we  set  2000  node  points 
to  make  the  substitution  effectively.  For  a  general  case  of 
an  electrochemical  reaction-diffusion  process  in  a  fuel  cell 
catalyst  layer,  just  two  iteration  loops  can  converge  the 
solution.  Using  this  method  to  study  the  effect  of  Thiele 
modulus  on  the  electrochemical  reaction-diffusion  process, 
some  intrinsic  phenomena  due  to  the  nonlinear  behavior  o 
the  system  is  found  by  numerical  solution  that  is  beneficial  to 
analyze  and  scale-up  a  fuel  cell  system. 


2.  Model  equation 


We  consider  for  the  isothermal  hydrogen  oxidation  reac¬ 
tion  occurring  in  a  porous  catalyst  layer  of  a  fuel  cell  anode. 

H2^2H+-l-2e" 

A  schematic  representation  of  it  is  shown  in  Fig.  1.  The 
catalyst  layer  is  viewed  as  a  continuum  of  two  phases,  each 
phase  either  a  pure  ionic  or  electronic  conductor.  The  lomc 
phase  is  considered  a  cationic  selective  polymer,  such  as 
Nafion.  Therefore,  the  proton  concentration  throughout  the 
catalyst  layer  is  fixed.  This  uniform  electrolyte  concentra¬ 
tion  means  the  superficial  current  density  in  the  pore  ionic 
phase  is  due  only  to  ion  migration.  This  is  represented 
mathematically  as 


The  rate  at  which  ionic  current  enters  the  pore  solution  di2/ 
dX  is  proportional  to  the  reaction  rate  on  a  volumetnc  basis 
that  is  expressed  by  Butler-Volmer  expression  (defining 
anodic  current  as  positive)  [8] 

^  =  nFRn,  = 
dX  ' 

(3) 

The  movement  of  electrons  in  the  solid  matrix  phase  of  the 
porous  electrode  is  governed  by  Ohm’s  law 


The  consequence  of  electroneutrality  is  that  the  divergence 
of  the  total  current  density  is  zero. 


(5) 

dX~dX'^dX 

The  flux  of  dissolved  hydrogen  in  the  porous  anode  is 
determined  by  diffusion 


dCH, 


(6) 
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Current  Collector 
Porous  Catalytic  Anode 

Membrane 

Porous  Catalytic  Cathode 

Current  Collector 

Fig.  1.  Hydrogen  oxidation  in  a  fuel  cell  anode. 


Material  balance  for  this  reaction-diffusion  process  c^  be 

modeled  according  to  equaUonofcontinuity  at  steady  stat  [  ] 

cWh: 


=  a  —  (Thiele  modulus) ,  V  =  n/E®' , 


dX 


=  Rhi 


(7)  <^]=- 


an 


H2 


4>l  = 


an^fFk^i^C^H, 


(8) 


dX2 


Substituting  Eq.  (3).  (6)  into  Eq.  (7)  yields 

Differentiating  Eqs.  (2)  and  (4)  and  substituting  Eq.  (3)  into 
the  resulting  equation  gives 

AV2  ^ 

Boundary  conditions: 

1.  @X==0, 

Ch:  =  r^2 

d^j  fcell 

"dX  ^  “  ff 

^2  =  0 

2.  @  X  =  /, 

dX 
d(^i 
dX 

d<l>2  _  fcell 

dX 


!k  =  o 

^  =  0 


(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 


(17) 

(18) 

(19) 

a  ^ 

Then  Eqs.  (8).  (10)  and  (1 1)  are  reduced  to  the  dimensionless 
form 


ffi  = 

dbc- 

dx^ 

dx- 

Boundary  conditions*. 

1.  @  x^O, 

c\  =  1 

dx 

a  =  0 

2.  @  2:=  1, 

dci 


dr 


=  0 


!fi  =  0 

dc 


dC3 

(be 


=  -k2 


(20) 

(21) 

.(22) 

(23) 

(24) 

(25) 


Setting  the  following  dimensionless  variables  and  parameters 
r*,.  ^ 


Ch+ 
^0  =  70“ 


Ch. 

^H2 


C2  =  n/(/>i,  C3— n/</>2f 


Although  the  above  model  Eqs.  (17)  (25)  are 
isothermal  system,  due  to  the  unique  nonlinear  character¬ 
istics  of  the^lectrochemical  reaction  term  that  is  similar 
to  the  effect  of  temperature  on  a  chemical  reaction  ra  e, 
r  solvinl  me-hod  tor 

reaction-riiffusion  equations  can  be  transformed  to  study 
non-isothcrmal  problems. 
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3.  Numerical  solution 

For  solving  the  model  Eqs.  (17)-'(25),  they  are  cast  into 
sets  of  nonlinear  algebraic  equations  approximated  by 
finite  difference.  Then  an  iterative  algorithm  is  constructed 
based  on  the  Taylor  expansion.  However,  if  one  adds  the 
Eqs.  (17H19)  together  to  try  to  eliminate  the  nonlinear 
term,  the  computation  is  not  stable.  Analyzing  the  com¬ 
puting  process  on  a  microcomputer  with  eight  word  bits 
for  a  double  precision  real  showed  that  the  block 
matrixes  were  ill-conditioned.  The  condition  numbers 
are  too  large  to  inverse  the  coefficient  block  matrixes 
correctly  due  to  the  error  accumulation  even  for  the  second 
node  point. 

Therefore,  at  each  node  point  except  the  boundary  points, 
the  second-order  derivatives  can  be  approximated  by  three 
point  central  difference  accurate  to  0{h^).  For  example, 
Eq.  (17)  is  approximated  as 

ci{j-  l)-‘2ci{j)+C]{j+  1) 

(26) 

If  there  are  enough  node  points,  we  can  suppose  that  the 
difference  equation  at  node  pointy,  such  as  Eq.  (26),  depends 


only  on  the  variables  at  three  points.  If  the  left-hand  term 
is  expressed  as  Fi,j<  then  the  Taylor  series  expands  it  to 
accuracy  0(h^)  as 


=  c°(y-l)-2c°(y)-)-c?{;+l) 

I  ,  f_  r  ..0/ 


+  h(i)-^SO)l|^L-o 


(27) 


Substituting  corresponding  partial  derivatives  and  intro¬ 
ducing 


ro  = 

rco  =  (1  -a)c?(/>('-“)1^5W-c?(/)-£/l  +acoe-“>^5(/)-5W-yi 

(28) 


it  can  be  simplified  as  a  more  compact  matrix  form. 


;cu-.\  (‘'A 

[1  0  0)|  C2y_i  I  -i-[-(2-fh^<fVo)  ]  I  C2,y 

VijJ 


-f(l  0  0] 


C\,J+l 

C2J+1 

Vc3,;+1, 


=  [-h2<?Vo  -h^‘P^rc, 


(‘h\ 

vIjJ 


+  0<j<N 


(29) 


In  the  same  way,  model  Eqs.  (18)  and  (19)  can  be  approximated  as  Eqs.  (30)  and  (31),  respectively. 


[0  1  0] 


7-1 


ci.7 


+[-h^^iro  -(2-f h'^^]rcj  C2,j  +(0  1  0] 


C2 


\^3,7-1, 


=  [-A2^i>o  -hH\rc,  hH]rJ 


(<j\ 

VIjJ 


.‘^3,7, 


+  h^'J>]rjo,  0<j<N 


7+1 


^<^1, 
C2,7+l 
\^3,7+i, 


(30) 


[0  0  1]|  C2,;-i 

.<^3,7-1/ 


+  [h'^<Plro  hH\rc,  -{p.  +  h^^W)] 


C2.S 

VxjJ 


-f[0  0  1) 


/  Cl, 7+1 
C2,7+l 
\^3,7+1. 


=  [h'^^lro  hH\rc,  -hH\r,,\ 


l<,\ 

Vh) 


-h^^lrjo,  0<j<N, 


(31) 
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Letting 

CU)  =  i<^iJ  <^2.^  0<j<N 

then  Eqs.  (29H31)  become 

AU)CU  -  1)  +  BU)C(j)  +  DU)C{j  +  1)  =  G(y), 

0<j<N, 


(32) 


(33) 


Eqs.  (33),  (35)  and  (36)  can  be  rewritten  as  a  matrix  form 

•B(0)  D{0)  X 
A(l)  B(l)  I)(l) 

A(J)  B(j)  D{j) 


in  which 


50')  = 


(-{2  +  h^^'^ro) 

/i^^fro 

/I  0  0 


-{2+hH]rco) 

-{2-Vh^^\rco)j 


h^4>or, 


'2^co 


-hH'^rco 

h^^\i'co 

\ 


-h^^lrco . 


(c\j\ 

^Zj 

V^J 


(34) 


D{j)  =  I  0  1 
\0  0 

/-h^^'^ro 

Gij)  =  -h^^]ro 

\  h'^'Plro 

(  h-^$^rjo  . 

+ 

Two-point  boundary  conditions  are  approximated  by 
three  point  forward  or  backward  difference  accurate  to 
OQ?)  and  are  also  expressed  as  a  compact  matrix  form 

B(0)C(0)  +  D(0)C(1)  +  XC(2)  =  G(0),  y  =  0  (35) 

and 

YC{N  -2)+  A{N)C{N  -  1)  +  B{N)C{N)  =  G{N), 

j  =  N 


D{0)  = 


(36) 


G(0)  = 


A(Af)  = 


G{N)  = 


(37) 


(38) 


A{N-l)  B{N-l)  D{N-l) 
Y  A{N)  B{N) 


■  C(0)  ■ 

G(0)  ■ 

C(1) 

G(l) 

CO) 

= 

GO) 

C{N-l) 

G{N-l) 

.  C{N)  . 

.  GiN)  . 

The  block  matrixes  BQ')  and  vectors  G(J)  at  node  points 

0  <7  <  Af  are  not  determined  due  to  the  variables  at  expanded 

points  (with  superscript  0)  being  unknown.  As  a  result,  the 
tridiagonal  matrix  algebraic  Eq.  (39)  cannot  be  solved  directly. 
But  if  these  unknown  values  at  expanded  points  are  guessed  as 
trial  values,  then  Eq.  (39)  can  be  solved.  For  a  tridiagonal 
matrix  algebraic  equation,  a  compact  forward  and  backward 
substitution  algorithm  can  solve  it  quickly.  This  algorithm,  at 
first,  decomposes  the  tridiagonal  matrix  into  a  LU  form.  Then 
in  the  forward  substituting,  an  intermediate  vector  will  be 
acquired,  and  in  the  backward  substituting,  the  variable  vector 
can  be  obtained.  In  all  the  process  of  LU  decomposition  and 
forward  and  backward  substitution,  the  Gauss-Jordan  full 
pivoting  algorithm  is  needed  to  inverse  the  block  matrixes. 
Then  the  block  matrix  in  the  lower  and  upper  matnx,  the 
intermediate  vector,  and,  at  last,  the  variable  vector  can  be 
solved.  These  solved  variables  are  based  on  guessed  values, 
they  are  not  the  real  solution  for  the  problem.  They  are  set  as 
new  guessed  values  in  the  block  matrixes  B(J)  and  GO)  and 
reproduce  the  whole  procedure.  The  iteration  will  not  stop 
until  the  setting  accuracy  is  reached,  which  makes  the  solu¬ 
tions  satisfy  the  whole  difference  equations  at  each  node  point. 

The  procedure  is  illustrated  as  follows.  When  the  values 
for  determining  the  B(j)  and  GO)  are  guessed,  then  the 
coefficient  matrix  in  Eq.  (39)  can  be  decomposed 

rB(0)  D(0)  X 
A(l)  B(l)  D(l) 


AO)  BO)  DO) 

A{N-l)  B{N-\)  D{N-l) 
Y  A{N)  B{N) 
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m 

A(l)  B(l) 


A(/)  B{J) 


A{N-l)  B{N-l) 

YiN)  AiN)  B{N) 


W(0)  ■ 

G(0)  1 

W(l) 

G(l) 

• 

WO) 

= 

GO) 

WiN  -  1) 

GiN-1) 

V. 

.  WiN)  . 

.  GiN)  . 

Through  the  shift  of  matrixes,  the  block  matrix  in  the  L  and 
U  matrixes  can  be  solved.  For  j  =  0,  1,  we  have 


B(0)=B(0), 


x  =  r'(o)x, 


5(o)  =  r'(o)D(o) 

(41) 


A(1)=A(1),  B(1)  =  B(1)-A(1)D(1), 

5(l)=r'(l)D(l)-B''(l)A(l)X  (42) 

As  the  same  way,  for  other  node  points  at/  =  2, . . . , 
N-l, 

A(J)  =  A{j),  B{j)  =  B[j)  -  A(j)D(j  -  1), 

can  be  solved.  At  the  last  node  pointy  =  N,  we  obtain 

Y{N)  =  Y,  AiN)  =  AiN)  -  YiN)DiN  -  2), 

BiN)  =  BiN)  -AiN)DiN  -  1)  (44) 

Then  Eq.  (41)  can  be  solved  through  the  solution  of  two 
sets  of  matrix  algebraic  equations. 

5(0) 

A(l)  5(1) 


A(/)  Bij) 


AiN-l)  BiN-l) 

YiN)  AiN)  BjN) 


(45) 


w 


and 

I  D(0) 
I 


X 

D(l) 


/  m 


1  DiN-l) 
I 


G(0)  1 

W(0)  ■ 

C(l) 

W(l) 

CO) 

■  = 

WO) 

C(/V-l) 

WiN  -  1) 

.  CiN)  . 

WiN) 

(46) 


w 


The  intermediate  vector  Wij),  y  =  0, . . . ,  ^  in  Eq.  (45)  can 
be  determined  by  forward  substituting 


W'(O)  =  S”'(0)G(0) 

W(l)  =  r'(l)G(l)  -  5"'(1)A(1)IV(0) 

WO')  =  B~'ij)Gij)  -  S"'0)A0)W0  -  1) 

WiN)  =  B~\N)GiN)  -  B~\N)YiN)WiN  -  2) 
-B~'iN)AiN)WiN-l) 

With  backward  substituting,  the  solution  vector  C(j),j  =  0, 
N  are  solved. 


(47) 

(48) 

(49) 

(50) 


C(A)  =  WiN) 

dN  -  1)  =  WiN  -  1)  -  DiN  -  l)CiN) 
Cij)  =  Wij)-Dij)Cij  +  l) 

C(0)  =  W(0)-G(0)C(1)-XC(2) 


(51) 

(52) 

(53) 

(54) 
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The  solved  C(J)  are  as  new  guessed  values  to  determine 
the  block  matrixes  B(j)  and  G(j),  j  =  -  I  and  to 

continually  carry  out  the  iteration,  until  the  approximated 
solutions  are  obtained. 


4.  Results  and  discussion 
4.1.  Convergence 

To  test  the  solving  technique,  we  chose  a  set  of  physical- 
chemical  parameters  shown  in  Table  1  that  are  common  used 
for  hydrogen  oxidation  in  a  fuel  cell  anode,  then  the  follow¬ 
ing  dimensionless  model  parameters  are  obtained 

Co  =0.5714.  a  =  0.5,  U  =  0, 

=  0.1406,  =  =  0.8277 

Different  currents  density  (0.2-1 .6  A  cm~^)  at  two-point 
boundary  that  make  ki  and  h  change  between  0.2191- 
1.7526  were  adopted.  A  Fortran  95  code  was  developed 
to  do  the  numerical  solution  on  a  micro  computer  (Dell, 
Dimension  4100,  Pentium  III,  260  MB  memory).  After  the 
number  of  node  points  was  determined  and  a  set  of  proper 
initial  values  were  guessed  for  each  variables  at  all  the  node 
points,  the  iterating  process  was  carried  out.  It  converged 
very  quickly  in  just  two  loops.  The  numerical  results  are 
illustrated  as  Figs.  2-5. 

It  was  found  that  a  large  number  of  node  points,  for 
example,  >1000,  is  required  to  converge  the  iteration.  We  set 
it  as  2000.  Setting  too  many  node  points  is  the  need  for  our 
assumption  in  the  Taylor  expansion  that  the  nonlinear 
difference  equation  at  node  point  j,  depends  only  on  vari¬ 
ables  at  the  three  points.  The  iteration  precision  was  chosen 
as  10”^.  The  zero  in  the  right-hand  side  of  the  sets  of 
nonlinear  algebraic  equations  in  the  form  of  Eq.  (26)  is 
approximated  by  this  value.  If  all  the  absolute  values  of 
computation  for  the  terms  in  the  left-hand  side  are  not 
>10““',  then  the  iteration  processes  stop.  For  the  case  of 
2000  node  points,  correspondingly,  there  are  6000  algebraic 


Table  1  j 

Common  physical-chemical  parameters  in  modeling  the  hydrogen 


oxidation  in  a  fuel  cell  anode 


Parameter 

Value 

a  (cm"*) 

250 

Cu^  (mol  cm"^) 

0.4  X  10"“ 

(mol  cm"’^) 

0.7  X  10"^  0.38  X  lO"'*  for  =  0.15 

Dh,  (erne's"*) 

2  X  10"“ 

/cell  (A  cm"*) 

0.2-1. 6 

k°(cms-') 

0.45 

1  (cm) 

5  X  10"“ 

r(K) 

353.15 

V 

0 

a 

0.5 

K  (0"*  cm"*) 

0.3  X  10"' 

a  (Q"*  cm"*) 

0.3  X  10"' 

5000  Y*-'" 


Fig.  2.  Distribution  of  electrochemical  reaction  rate  of  hydrogen  along  x  at 
different  ki  (current). 


equations  like  Eq.  (26)  except  those  at  the  two-point  bound¬ 
ary.  Like  all  trial  and  error  iteration  algorithm,  different 
initial  estimations  yield  different  results,  in  order  to  make 
the  solution  meet  the  physical  meanings  in  the  system, 
following  current  criteria  are  introduced  in  our  computing 
to  determine  the  initial  trial  solution  and  check  the  results 


ii  =  -ff 


ii  =  -<r 


l2  =  -K 


d$i 

dX 

dX 

d^2 


=  /cell 


x=o 


=  0 


x=/ 


—  /cell 


X=/ 


(55) 

(56) 

(57) 


Fig.  3.  Distribution  of  dimensionless  hydrogen  concentration  cj  along  x  at 
different  ki  (current). 
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0  0 


Fig.  4.  Distribution  of  ^limcnsionless  potential  in  matrix  along  x  at 
different  kx  (current). 


Fig.  5.  Distribution  of  dimensionless  potential  in  pore  solution  C3  along  x 
at  different  (current). 


ii  =  - 

d^2|  n 

^  lx=o 

(58) 

/cell  = 

J> 

(59) 

For  Eqs.  (55H58),  the  five  point  formula  of  numerical 
differential  are  adopted  to  calculate  the  potential  gradient. 
For  Eq.  (59),  complex  Simposon’s  method  is  used  to  calcu¬ 
late  the  integration  along  the  catalyst  layer.  All  the  calcula¬ 
tions  in  these  five  equations  are  based  on  the  data  obtained 
from  the  numerical  solution  of  the  model  equations. 

4.2.  Characteristics  of  electrochemical 
reaction-diffusion 

The  characteristics  of  reaction-diffusion  problems  are 
usually  analyzed  by  the  Thiele  modulus  [10].  We  changed 
the  model  parameters  that  is  shown  in  Table  2  to  get  different 
Thiele  modulus  tP  in  a  range  of  0.008-10,  and  corresponding 
electrochemical  Thiele  modulu  and  ^2  in  n  range  of 
0.01-8.06  and  0.02-10.95,  respectively.  The  numerical 
solutions  converge  to  the  results  very  fast  too,  two  loops 
of  iteration  for  all  case.  The  results  are  showed  in  Figs.  6-9. 


0  0.2  0.4  0.6  0.8  1 

Fig.  6.  Dependence  of  reaction  rate  on  Thiele  modulus,  the  parameters 
corresponding  to  these  values  of  are  shown  in  Table  2. 

When  the  Thiele  modulus  is  low,  such  as  ^  =  0.008,  the 
resistance  of  diffusion  can  be  neglected,  the  distributions  of 
reaction  rate  and  hydrogen  concentration  are  nearly  uniform. 
When  the  Thiele  modulus  becomes  larger,  the  diffusion 


Parameters  for  studying  the  dependence  of  hydrogen  electro-oxidation  reaction-diffusion  on  Thiele  modulus  m  Figs.  6^ 


0 


^2 


Dh2  (cm^  s  ’) 


(cm  s"*) 


/  (cm) 


a  (fi”‘  cm“  ) 


K  (O"'  cm*"*) 


0.008 
0.8 
3 
5 
10 


0.012 

0.41 

1.41 

1.73 

8.06 


0.018 

1.85 

4.71 

6.42 

10.95 


2.0  X  10 

O.I  X  10 
0.08  X  10 
0.08  X  10 
0.08  X  10 


-4 

4 

•4 

-4 


0.45 

0.4  X  10^ 
0.8  X  10^ 
0.8  X  10^ 
0.8  X  10“ 


0.1  X  10 

0.25  X  10 
0.6  X  10 

1.0  X  10 

2.0  X  10 


-4 

,-4 


0.6  X  10 
0.3  X  10 


0.55  X 
0.1  X  10 


-1 

O-i 

I0-* 


-2 


0.3  X  10 
0.15  X  10 
0.268  X  10“ 
0.40  X  10’^ 
0.55  X  10“^ 


•  Non-changing  parameters:  C„.  =  0.4  x  lO-^*  mol  cm'’;  C„.  =  0.65  x  lO'^  mol  cm"’:  /„«  =  1-6  A  T  =  353.15  K;  t/  =  0;  a  -  0.5. 
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Fig.  7.  Dependence  of  dimensionless  hydrogen  concentration  on  Thiele 
modulus,  the  parameters  corresponding  to  these  values  of  <P  are  shown  in 
Table  2. 


Fig.  8.  Dependence  of  dimensionless  potential  in  matrix  on  Thiele 
modulus,  the  parameters  corresponding  to  these  values  of  are  shown  in 
Table  2. 


Fig.  9.  Dependence  of  dimensionless  potential  in  solution  on  Thiele 
modulus,  the  parameters  to  these  values  of  ^  are  shown  in  Table  2. 


resistance  increases,  as  a  result,  the  reaction  is  not  uniform 
along  the  catalyst  layer,  the  larger  the  Thiele  modulus  is,  the 
more  non-uniform  the  reaction  rate  is.  When  the  Thiele 
modulus  is  equal  to  10,  the  reaction  rate  and  hydrogen 
concentration  decreases  rapidly  in  the  region  near  the  sur¬ 
face  of  porous  catalyst  layer.  It  should  be  pointed  out  that  in 
high  Thiele  modulus  region,  the  rate  of  mass  transfer  by 
diffusion  cannot  meet  the  rate  of  electrochemical  reaction, 
the  process  will  lose  its  stability  and  become  a  dissipative 
system.  As  the  same  case  in  the  stability  issues  in  chemical 
reaction  engineering,  our  numerical  computing  also  found 
sensitive  dependence  of  the  solution  results  on  some  para¬ 
meters  when  the  Thiele  modulus  is  high.  It  has  been  reported 
by  a  large  amount  of  literature  that  in  open  chemical  reaction 
system,  especially  coupled  nonlinear  chemical  reaction  and 
diffusion  occur  simultaneously,  many  complex  phenomena 
will  take  place,  such  as  multiple  steady  states,  unstable  states 
and  self-generated  sustained  oscillations.  Our  present  work 
just  makes  a  beginning  in  this  field,  further  attention  should 
be  paid  to  it  because  the  electrochemical  oxidation  reaction 
rate  of  hydrogen  is  very  fast  under  the  catalysis  of  platinum 
catalyst.  When  a  fuel  cell  stack  is  needed  to  scale  up,  it  is 
necessary  to  know  the  dissipative  behavior. 


5.  Conclusion 

A  numerical  solution  method  is  developed  for  coupled 
electrochemical  reaction-diffusion  equations.  Even  though 
the  convergence  of  this  technique  depends  mainly  on  the 
large  number  of  node  points,  computation  experiences  on  a 
microcomputer  with  double  precision  real  only  eight  word 
bites  shows  that  the  procedure  converges  very  fast.  After  the 
Taylor  series  is  expanded  and  the  corresponding  iterating 
procedure  is  constructed,  the  computation  process  is 
stable.  Arranging  the  elements  in  the  coefficient  matrix 
to  block  matrix  form,  the  whole  coefficient  matrix  is  easily 
decomposed  to  lower  and  upper  matrix,  and  the  compact 
forward  and  backward  substitution  algorithm  based  on  the 
shift  of  block  matrixes  with  Gauss-Jordan  full  pivoting 
method  can  perform  the  numerical  calculation  quickly. 
Local  convergence  depends  on  the  first  trial  solution, 
current  criteria  are  required  to  make  the  solution  converge 
to  the  correct  results.  It  is  suggested  by  the  model  solutions 
that  dissipative  behaviors  in  a  electrochemical  reaction- 
diffusion  system  might  occur  when  the  Thiele  modulus  is 
high. 
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SDiitter  deposition  has  been  investigated  as  a  tool  for  manufacturing  proton-exchange  membrane  fuel  cell  (PEMFC)  electrode 

r"d  s  “ " 

Tee  and  six  layer  Pt/NCI  membrane  electrode  assemblies  (MEAs)  resulted  in  Pt  activities  double  tot  ^  ^ 

,/hLed  bv  the  ink-based  baseline.  Decreasing  the  thickness  of  each  NCI  layer  increased  the  performance  of  the  s  x-layered 
pJlIlCl  MEA ^0^0.13^  at  0.6  V.  providing  an  activity  of  2650  A/g  at  0.6  V.  It  is  likely  that  by  further  decreasing 

the  ratio  of  NCI  to  Pt  in  these  electrodes,  Pt  activity,  and  PEMFC  electrode  performance  can  be  increased. 
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Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  gaining 
popularity  due  to  their  high  operating  efficiency  and  environmental 
friendliness.  One  of  the  barriers  to  commercialization  is  the  prohibi¬ 
tive  cost  of  this  technology.'  In  a  recent  solicitation,  the  U.S.  De¬ 
partment  of  Energy  set  long-term  goals  for  PEMFC  performance  in 
a  50  k\V  stack  that  included  operation  with  cathode  loadings  of  0.05 
mg/cm^  or  less  of  precious  metals.^ 

Typical  methods  of  manufacturing  membrane-electrode  assem¬ 
blies  (MEAs)  for  PEMFCs  involve  painting,  spraying,  or  printing  of 
catalyst  inks  that  contain  a  matrix  of  electrolyte  and  carbon- 
supported  catalyst,^'®  It  is  this  three-phase  interface  of  catalyst,  car¬ 
bon,  and  electrolyte  (typically  Nafion®)  that  allows  effective  gas  and 
water  diffusion  and  proton  transport  and  electron  transport  to  and 
from  the  catalyst  sites.  Refinements  of  this  process  have  involved 
optimizing  the  ratios  of  Pt.  C,  and  Nafion  present  in  this  three-phase 
interface.^’^  There  are  limitations  on  the  catalyst  activity  imposed  by 
the  particle  size  of  Pt  on  activated  carbon. 

As  alternatives,  electrodeposition  and  sputter  deposition  have 
been  used  to  manufacture  MEAs  of  loW  catalyst  loadings.  ‘  '  Both 
pulse  and  direct  current  (dc)  electrodeposition  have  been  used  to 
localize  a  thin  layer  of  Pt  near  the  surface  of  the  MEA,  *  resulting  m 
the  development  of  electrodes  on  the  order  of  0.05  mg  Pt/cm  . 
Sputter  deposition  is  widely  used  for  integrated  circuit  manufactur¬ 
ing  and  has  been  investigated  for  the  preparation  of  more  effeefive 
fuel  cell  electrodes  for  more  than  a  decade.  Srinivasan  et  al  ’ 
applied  a  50  nm  Pt-sputtered  film  to  the  front  surface  of  a  catalyzed 
gas  diffusion  layer  (GDL)  to  reduce  the  loading  tenfold  (4  mg/enj^ 
to  0.4  mg  Pt/cm^)  without  reduction  in  performance.  Hirano  et  al  “ 
sputter  deposited  platinum  on  uncatalyzed  GDLs  resulting  in  cell 
performances  at  loadings  of  O.IO  mg  Pt/cm^  equivalent  to  those  of 
standard  methods  at  loadings  of  0.40  mg  Pt/cm^.  Sputter  deposition 
has  been  used  to  reduce  the  amount  of  anode  catalyst  required  for 
direct  methanol  fuel  cell  (DMFC)  anodes  as  well,  Witham  et  al 
achieved  DMFC  anode  catalyst  activities  one  to  two  orders  of  mag¬ 
nitude  higher  than  those  of  conventional  ink-based  catalysts,  sug- 
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gesting  that  DMFC  anodes  could  be  manufactured  containing  less 
than  one-tenth  the  amounts  presently  used  (2.5-4  mg  Pt/cm  )  with¬ 
out  loss  in  performance. 

Cha  and  Lee'"  alternated  sputtering  microthm  (5  nm)  Pt  layers 
and  painting  layers  of  Nafion  and  carbon  ink  directly  onto  the  mem¬ 
brane.  By  reducing  the  amount  of  Pt  on  each  layer,  they  were  able  to 
achieve  cell  performances  at  extremely  low  loadings  v 
mg/cm^)  that  were  nearly  equivalent  to  that  of  higher  loadings  (—0.5 
mg/cm^)  under  similar  conditions.  Very  efficient  usage  of  Pt  is  dem¬ 
onstrated  in  this  method. 

The  goal  of  this  work  was  to  examine  the  sputter-deposition 
technique  as  a  means  to  improve  performance  and/or  reduce  the 
catalyst  loading  of  proton  exchange  membrane  (PEM)  fuel  cells. 
First  Pt  was  sputtered  on  the  different  substrates  that  comprise  indi¬ 
vidual  fuel  cells.  These  were  built  into  MEAs  and  compared  to 
MEAs  made  through  ink-based  methods  and  ink-based  MEAs  aug¬ 
mented  through  sputter  deposition.  Based  on  these  results,  GDLs 
and  membranes  were  then  subjected  to  the  layered  technique  first 
developed  by  Cha  and  Lee'"  with  the  goal  of  reducing  the  amount  of 
Pt  catalyst  used  and  increasing  fuel  cell  performance.  Based  upon 
the  findings,  improvements  were  made  on  Cha  and  Lee  s  technique 
allowing  for  the  manufacture  of  thinner,  more  effective  electrodes 
for  MEAs. 

Experimental 

Catalyst  mfa.— The  method  for  the  catalyst  ink  preparation  and 
MEA  fabrication  performed  in  this  project  have  been  described 
elsewhere.'*  The  following  catalyst  inks  were  prepared:  (i)  Nafion 
-1-  carbon  only  (Nafion-earbon  ink  or  NCI)  and  (it)  Nafion 
-f  20%  Pt  on  carbon. 

The  inks  were  prepared  for  Pt  by  adding  the  E-TEK  catalyst 
(20%  catalyst  on  XC-72  carbon)  to  a  solution  of  5  wt  %  Nafion 
(DuPont).  In  the  case  of  NCI,  XC-72  carbon  was  added  to  a  solution 
of  5  wt  %  Nafion. 

Ink-based  Decals  (Teflon,  10  cm^  three  ply)  were 

weighed  prior  to  application  of  the  catalyst  ink.  The  ink  was  drawn 
across  the  surface  of  the  decals  using  a  Meyer  rod.  Anode  and  cath¬ 
ode  target  loadings  for  the  baseline  ink-based  MEA  were  both  0.15 
mg  Pi/cm^  for  a  combined  MEA  loading  of  0.30  mg  Pt/cm  .  A  low 
Pt  loading  was  chosen  as  the  baseline  because  low  catalyst  loading 
is  a  goal  for  commercializing  PEMFCs.*  The  coated  decals  were 
dried  in  an  oven  at  lOS-C  under  ambient  pressure  for  10  min. 
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To  fonn  a  MEA  with  ink-coated  decals,  appropriate  decals  were 
placed  on  either  side  of  the  PEM  (Nafion  117,  protonated  fonn). 
This  assembly  was  hot-pressed  to  ensure  bonding.  It  was  then 
cooled  to  room  temperature,  before  the  decals  were  carefully  peeled 
from  the  assembly.  An  uncatalyzed  gas  diffusion  layer  (Toray)  was 
placed  on  either  side  of  the  MEA  to  complete  the  membrane- 
electrode  unit  (MEU). 

Plasma  treatment  and  Pt  sputter  deposition.— Nation  1 17  mem¬ 
branes,  uncatalyzed  GDLs  (Toray)  and  MEAs  were  chosen  as  sub¬ 
strates  for  sputter  deposition.  Plasma  modifications  and  sputter- 
deposition  augmentations/additions  were  both  completed  using  an 
Anatech  Hummer  10.2  sputter-coating  tool.  A  modified  sample  stage 
was  used  to  support  PEM,  MEA,  and  GDL  substrates  of  sizes  up  to 
6  X  6cm  while  masking  1.5  cm  about  the  membrane’s  pennieter. 
An  aluminum  target  was  used  for  ac  and  dc  plasma  modifications, 
while  the  Pt  (Anatech)  targets  as  well  as  a  carbon  evaporation  sys¬ 
tem  (Anatech)  were  used  for  sputter-deposition  augmentations/ 

additions.  ,  .  . 

All  PEMs,  MEAs,  and  GDLs  subject  to  sputter-deposition  were 
first  ac  plasma  cleaned  for  a  period  of  5  min  at  5  m.A  and  1 .2  kV  to 
remove  residual  buildup  from  the  target  as  well  as  roughen  the  sub¬ 
strate  surface.  All  treatments  were  completed  at  a  pressure  of  ~62 
mTorr  A  separate  vacuum  chamber  was  used  to  evacuate  each  sub¬ 
strate  to  a  pressure  of  ~45  mTorr,  before  it  was  placed  m  the 
sputter-coating  tool  to  minimize  contaminant  out-gassing  in  the 
deposition  chamber. 

Since  all  plasma  and  sputter-deposition  treatments  must  be  per¬ 
formed  under  vacuum,  the  stability  of  Nafion  117  under  vacuum  was 
evaluated.  Each  PEM  material  was  slowly  evacuated  to  a  pressure 
below  1  mTorr.  Nafion  117  remained  very  pliable  following  ^acu- 
ation  Mass  measurement  of  Nafion  before  and  after  seven  different 
evacuation  and  processing  sequences  indicated  that  the  rnembrane 
loses  approximately  2.3%  of  its  total  mass  as  a  result  of  dehydration 
during  evacuation.  In  each  case  however,  Nafion  quickly  (~30  mm) 
rehydrates  upon  exposure  to  ambient  air.  The  exposed  membranes 
were  then  built  into  MEAs  and  subjected  to  performance  testing. 
They  performed  analogously  to  nonevacuated  membranes,  indicat¬ 
ing  that  Nafion  was  a  viable  candidate  for  plasma  and  sputter- 
deposition  treatments.  The  performance  of  a  vacuum-treated  MEA 
also  showed  no  adverse  effects  when  tested.  .  j  , 

A  potential  of  1.8  kV  and  a  cuaent  of  8  mA  was  mamtamed  to 
control  the  deposition  rate  for  Pt.  A  SiO^  sample  was  sputter- 
deposited  in  situ  with  each  PEM,  MEA,  or  GDL.  The  resultant 
metal/SiO,  stack  was  subjected  to  cross-sectional  view  SEM  imag¬ 
ing  to  verify  the  sputter-deposited  film  thickness  and  top-view  SEM 
imaging  to  determine  surface  characteristics  of  the  film. 

Multilayered  MEAs  were  prepared  by  first  spraying  the  anode 
and  cathode  sides  of  PEM  with  NCI.  The  MEA  was  then  subject  to 
a  vacuum  of  30  mTorr  before  the  appropriate  catalyst  was  sputter 
deposited.  This  sequence  was  repeated  until  the  desired  number  of 
layers  was  achieved.  Multilayered  GDLs  were  prepared  by  an  iden¬ 
tical  process  except  that  only  one  side  of  the  GDL  was  treated. 

MEAs  containing  sputter-deposited  catalyst  layers  were  built  into 
MEUs  through  methods  similar  to  that  of  the  baseline.  In  the  case  of 
multilayered  GDLs,  appropriate  GDLs  were  placed  on  either  side  of 
a  blank  Nafion  117  membrane  and  hot-pressed  to  ensure  a  well- 
bonded  MEU. 

Cell  assembly  and  to/mg.-The  MEUs  were  placed  in  a  10  cm^ 
cell  assembly  and  incubated  for  4  to  8  h  at  ambient  pressure,  cell 
temperature  of  WC.  stoichiometric  ratio  ([actual  flowj/ 
fstoichiometric  flow]  required  for  a  1.0  A/cm^  current)  of  1.5  at  the 
anode  and  2.0  at  the  cathode.  Fuel  cell  performance  curves  were 
obtained  under  the  conditions  set  out  in  Table  1. 


Table  I.  Fuel  cell  test  conditions. 


Pressure 

Cell  temperature 
Stiochiomctric  ratio  (at  I  A/cm") 

Feedstrcams 

Humidification 


I  atm 
70°C 

1.5  Hydrogen 
2.0  Air 

Anode:  hydrogen 
Cathode:  air 

Complete  humidification  of  anode 
and  cathode  gas  streams  for  all  trials, 


plasma  and  dc  plasma  substantially  roughen  the  surface  of  the  PEM. 
Analogous  roughening  by  ac  plasma  or  dc  plasma  did  not  signifi¬ 
cantly  increase  the  surface  of  MEAs  or  GDLs,  since  the  surfaces  of 
these  samples  are  much  rougher  in  their  original  state.  The  SEM 
images  in  Fig.  2  document  the  influence  of  ac  and  dc  plasma  modi¬ 
fication  of  MEAs.  The  high  and  medium  magnification  SEM  images 
both  suggest  ac  and  dc  plasma  minimally  alter  the  surface  or 
MEA.  However,  the  low  magnification  SEM  images  indicate  that  dc 
plasma  treatments  create  cracks  in  the  surface  layer  of  the  MEAs. 
Therefore,  ac  plasma  treatments  were  used  throughout  the  remainder 
of  this  program.  It  was  hoped  that  roughening  the  surface  of  the 
PEM  or  MEA  would  enhance  fuel  cell  performance  by 
the  number  of  active  catalytic  sites.  However,  performance  of  MEUs 
made  from  roughened  PEMs  and  MEAs  showed  no  improvement 
over  a  baseline  MEU. 

Determination  of  sputter-deposition  rates  and  catalyst 
loading.— The  top  view  SEM  images  in  Fig.  3  were  analyzed  to 
determine  the  surface  coverage  of  the  sputter-deposited  film  (Pt 
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Results 

AC  and  dc  plasma  roughening. — The  AFM  images  in  Fig.  1 
show  the  effects  of  an  ac  plasma  treatment  on  Nafion.  Both  ac 


eure  1.  AFM  images:  (a)  as-received  Nafion  117,  (b)  ac  plasma  rough¬ 
ed  Nafion  117  (potential,  1350  V;  current,  II  mA;  duration.  10  mm), 
ale:  1.0  jim  division. 
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As-received  MEA 


Hieh  Magnification  Medium  Magnification 


DC-plasma 
Roughened  MEA 
Target:  Al,  Potential: 


Low  Magnification 


~  65%).  Pt  did  not  form  continuous  films  on  the  SiOi  subshate, 
but  rather  agglomerated.  This  is  consistent  with  the  literature.  The 
surface  coverage  was  used  in  conjunction  with  the  bulk  density 
(Pt  =  21 .4  g/mL)  and  the  film  thicknesses  from  the  cross-sectional 
SEM  images  to  calculate  the  subsequent  Pt  loadings.  Based  on  this 
method,  the  Pt  sputter-deposition  rate  was  found  to  be  5.6  p.g 
Pt/cm^/min.  Pt  loading  analyses  of  all  sputtered  PEMs,  GDIs,  and 
MEAs  from  Adirondack  Environmental  Services  on  yielded  Pt  load¬ 
ings  that  were  generally  consistent  with  the  deposition  rate  deter¬ 
mined  from  cross-sectional  and  top  view  SEM  images. 

This  sputter  deposition  rate  (ixg/cmvmin)  for  Pt  was  found  to  be 
constant  with  respect  to  time  for  durations  of  5  min  or  more.  The 
cross-sectional  view  in  Fig.  3  shows  that  the  thickness  of  30  and  45 
min  sputter-deposited  Pt  is  roughly  two  and  three  times  the  thick¬ 
ness  of  a  15  min  deposition  of  Pt,  respectively.  Semiquantitative 
analyses  via  energy  dispersive  X-ray  spectroscopy  (EDXS)  and  Ru¬ 
therford  backscattering  spectrometry  (RBS)  confirm  these  Pt  depo¬ 
sition  rates.  These  methods  in  combination  with  the  analysis  pro¬ 
vided  by  Adirondak  Environmental  Services  were  also  used  to 
determine  loadings  of  individual  MEUs  composed  of  Pt  sputter- 
deposited  layers  of  less  than  5  min  duration. 

Sputter-deposition  ausmentation  and  addition.  The  sputter 
deposition  treatments  performed  for  this  experiment  are  listed  in 

Table  n.  ,  r  a 

Pt  was  sputter-deposited  on  GDLs,  PEMs,  and  MEAs.  figure  4 

compares  the  cell  performance  of  MEUs  made  from  GDLs  upon 
which  15,  30,  45,  and  90  min  of  Pt  was  sputter  deposited  as  a  single 
layer; in  all  cases,  tlie  equivalent  amounts  of  Pt  were  added  to  the 


anode  and  cathode.  The  MEU  comprised  of  30  min  of  sputter- 
deposited  Pt  is  closest  to  the  baseline  MEU  in  loading  (0.168  and 
0.150  me  Pt/cm^,  respectively).  Under  conditions  of  hydrogen/ah 
(anode/cathode)  feed  at  70  C  and  1  atm.  15  min  Pt  (0.084  mg/cm  ) 
deposited  on  the  anode  and  cathode  GDL,  0.147  A/cm'  w'as  ob¬ 
served  at  0.6  V  compared  to  0.276  A/cnr  for  the  baseline.  Doubling 
the  sputter-deposition  time  to  30  min  resulted  in  a  cell  performance 
comparable  to  the  baseline  MEA.  Figure  4  shows  that  further 
amounts  of  sputter-deposited  Pi  showed  no  appreciable  increase  in 

cell  performance.  m 

The  goal  of  the  sputter  depositing  an  additional  layer  of  Pt  on  the 
surface  of  the  anode  and  cathode  of  the  baseline  MEA  is  to  increase 
the  performance  of  the  electrodes  (primarily  the  cathode).  Oxygen 
kinetics  (at  high  potentials)  and  oxygen  diffusion  (low  potentials) 
are  two  factors  that  limit  the  performance  of  the  fuel  cell.  It  was 
hoped  that  the  application  of  this  additional  Pt  layer  would  increase 
the  rate  of  the  oxygen  reduction  (by  having  more  available  sites)  and 
that  by  having  a  layer  of  Pt  on  the  surface  of  the  MEA,  that  there 
would  be  less  diffusional  resistances  as  oxygen  would  not  have  to 
difhise  far  to  get  to  a  reaction  site.  Figure  5  compares  the  cell  per¬ 
formance  of  a  baseline  MEA  to  those  onto  which  Pt  has  been  sputter 
deposited.  Under  conditions  of  hydrogea'air  feed  at  70  C  and  1  atm, 
15  min  of  Pt  (0.84  mg/cm^)  deposited  on  the  anode  and  cathode  of 
a  baseline  MEA  provides  an  increase  in  perfomiance  at  high  volt- 
a^’es  (from  0.025  to  0.072  mg/cm^  at  0.8  V).  However,  the  added 
sputter-deposited  layer  of  Pt  clearly  reduced  the  limit  at  which  the 
gases  (oxygen,  hydrogen,  water  vapor)  could  diffuse  through  the 
electrodes.  This  elTect  resulted  in  a  reduction  in  performance  of  52% 


Figure  3.  Top  and  cross-sectional  SEM  images  (100,000  times)  of  sputtcr-deposited  Pt  on 
min  deposition  is  200  nm). 


Si/Si02  substrates.  Scale;  300  nm  (such 


that  the  thickness  of  the  45 


Table  II.  Sputter-deposition  treatments.  All  Pt  depositions  were 
performed  identically  on  the  anode  and  cathode.  In  the  case  of  a 
sputter-deposition  on  a  GDL  substrate,  identical  depositions  were 
performed  on  two  GDLs  that  would  then  form  the  anode  and 
cathode. 


MEA  (baseline) 

MEA  (baseline) 

Nafion  117 

Nafion  117 

Nafion  U7 

Nafion  117 

GDL 

GDL 

GDL 

GDL 

GDL 

GDL 

Nafion  117 
Nafion  117 
Nafion  117 
Nafion  117 
Nafion  117 
Nafion  117 
Nafion  117 


15  min  sputter-deposited  Pt 
30  min  sputter-deposited  Pi 
15  min  sputter-deposited  Pt 
30  min  sputter-deposited  Pt 
15  min  sputter-deposited  Pt  +  NCI 
30  min  sputter-deposited  Pt  +  NCI 
15  min  sputter-deposited  Pt 
30  min  sputter-deposited  Pt 
45  min  sputter-deposited  Pt 
60  min  sputter-deposited  Pi 
90  min  sputter-deposited  Pt 
NCI  +  (5  min  sputter-deposited  Pt  +  NCI)  X  3 
NCI  +  (15  min  sputter-deposited  Pt  +  NCI)  X  3 
NCI  +  (5  min  sputter-deposited  Pt  +  NCI)  X  3 
NCI  +  (2.5  min  sputter-deposited  Pt  +  NCI)  X  3 

NCI  +  (l.O  min  sputter-deposited  Pt  +  NCI)  X  3 

NCI  +  (0.5  min  sputter-deposited  Pt  +  NCI)  X  3 

NCI  +  (0.5  min  sputter-deposited  Pt  +  NCI)  X  6 

NCI  +  (0.5  min  sputter-depositedPt 
+  diluted  NCI)  X  6 


Substrate,  Addition 

- BASELINE 

-»-GOL  ISmInPtSD 
-^GOL,  30  min  Pt  SO 
-*-GDU45ininPtSD 

GOU  90  min  Pt  SO 


at  0.4  V  (from  0.56  A/cm'  to  0.27  cm')  for  the  case  where  a  15  min 
(0  084  nu’/cm')  Pt  layer  is  added  to  the  anode  and  cathode,  uas 
diffusion  is  all  but  stopped  in  the  case  where  30  min  of  Pt depos¬ 
ited  on  top  of  the  baseline  MEA,  and  the  result  is  negligible  cell 
performance.  This  lack  of  perfonnance  is  accompanied  by  an  in¬ 
crease  in  cell  resistance  from  25  to  200  mfl  as  shown  Table  11.  It 
is  likely  that  the  resistance  is  ionic  in  nature,  as  the  coating 


Current  Density  (A/cm*) 

Figure  4.  Performance  comparison  for  MEUs  prepared  from  GDLs  with 
various  sputter-deposited  Pt  loadings.  P  =  1  atm.  T  -  70"C. 


surface  of  the  MEA  fills  the  pores  necessao'  for  the  transport  of 
water  through  to  the  membrane  resulting  in  a  reduction  m  proton 

The  spiUter  deposition  of  Pt  onto  a  Nafion  117  substrate  jdelded 
an  MEU  that  showed  very  poor  P^tfojnwnce  as  shown  mPig^. 
Under  conditions  of  hydrogen/air  (anode, 'cathode)  feed  at  70  C  and 
1  atm,  15  and  30  min  of  sputter-deposited  Pt  deposued  on  ‘he  “nod 
and  cathode  side  of  Nafion  1 17  resulted  in  current  densities  of  0.044 
and  0.065  A/cm'  at  0.6  V,  respectively.  The  resistance 
open-circuit  voltage  was  very  high  compared  to  the  ‘’“^eline  (250  ^ 
25  mn).  In  contrast  to  the  sputter-deposition  augmented  MEAs.  it  is 
believed  that  the  high  resistance  is  the  result  of  poor 
duction  between  the  catalytic  Pt  sites  and  the  cathode.  As 
Fig.  3.  Pt  is  not  deposited  as  a  continuous  film,  but  an  agglomeration 
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Tabic  III.  Cell  resistance  of  various  MEUs  at  open-circuit  volt¬ 
age.  Resistances  were  measured  by  a  micrometer  set  at  1  kHz  ac 
voltage. 


Substrate 

Baseline 

Baseline 

Baseline 

Nafion  membrane 
Nafion  membrane 
Nafion  membrane 
Nafion  membrane 
Nafion  membrane 
GDL 
GDL 
Baseline 


Resistance 

Addition  ®  (mfl) 


15  min  Pt  SD 
30  min  Pt  SD 
15  min  Pt  SD 
15  min  SD  +  NCI 
30  min  Pt  SD 
30  min  Pt  SD  +  NCI 
NCI  +  3  X  (5  min  Pt  SD  +  NCI) 
15  min  Pt  SD 
30  min  Pt  SD 
Ru  filter  (ink-based) 


25 

27 

200 

250 

40 

83 

54 

27 

21 

23 

31 


Figure  6.  Performance  comparison  of  MEUs  prepared  from  PEMs  and 
GDLs  with  one  or  several  layers  of  sputter-deposited  Pt.  P  -  atm, 

=  70'’C. 


of  Pt  islands  with  a  columnar  microstructure.  The  addition  of  a  layer 
of  NCI  (50%  Nafion  solids  by  weight)  over  the  layer  of  sputter- 
deposited  Pt  decreases  cell  resistance  and  increases  perforrnance. 
Islands  of  sputter-deposited  Pt  are  connected  by  the  carbon  in  the 
NCI,  activating  the  sites  by  providing  a  pathway 
flow  fd  and  from  the  sites.  In  the  case  of  a  15  ram  Pt  deposited  on 


Figure  7.  Representation  of  a  single  layer  of  a  high  (left)  and 

nrsmitter-deoosited  Pt  between  layers  of  Nafion-carbon  ink  (NCI). 


Nafion,  Table  III  shows  the  drop  in  resistance  from  250  to  40  mil 
with  the  addition  of  the  NCI  layer,  while  Fig.  6  shows  ‘hat  there  is 
some  improvement  in  cell  performance  at  lower  voltages  (0  1 37  w. 
0.074  A/cm^  at  0.4  V).  The  Pt  islands  not  in  contact  with  the  GDL 
provide  no  conduit  for  the  electrons  and  are  thus  deactivated. 

For  an  equal  amount  of  sputter-deposited  Pt,  performance  of  the 
resultant  MEU  was  better  when  applied  to  the  GDL  as  opposed  to 
the  MEA  or  PEM.  The  sputter  deposition  of  Pt  directly  ‘h« 
brane  showed  far  less  performance  (0.044  vs.  0.138  A/cm  at  0.6  ) 
than  when  deposited  on  the  GDL.  The  rough  surface  of  the  GDL 
allows  for  the  generation  of  a  greater  Pt  active  area.  Also  the  high 
porosity  of  the  GDL  allows  gas  diffusion  to  and  from  the  electrodes 
even  after  Pt  deposition.  It  is  hypothesized  that  the  deposition  of  a 
single  Pt  layer  creates  regions  of  active  and  inactive  Pi  as  shown  m 
the  schematic  in  Fig.  7.  Pt  is  only  active  as  a  catalyst  when  it  is  in 
contact  with  the  electrolyte  and  a  conductive  support.  Pt  not  m  con¬ 
tact  with  Nafion  is  inactive  because  proton  transport  is  not  possible. 
The  inherent  roughness  of  the  GDL  increases  this  area  of  active  Ft 
compared  to  the  relative  smoothness  of  the  MEA  and  bare  mem¬ 
brane.  Even  on  the  GDL  the  deposition  of  Pt  beyond  a  given  amount 
(30  min  or  0.168  mg/cm^)  provides  no  added  perforniance,  oeca^e 
the  added  Pt  serves  only  to  increase  the  layer  of  inactive  Pt  de- 

*'"a  multilayer  electrode  technique''*  was  used  to  increase  the  re¬ 
gions  of  active  Pt  area  by  increasing  the  number  ot  layers  of  sputter- 
deposited  Pt.  MEAs  were  prepared  from  PEMs  with  multiple  layers 
of  sputter-deposited  Pi  and  spray-deposited  NCI  as  shown  in  Fig.  s. 
The  anode  and  cathode  were  built  in  an  identical  manner.  Figure  y 
compares  the  performance  of  cells  for  which  the  anode  and  cathode 
is  composed  of  three  sputter-deposited  Pt 
NCI.  Individual  Pt  layer  thicknesses  of  15,  5,  2.5,  1.0. 
resulted  in  electrode  loadings  of  0.2089,  0.07,  0.0795,  0.0407,  an 
0.027  mg  Pt/cm",  respectively.  As  the  thickness  of  each  layer  an 
resultant  Pt  loading  decreases,  the  cell  performance 
ofont  Thic  vprifif**;  the  theorv  that  the  fraction  of  inactive  Pt  shown 


Figure  8,  Representation  of  a  three  layer  Pt-NCI  MEA. 
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Fieurc  9.  Performance  comparison  between  the  baseline  and  three-layer 
MEAs  using  various  Pt  loadings.  =  I  atm,  7  =  70«C  Individual  Pt  layer 
thicknesses  of  15.  5,  2.5,  1.0,  and  0.5  min  resulted  in  electrode  loadings  of 
0.2089,  0.07,  0.0795,  0.0407,  0.027  mg  Pt/cm^  respectively. 


Figure  10.  Comparison  Pt  Activity  between  the  baseline  and  various  three- 
layer  anode  and  cathode  MEAs.  P  =  I  atm,  T  -  70‘’C. 


in  Fig.  7  is  eliminated  as  the  thickness  of  each  Pt  layer  is  reduced. 
The  amount  of  active  Pt  in  the  three-layer  electrode  did  not  change 
with  the  thickness  of  each  sputter-deposited  layer.  Furthermore,  all 
three-layer  MEAs  outperformed  the  single  Pt  layer  MEA 
(0.091-0.140  V5.  0.044  A/cm^  at  0.6  V),  suggesting  that  using  mul¬ 
tiple  layers  increased  the  active  Pt  area  available  for  reaction. 

Figure  10  shows  that  fuel  cell  activity  continuously  increases  as 
the  amount  of  sputter-deposited  Pt  in  each  of  the  multiple  i^yere  is 
reduced.  At  0.6  V,  the  activity  of  the  0.5  min  Pt/layer  MEU  is  1835 
A/g  Pt  vs,  905  A/g  Pt  for  the  baseline  MEA  and  670  A/g  Pt  for  the 
15  min  Pt  on  the  GDL.  Even  higher  activity  should  be  attainable 
since  a  sputter-deposited  Pt  layer  only  a  few  monolayers  thick  could 
provide  an  equivalent  number  of  active  catalytic  sites  (a  0.5  mm 
sputter-deposited  layer  of  Pt  is  roughly  equivalent  to  six  monolay- 
ers). 

Multilayer  treatments  were  of  no  benefit  when  applied  to  the 
GDLs  as  shown  in  Fig.  1 1.  This  is  consistent  with  the  high  degree  of 
roughness  and  porosity  of  the  GDL.  A  high  surface  area  allovys 
successful  distribution  of  sputter-deposited  Pt  throughout  the  GpL, 
As  a  result,  the  performance  of  the  single  15  min  sputter-deposited 
layer  of  Pt  on  the  GDL  is  greater  than  an  equivalent 
deposited  directly  on  the  membrane  in  three  layers  (0.147  vs.  0.091 
A/cm^  at  0.6  V). 


Figure  U.  Perfonnance  comparison  of  layered  PEM  and  GDL  electrodes. 
P  =  I  atm,  r  =  70®C. 


To  determine  the  effect  of  layering  on  the  performance  of  the 
MEA  when  the  layers  are  applied  to  the  membrane,  three  different 
layered  catalyst  structures  were  fabricated:  (1)  an  anode  and  cathode 
consisting  of  NCI  +  (1  min  of  sputter-deposited  Pt  +  NCI)  X  3; 
(2)  an  anode  and  cathode  consisting  of  NCI  +  (0.5  min  of  sputter- 
deposited  Pt  +  NCI)  X  6;  and  (3)  an  anode  and  cathode  consisting 
of  NCI  +  (0.5  min  of  sputter-deposited  Pt  +  NCI  diluted  1:5  with 

isopropanol)  X  6.  u  .  x<ca  n 

MEA  2  contains  the  same  amount  of  Pt  as  MEA  1,  but  MEA  / 
contains  twice  the  amount  of  Pt-NCl  interfaces  as  a  result  of  dou- 
blino  the  number  of  Pt  layers  from  three  to  six.  Comparing  the 
perfonnances  of  the  single-layer  15  min  sputtered  MEA  shown  in 
Fig.  6  and  the  three-layer  sputtered  MEAs  shown  in  Fig.  9,  the 
active  area  and  hence  cell  performance  is  a  function  of 
of  sputtered  layers  and  not  the  amount  of  Pt.  Therefore,  MEA  2  has 
roughly  double  the  active  area  of  MEA  1.  MEA  3  contains  the  same 
number  of  Pt  layers  as  2,  but  the  NCI  used  in  3  is  diluted  in  order  to 
reduce  the  thickness  of  the  NCI  region  between  each  Pt  layer.  Figure 
12a  shows  a  cross-sectional  view  SEM  image  of  a  spray-deposited 
NCI  layer  on  an  SiOj  substrate.  The  spray  deposited  NCI  layer  is 
~12  am  thick.  Therefore  for  MEA  2,  the  anode  and  cathode  stacks 
are  roughly  --70  p-m  thick  (the  thickness  of  the  Pt  is  negligible). 
This  is  much  thicker  than  the  anode  and  cathode  of  the  baseline 
MEU,  which  are  typically  5-10  p.m  thick.  Like  the  layers  of  sputter- 
deposited  Pt,  it  is  not  necessary  for  the  NCI  layers  to  be  so  thick. 
Based  on  the  fact  that  a  layer  of  15  min  of  sputter-deposited  layer  is 
75  nm  thick  (see  Fig.  3),  that  the  0.5  min  sputter-deposited  layer  of 
Pt  is  a  fraction  of  the  thickness  of  the  15  min  layer,  and  that  the  NL 
layer  is  roughly  12  p.m  thick,  it  can  be  concluded  that  the  NCI 
comprises  at  least  99.9%  of  the  thickness  of  the  multilayer 

electrodes.  ,  ,  ,  . 

To  reduce  the  thickness  of  the  NCI  layers  and  thus  reduce  these 
resistance  losses,  the  NCI  was  diluted  with  isopropanol.  The  thick¬ 
ness  of  a  spray-deposited  layer  of  NCI  diluted  1:1  by  volume  with 
isopropanol,  shown  in  Fig.  12b.  is  roughly  ~4  p.m  thick,  compared 
to  the  12  p.m  thickness  of  the  pure  NCI.  To  achieve  a  multilayered 
electrode  of  less  than  10  iJtm  (MEA  3),  a  dilution  of  1:5  NCI  to 
isopropanol  was  used.  . 

Figure  13  shows  that  as  the  number  of  electrode  layers  u  m- 
creased  from  three  (MEA  1)  to  six  (MEA  2).  the  cell  performance 
increases  (0.105  to  0.132  aW  at  0.6  V).  However,  there  is  not  a 
doubling  of  performance  that  accompanies  the  doubling  of  the  ac¬ 
tive  area.  The  increased  diffusional,  ionic,  and  electronic  resistances 
caused  by  the  increased  thickness  of  the  electrodes  of  MEA  2  vs. 
MEA  1  are  likely  causes  for  the  limited  performance  increase  seen 
in  Fig.  13,  Versus  the  three-layer  electrode,  gases,  protons,  and  elec¬ 
trons  must  travel  further  in  the  six-layer  electrode  to  get  to  and  from 
all  available  Pt  sites. 
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a)  Pure  NCI 


b)  Diluted  NCI 


Fisure  12.  Cro.ss-scctional  view  SEM  images  of  coating  produced  by  (a) 
pure  NCI,  (b)  diluted  NCI  (1:1  by  volume  with  isopropanol). 


Figure  13.  Performance  comparison  of  MEAs  comprised  of  three  and  six 
layers  (anode  and  cathode)  using  pure  NCI.  P  =  1  atm,  T  =  70°C. 


As  shown  in  Fig.  13,  MCA  3  outperforms  MCA  2  0.170  ‘o  0^*32 
A/cm-  at  0  6  V.  This  value  is  closer  to  double  that  of  MEA  1, 
containing  three-layer  electrodes.  This  further  emphasizes  the  fact 
that  as  the  thickness  of  the  NCI  decreases,  performance  increases. 
The  performance  of  MEA  3  is  also  greater  than  the  MEU  containing 
15  mih  of  Pt  was  sputtered  onto  the  anode  and  cathode  ODL,  but 


Fioure  14.  Representation  of  a  continuous  tlirec-phase  interface  (Nafion, 
cartion.  catalyst)  prepared  by  applying  Nafion.  carbon,  and  Pi  simulta- 


nf'iTU'slv. 


less  than  the  MEU  containing  30  min  of  Pt  on  the  anode  and  cathode 
GDL  and  less  than  the  performance  of  the  baseline  " 

ing  for  each  electrode  for  the  MEA  3  was  found  to  be  -0.0324  mg 
Pt/cm-  anode  /cathode  (0.0648  mg  Pt/cm'  in  the  entire  MEA)  result¬ 
ing  in  an  activity  of  2650  A/g  Pt  at  0.6  V.  the  highest  value  of  all 
MEAs  tested  and  three  times  greater  than  the  baseline  (90:)  A/g). 

Conclusions 

1  Of  the  three  substrates  studied  (membrane,  GDL,  and  MEA), 
sputter-depositing  Pt  on  the  GDL  showed  the  best  performance 
filling  that  of  the  baseline  MEA  for  an  equivalent  amount  of  P  . 
However,  sputter-depositing  multiple  layers  of  Pt  on  the  OU 
showed  no  improvement  in  performance  over  an  equivalent  amoun 

of  Pt  sputtered  as  a  single  layer. 

2.  MEAs  built  from  multilayered  sputter-deposited  Pt  and  spray- 
deposited  NCI  demonstrated  improved  performance  over  single- 
layer,  sputter-deposited  MEAs  of  equivalent  or  greater  Pt  loadings 
By  reducing  the  amount  of  Pt  and  NCI  in  each  layer,  a  'ovel  o 
activity  higher  than  that  of  the  baseline  has  been  achieved.  This 
level  of  activity  is  maintained  as  the  number  of  layers  increase  from 
three  to  six,  provided  the  NCI  layers  are  sufficiently  thim  The  prepa¬ 
ration  of  MEAs  with  extremely  high  activity  suggests  future  study. 
This  is  hypothesized  to  be  the  result  of  increasing  the  three  phases  ot 
Pt  Nafion,  and  carbon  necessary  for  an  active  catalyst. 

'  3.  The  optimal  performance  achieved  from  the  MEA  containing 
six-layer  Pt  +  dilute  NCI  anode  and  cathode  was  0.17  A/cm  itO-o 
V.  This  is  less  than  the  performance  achieved  by  Cha  and  Lee 
(roughly  0.32  A/cm^  at  0.6  V).  However,  all  fuel  cells  in  this  expen- 
ment  were  fed  Hj/air  at  1  atm  and  70°C,  compared  to  an  H,  /Oj  feed 
at  1  atm  and  60°C  used  by  Cha  and  Lee.  u-  j  . 

4.  Catalyst  activities  of  greater  than  2650  A/g  were  achieved  at 
0.6  V  and  5500  A/g  at  0.4  V  from  the  MEA  containing  six-layer 

Pt  -I-  dilute  NCI  anode  and  cathode.  . 

5.  Further  research  can  be  done  to  optimize  this  three-phase  in¬ 
terface  area  and  eliminate  the  unused  portion  of  the  electrode.  Even 
diluting  the  NCI  between  each  sputter-deposited  layer  5:1  with  iso- 
propanol  results  in  a  catalyst  electrode  in  which  NCI  accounts  for 
99%  of  the  width.  However,  to  generate  an  electrode  with  even  more 
than  six  Pt  4-  NCI  layers  is  neither  the  most  economical  (due  to  the 
time  required  to  generate  such  a  multilayered  MEA)  nor  most  effec¬ 
tive  approach.  The  generation  of  a  continuous  three-phase  interface 
is  the  ultimate  goal  of  this  method,  and  this  is  what  shou  d  be  pur¬ 
sued  using  the  method  of  sputter  deposition.  Simultaneously  sputter- 
depositing  Pt  and  spray-depositing  NCI  could  produce  a 
three-phase  interface  region  such  as  that  shown  m  Fig.  14  This 
would  result  in  an  extremely  thin  (-1  micrometer)  electrode  that  is 
conceptually  almost  identical  to  the  ink-based  catalyst  electrode 
used  in  the  baseline.  The  only  difference  is  that  the  particle  size  of 
the  sputter-deposited  Pt  would  be  much  smaller,  resulting  in  greater 
Pt  surface  area,  higher  Pt  activity,  higher  Pi  utilization,  and  hence  a 
superior  performing  electrode  compared  to  a  conventional  electrode 
prepared  from  colloidal  catalyst  inks.  Also  by  applying  this  Pt/C/ 
Nafion  electrode  in  a  single  application  to  the  proton-exchange 
membrane,  the  process  is  less  time  consuming  and,  thus,  more  eco¬ 
nomical. 


Journal  of  The  Electrochemical  Society,  149  (3)  A280-A287  (2002) 


A287 


Acknowledgments 

The  authors  acknowledge  the  financial  support  from  the  National 
Institute  of  Standards  and  Technology  under  cooperative  agreement 
no.  70NANB8H4039. 

University  of  South  Carolina  assisted  in  meeting  the  publication  costs  of 
this  article.  ^ 

References 

1.  Arthur  D.  Little  Inc..  Cost  Analysis  of  Fuel  Cell  Systems  for  Transportation:  Base^ 
line  System  Cost  Estimate,  Final  Report  to  Department  of  Energy  (2000). 

2.  Solicitation  for  Financial  Assistance  Applications  (SFAA)  no.  pE-RP04- 
0 1 AL67057  Research  and  Development  and  Analysis  for  Energy  Efficient  Tech¬ 
nologies  in  Transportation  and  Buildings  Applications,  Nov  21,  2W0. 

3.  E.  A.  Ticianclli,  C.  R.  Derouin.  and  S.  Srinivasan,  J.  Electroanai  Chem..  251.  275 

’  (iWs). 

4.  M.  S.  Wilson.  U.S.  Pat.  5.211,984  (1993).  ^  .  zr  / 

5  V  A  Paganin  E.  A.  Ticianelli.  and  E.  R.  Gonzalez,  tn  Proton  Conducting  Fuel 


Celb  A  S.  Gottesfcld,  G.  Halpert,  and  A.  Landgrebe.  Editors,  PV  95-23.  p.  102,  The 
Electrochemical  Society  Proceedings  Series,  Pennington,  NJ  (1995). 

6.  C.  K.  Wiiham,  W.  Chun,  T.  1.  Valdez,  and  S.  R.  Narayanan,  Electrochem.  Solid- 
State  Lett.,  491  {2000). 

7.  E.  J,  Taylor.  E.  B.  Anderson,  and  N.  R.  K.  Vilambi.y.  Electrochem.  Soc..  139,  L45 
(1992). 

8.  K.  H.  Choi.  H.  S.  Kim,  and  T.  H.  Lee,  J.  Power  Sources.  75,  230  (1998). 

9.  S,  Srinivasan,  D.  J,  Manko,  J.  Koch,  M.  A.  Enayetullah,  and  A.  J.  Appleby,  J. 
Power  Sources.  29,  367  (1990). 

10.  E.  A.  Ticianelli,  C.  R.  Derouin,  A.  Redondo,  and  S.  Srinivasan,  J.  Electrochem. 
Soc..  135.  2209  (1988). 

11.  S.  Mukeijee,  S.  Srinivasan,  and  A.  J.  Appleby,  Electrochim.  Acta.  38,  1661  (1993). 

12.  S.  Hirano,  J.  Kim,  and  S.  Srinivasan,  Electrochim.  Acta.  42,  1587  (1997), 

13!  C.  K.  Witham.  W.  Chun.  T.  I.  Valdez,  and  S.  R.  Narayanan,  Electrochem.  Solid- 
State  Lett..  3,  497  (2000). 

14.  S.  Y.  Cha  and  W.  M.  Lee,  J.  Electrochem.  Soc.,  146,  4055  (1999). 

15.  J.  A.  Thornton,  Deposition  Technologies  for  Films  and  Coatings,  p.  170,  Noyes 
Publications,  Park  Ridge,  NJ  (1982). 

16.  J.  A.  Poirier  and  G.  E.  Stoner.  J.  Electrochem.  Soc..  141,  425  (1994). 


